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I  j  The  objective  of  tble  progrsa  wee  to  conduct  en  Analytical 

I  iH  lnvee  get Ion  of  the  feeelblllty  of  ualng  e  spec lei  type 

ttUrf  of  transmission  theft  coupling  to  achieve  control  by  the 

— w  \j  i  designer  over  one  helicopter  drive  train  parent tar  -  torsional 

"  1  '  — -erifl/ieae  -  and  poeelbly  another  -  shaft  dashing. 


A  shaft  coupling  can  be  designed  to  Met  the  horsepower 
and  speed  requlrtMnt  of  modern  helicopter  drive  trains 
with  assentlally  sero  torsional  stiffness  et  the  design 
speed  and  power.  Two  different  types  of  aero  torsional 
stiffness  (ZT8)  couplings  were  designed  to  trensmlt  5270 
horsepower  at  13,820  rpn.  For  one  type  of  coupling, 
simple  design  equations  were  derived  to  else  a  coupling 
for  any  helicopter  application.  A  ZTS  coupling  in  a 
helicopter  drive  train  acts  as  a  torsional  vibration 
Isolator  and  can  provide  significant  damping  that  is  not 
otherwise  available  to  the  designer. 

The  technical  monitor  for  this  contract  was  Mr.  James  Gomes, 
Technology  Applications  Division. 
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An  Investigation  la  made  of  the  feasibility  of  designing  a  nonlinear  shaft 
coupling,  with  zero  or  very  low  torsional  stiffness  at  predetermined  oper¬ 
ating  conditions,  for  helicopter  drive  tr.ilns.  The  purpose  of  such  a  cou¬ 
pling  Is  to  Isolate  torsional  vibration  and  eliminate  resonarces. 


It  is  found  that  a  coupling  can  be  designed  to  have  zero  torsional  stiffness 
at  high  rotational  speeds  by  using  the  centrifugal  force  field  on  connecting 
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link*  from  the  driving  to  tho  drtm  flange  of  the  coupling.  Prevloua  te* 
■lfna  of  such  coupling*  for  lovi(>Hd  application*  km  boon  uaad  auffcaaa- 
fulljr  In  avlu  dloool  drive*. 

IM1|B  aquation*  for  too  typoo  of  s*ra-*tlffa**a  ooupllaga  fultabl*  for  high 
n>aad  turboahaft-powerod  haliooptara  a ra  developed. 

Coeputer  alnulatlona  of  tho  dynamic  performance  of  drive  traina,  with  and 
without  a  aero -tor  atonal  -otlffaaaa  ooupllag,  are  parfomad.  Both  open-loop 
and  oloaod-loop  (with  autonatlo  apaad  control)  ay  a  tana  are  ainulatad.  The 
result*  indioato  that  a  aero "toralonal -atiffnaaa  ooupllag  earn  prodoeo  three 
major  advantage*  In  a  helicopter  p'opulalon  ayatan: 

1.  Tho  ooupllag  can  ioolata  tha  angina  and  fual  oeatrol  fron  unavoid- 
abla  tor*  tonal  axel  tattoo  at  tha  rotor. 

3.  Tha  coupling  can  make  tha  achiavnaant  of  cloaad-loop  at  ability, 
through  proper  governor  daaign,  a  leaa  difficult  tank. 

3.  Tho  coupling  umkoa  it  poaalbla  to  incorporate  aignlficant  aaauata 
of  dangling  into  »  helicopter  drive  train,  thereby  lap  roving  tor- 
atonal  stability  and  further  euppreaslng  raaonanca. 
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This  report  presents  the  results  of  *  study  conducted  by  the  University  of 
Plorids  Engineering  end  Industrial  Experiment  Station  for  the  lustis  Dirac- 
tor  ate,  U.8.  Anqr  Air  Uobllity  Research  and  Bevel  opoent  Laboratory  under 
Contract  DAAJ02-72-C-0076,  Project  1 0192207 AA72. 

The  study  consisted  of  nstheauitical  analyses  of  sero-toralonal-etiffness 
shaft  coupling  designs  described  in  the  foreign  literature,  as  well  as  of 
now  design  concepts  generated  in  the  Department  of  Mechanical  Engineering 
at  the  University  of  Florida,  for  potential  application  in  helicopter  drive 
trains. 

The  principal  investigator  at  the  University  of  Florida  was  Dt*.  John  M. 
Vance.  Mr.  Roger  Brown  and  Mr.  Mark  Darlow  were  research  assistants. 

Part  I  of  the  report,  which  deals  with  design  of  the  couplings  and  analy¬ 
sis  of  their  effect  on  torsional  vibration,  was  primarily  authored  by 
Mr.  Brown.  Part  II,  which  describes  stability  analyses  of  closed-loop 
drive  trains  to  determine  the  effect  of  ZT8  couplings  on  torsional  stabil¬ 
ity,  was  primarily  authored  by  Mr.  Darlow.  The  Introduction,  Conclusions, 
and  Recoamsndations  were  written  by  Dr-  Vance,  who  also  edited  the  entire 
report. 

Dr.  Will  lam  H.  Boykin  of  the  Engineering  Science  and  Mechanics  Department 
at  the  University  of  Florida  served  as  an  Invaluable  consultant  on  questions 
of  stability  and  control  analysis. 
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INTRODUCTION 


TorsionnI  oscillation  is  a  significant  source  of  dynamic  loads  in  helicopter 
drive  trains,  which  must  transmit  ever-increasing  amounts  of  horsepower  with 
very  high  overnll  speed  reduction  ratios  over  significant  distances  between 
engines,  transmissions,  and  rotors. 

The  automatic  speed  governor  has  been  a  significant  development  for  heli¬ 
copter  flight  control,  :s  it  relieves  the  pilot  of  the  requirement  to  main¬ 
tain  constant  rotor  spied  during  changes  of  load  and  flight  conditions. 

A  special  problem,  which  has  been  encountered  in  almost  every  recent  heli¬ 
copter  development  urogram,  is  that  of  maintaining  torsional  stability  of 
the  engine  and  drive  train  while  at  the  same  time  providing  n  sufficiently 
rapid  response  to  demand  for  power  and  speed  changes.  The  helicopter  drive 
train  with  closed-loop  speed  control  is  a  system  which  is  difficult  to  opti¬ 
mize  without  inducing  torsional  instability  at  or  near  a  resonant  frequency. 

In  the  past,  whenever  the  problem  of  torsional  stability  has  been  encoun¬ 
tered.  a  typical  solution  has  been  to  reduce  the  gain  of  the  fuel  control, 
thereby  sacrificing  desirable  speed  control  properties  of  the  governor.1 
However,  there  are  other  critical  parameters  of  the  system,  such  as  drive 
train  torsional  stiffness  and  damping,  drive  train  inertia,  etc. ,  which 
might  be  used  as  design  variables  to  achieve  torsional  stability  without 
sacrificing  desirable  control  properties. 

Torsional  stiffness  has  always  been  a  difficult  design  parameter  to  control, 
since  it  is  effectively  reduced  by  the  square  of  the  very  high  gear  reduc¬ 
tion  ratio  of  the  main  transmission  (see  Appendix  II).  This  very  low  result¬ 
ing  torsional  stiffness  produces  a  correspondingly  low  natural  frequency  of 
the  drive  train  (typically  3-9  Hz),  which  makes  it  difficult  to  avoid  res¬ 
onant  peaks  in  the  closed-loop  (governed)  system. 

Shaft  damping  is  perhaps  even  more  difficult  to  control  in  a  drive  train, 
since  it  requires  large  amounts  of  vibratory  energy  dissipation  which  usu¬ 
ally  involves  friction,  wear,  or  temperature  dependence.  As  a  result, 
typical  values  of  shaft  damping  normally  found  in  helicopter  drive  trains 
amount  to  only  a  few  percent  of  the  critical  value,  producing  a  vibratory 
system  which  is  almost  completely  undamped  between  engine  and  rotor. 

The  major  objective  of  this  report  is  to  document  an  analytical  investiga¬ 
tion  of  the  feasibility  of  the  use  of  a  special  type  of  shaft  coupling  to 
achieve  control  by  the  designer  over  one  of  these  drive  train  parameters — 
torsional  stiffness — and  possibly  another — shaft  damping. 

Shaft  couplings  have  been  developed  independently  in  Kngland  and  in  Russia 
which  use  centrifugal  force  on  link  masses  (see  Figures  1  and  2)  to  produce 
torsional  stiffnesses  which  are  nonlinearly  dependent  on  speed  and  trans¬ 
mitted  torque.  In  particular,  these  couplings  have  been  designed  to  have 
zero  torsional  stiffness  at  predetermined  conditions  of  speed  and  torque. 


1 


2 

The  designs  developed  in  England  by  Chapman  are  being  used  to  eliminate 
torsional  resonances  In  marine  drives  with  diesel  propulsion.  Very  little 
information  is  available  on  the  Russian  designs,  other  than  Reference  3, 
and  the  specific  applications,  if  any,  are  unknown. 

The  marine  drive  applications  with  diesel  propulsion  involve  low  shaft 
speeds  (100-2000  rpm) ,  and  consequently  the  couplings  developed  for  such 
use  are  large  and  heavy  to  handle  the  high  torques  and  to  produce  the  re¬ 
quired  centrifugal  force. 

At  the  beginning  of  this  study,  it  was  postulated  that  a  zero-tors  Iona  1- 
stiffness  (ZTS)  coupling  could  be  developed  for  high-speed  application  in 
a  helicopter  drive,  and  could  eliminate  Instability  at  resonance  by  decou¬ 
pling  the  inertias  contributing  to  the  resonance.  denlly,  a  ZTS  coupling 
would  not  transmit  any  oscillatory  deviations  of  torque  because  of  its  low 
stiffness,  and  would  therefore  reduce  the  natural  frequency  of  the  system 
to  zero.  It  should  be  noted  that  this  is  opposite  to  recommendations'*  and 
efforts  in  the  past  toward  stiffening  the-  system  and  raising  the  natural  fre¬ 
quency  above  the  range  of  excitations,  a  procedure  which  has  been  found  to 
be  impractical  because  of  the  effect  of  the  gear  ratio  noted  above. 

In  practice,  a  ZTS  coupling  has  zero  stiffness  only  over  a  finite  range  of 
angular  twist,  and  large  excitations  drive  the  coupling  into  nonlinear  re¬ 
gions  with  some  effective  nonzero  stiffness.  Thus,  one  of  the  questions  to 
be  answered  in  Part  I  of  this  report  is  how  a  ZTS  coupling  actually  nffects 
the  natural  frequency  and  vibratory  response  of  a  helicopter  drive  train 
with  realistic  levels  of  excitation.  Another,  more  basic  question  addressed 
in  Part  I  is  whether  a  ZTS  coupling  can  be  designed  to  meet  modern  helicop¬ 
ter  requirements  within  the  constraints  of  size,  weight,  and  reliability. 
Finally,  results  of  stability  analyses  are  presented  in  Part  II  to  indicate 
the  predicted  effect  of  a  ZTS  coupling  on  closed-loop  torsional  stability. 

Before  proceeding  to  Parts  I  and  II  of  the  report  for  the  detailed  answers 
to  these  and  other  questions,  the  reader  may  find  it  helpful  to  consider 
how  zero  torsional  stiffness  can  be  achieved  in  a  shaft  coupling,  and  thus 
gain  some  insight  into  the  basic  concept. 

The  pinned  link  design  shown  in  Figure  1  lends  Itself  well  to  a  physical 
description  of  the  principle  of  operation.  Note  that  only  one  set  of  links 
is  shown.  A  symmetric  set  of  links  with  another  counterweight  would  be 
required  for  high-speed  balance. 

At  zero  speed  the  coupling  has  no  stiffness,  since  the  driving  flange  can 
rotate  freely  with  respect  to  the  driven  flange  with  no  torque  transfer 
(until  the  linkage  straightens  out). 

At  operating  speed,  centrifugal  force  holds  the  link  mass  (counterweight) 
out  at  some  equilibrium  value  of  a  which  is  determined  by  the  torque  being 
transferred.  Now  when  the  driving  flange  is  twisted  with  respect  to  the 
driven  flange,  two  effects  occur  simultaneously.  First,  the  link  mass 
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moves  inward  to  a  reduced  radius.  This  effect  results  in  a  reduced  cen¬ 
trifugal  force  and  therefore  tends  to  reduce  the  transmitted  torque. 

Second,  the  linkage  geometry  Is  changed  (the  angle  or  increases),  which  can 
cither  increase  or  decrease  the  transmitted  torque,  depending  on  the  link 
lengths. 

The  first  effect  produces  o  negative  torsional  stiffness  (positive  twist 
produces  a  negative  change  in  torque).  The  coupling  can  be  designed  so 
that  the  second  effect  produces  r,  positive  torsional  stiffness  which  is 
equal  and  opposite  to  the  negative  stiffness  from  the  first  effect.  The 
net  result  is  zero  torsional  stiffness  at  a  particular  design  speed  and 
angle  of  twist.  The  link  mass  is  sized  to  produce  the  uesired  torque  under 
the  design  conditions. 

In  the  design  of  Figure  2,  centrifugal  force  is  augmented  by  elasticity,  so 
the  coupling  stiffness  can  be  finite  at  zero  speed. 


PART  I: 


REDUCTION  OF  TORSIONAL  VIBRATION  THROUGH  THE  USE 
CF  ZERO-TORSIONAL-STIFFNESS  COUPLINGS 


PRELIMINARY  CONSIDERATIONS 

Rotating  machinery  is  inevitably  subjected  to  exciting  forces  which  set  up 
torsional  vibrations  throughout  the  system.  It  is  of  utmost  Importance 
that  the  exciting  frequency  be  different  from  the  critical  frequency  of  the 
system.  If  a  system  is  excited  at  or  near  its  critical  frequency,  the  sys¬ 
tem  energy  will  increase  with  each  successive  cycle.  The  limiting  factor 
for  the  energy  buildup  is  internal  or  external  damping  which  dissipates 
energy  in  the  form  of  heat. 

Excitation  at  the  critical  frequency  of  a  system  normally  designed  to  oper¬ 
ate  away  from  the  critical  can  come  about  through 

1.  Start-up  or  shutdown 

2.  Operation  at  other  than  design  speeds 

3.  Unexpected  exciting  forces 

4.  Unexpected  critical  frequencies. 

Figure  3a  is  an  example  of  an  open-loop  system  where  the  torque  input  is 
Independent  of  the  system  response.  Open-loop  systems  of  this  type  will 
not  become  torsionally  divergent  but  can  experience  severe  oscillation  when 
externally  excited  at,  or  near,  a  critical  frequency.  For  a  closed-loop 
system  (Figure  3b)  where  control  action  is  dependent  upon  system  response, 
instabilities  can  be  caused  by  internal  or  external  excitation.  Turbo- 
powered  machinery  such  as  helicopters  and  electrical  generators  are  closed- 
loop  systems  where  the  operating  speed  is  controlled  by  a  governor.  The 
speed  governor  causes  a  dependence  between  driving  torque  and  output  speed. 
Swlck  and  Akarvan,5  Poczkowski,4  Canale  and  MaCabe6  and  others  describe 
governor  systems  for  free  turbine  engines.  Part  I  considers  only  the 
open-loop  characteristics  of  torsional  systems,  such  as  free  turbine  heli¬ 
copter  systems,  which  are  subjected  to  external  excitation. 

Removing  critical  speeds  from  the  operating  region  in  the  open-loop  system 
will  generally  Improve  closed-loop  stability.  In  addition,  fatigue  fail¬ 
ure  due  to  cyclic  torque  variation  from  external  excitntionwill  be  lessened. 

Design  Considerations 

Critical  design  parameters  such  as  strength,  weight,  material,  and  size 
often  dictate  the  mass  distribution  of  a  system.  It  may  be  difficult  to 
change  the  mass  or  inertia  profile  of  a  system  without  adversely  affecting 
machine  performance.  Therefore,  the  designer  may  only  be  left  with  shaft 
stiffness  and  damping  as  free  parameters  with  which  to  ensure  safe  oper¬ 
ation  at  all  expected  operating  speeds. 


Figure  3.  Torsional  Systems. 
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With  proper  damping,  passage  through  critical  frequencies  during  start-up 
or  shutdown  is  usunlly  not  a  problem.  This  is  due  to  the  short  time  dura¬ 
tion  that  the  machine  is  driven  at  a  critical  speed.  However,  large 
amounts  of  damping  are  needed  when  there  is  continuous  excitation  at  or 
near  a  critical  frequency.  The  damping  will  prevent  the  system  from  becom¬ 
ing  divergent  and  magnifying  the  exciting  force  out  of  safe  bounds.  Unfor¬ 
tunately.  this  type  of  damping  is  generally  hard  to  achieve  in  mechanical 
drive  trains  in  addition,  a  highly  damped  system  will  transmit  more 
exciting  force  than  a  lesser  damped  system  when  operation  is  above  critical 
speed. 

If  the  exciting  frequency  is  known,  it  may  be  possible  to  design  a  shaft 
that  has  a  stiffness  which  gives  the  system  a  critical  frequency  much 
higher  or  lower  than  the  frequency  of  excitation.  However,  it  may  often 
be  difficult  to  design  a  shaft  with  the  desired  stiffness  and  still  have 
acceptable  strength,  weight,  and  deflection  characteristics. 

Nonlinear  Couplings 

A  low-stiffness  shaft  coupling  may  be  used  to  solve  the  problem  of  excita¬ 
tion  at  a  critical  frequency.  This  is  done  by  lowering  the  critical  fre¬ 
quencies  of  a  system  to  nearly  zero,  while  at  the  same  time  transmitting 
the  system  torque  requirements.  Figures  1  and  5  show  the  characteristics 
of  some  linear  and  nonlinear  couplings.*  it  is  apparent  from  Figure  1  that 
a  linear  torsional  spring  with  near  zero  stiffness  would  require  large 
amounts  of  twist  in  order  to  transmit  the  design  torque.  The  nonlinear 
spring,  however,  can  transmit  design  torque  with  zero  torsional  stiffness 
at  a  small  angle  of  twist.  One  of  the  objectives  of  this  report  is  to 
determine  if  such  a  nonlinear  "spring”  (coupling)  is  a  practical  possibil¬ 
ity  for  high  torques  and  speeds.  Present  designs  of  nonlinear  couplings 
which  have  zero  or  low  stiffness  for  certain  operating  conditions  are  in¬ 
vestigated  in  the  following  section. 

For  linear  systems  the  resonant  frequencies  are  inherent  to  the  system  and 
are  not  functions  of  the  exciting  forces.  The  resonant  frequencies  of  non¬ 
linear  systems  are  functions  of  the  amplitude  and  period  of  the  exciting 
force,  as  well  as  the  system  state  when  excitation  occurs.  As  can  be  seen 
in  Figure  5.  if  the  torque  is  varied  over  a  large  range,  the  nonlinear 
spring  stiffness  increases  on  either  side  of  the  low  stiffness  operating 
region.  In  order  to  obtain  an  "average"  stiffness  over  a  cycle,  the  time 
variation  of  the  torque  must  be  considered. 

The  Helicopter  Problem 


One  objective  of  this  report  is  to  develop  a  procedure  for  the  design  and 
analysis  of  a  zero-torsional -stiffness  coupling,  with  particular  emphasis 


*Torsicnal  stiffness  is  defined  as  the  ratio  of  the  change  in  coupling 
torque  to  the  change  in  angle  of  twist  across  the  coupling,  i.e. .  the 
slope  of  the  torque  vs.  twist  angle  curve. 


K 


0  - 

Twist  Angle 


Figure  4.  Coupling  Characteristics 


Figure  5.  Coupling  Stiffness. 
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on  helicopter  applications.  The  values  used  for  the  dynamic  characteris¬ 
tics  of  the  helicopter  drive  train  were  obtained  from  a  helicopter  sta¬ 
bility  analysis  by  Swick  and  Skarvan.5  Figure  6  shows  the  basic  helicop¬ 
ter  drive  systen. 

It  has  been  found  in  past  helicopter  designs  that  the  closed-loop  systen 
can  become  unstable  under  certain  flight  conditions.  When  this  happens, 
sel f-exci t at ioi.  from  a  snail  perturbation  can  quickly  cause  destructive 
vibrations.  It  is  hoped  that  by  Including  a  zero-torsional-stiffness  cou¬ 
pling  as  part  of  the  overall  systen,  all  torsional  Instabilities  can  be 
elinlnated.  The  problen  of  fatigue  failure  of  elastic  aenbers  can  also 
be  reduced  by  lowering  the  amplitude  and  frequency  ol'  self-excited  vibra¬ 
tions  as  well  as  the  anplitude  of  vibration  caused  b"  external  excitation. 
The  major  external  sources  of  torsional  excitation  :  e: 

1.  Pilot  excitation  of  control  paraneters 

a.  collective  pitch 

b.  engine  speed  setting 

2.  Aerodynamic  excitation  of  rotor  blades. 

As  helicopter  ground  speeds  Increase,  the  whip  effect  caused  by  aerodynamic 
excitation  of  the  rotor  blades  can  be  expected  to  increase.  This  is  due  to 
the  changing  relative  air  velocities  the  rotor  blade  sees  as  it  enters  and 
leaves  the  direction  of  flight.  Figure  7  shows  that  the  maximum  velocity 
the  rotor  blades  see  during  each  revolution  is 


and  the  minimum  velocity  is 


This  effect  combined  with  the  cylic  variation  of  the  rotor  blades  to 
achieve  forward  motion  causes  substantial  torque  variations  on  the  rotor 
and  drive  train. 

An  added  incentive  for  the  development  of  a  versatile  coupling  for  hell' 
copters  comes  from  the  nature  of  helicopter  design  and  development. 

Nearly  every  major  component  of  a  helicopter  drive  train  is  made  by  a  dif¬ 
ferent  company.  It  is  usually  the  responsibility  of  the  airframe  manufac- 
tuier  to  couple  all  the  drive  train  components  together.  Thus  it  is  very 
possible,  unless  there  is  exceptional  coordination  during  the  component 
design  phase,  that  the  assembled  system  will  have  critical  frequencies 
which  are  unacceptable. 
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Figure  e.  Closed-Loop  Helicopter  Drive  Train. 
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helicopter  velocity 


Rotor  Blade  Velocities. 


Tvc  ~'ncept«  of  zero-torslonal-atiffness  couplings  for  high-speed  helicop¬ 
ter  applications  are  analyzed  later  in  this  report.  In  addition,  a  Lagran- 
glan  formulation  for  a  dynamic  simulation  of  an  open-loop  helicopter  drive 
train  is  Included.  The  simulation  is  used  to  determine  the  effectiveness 
of  a  zero-torsional -stiff ness  coupling  in  reducing  torsional  vibrations. 

The  procedure  used  for  analyzing  torsional  vibrations  of  the  open-loop  sys¬ 
tem  readily  lends  Itself  to  the  inclusion  of  governor  dynamics  for  closed- 
loop  analysis. 

A  synthesized  design  of  a  zero-torslonal-st  iff  ness  coupling  for  helicopter 
applications,  and  the  results  and  discussion  of  the  drive  train  simulation 
study  using  the  synthesized  coupling  are  presented  in  subsequent  sections. 

SURVEY  OF  EXISTING  LOW -TORSIONAL  -STIFFNESS  COUPLINGS 

Low -torsional -stiffness  couplings  have  been  designed  to  place  the  critical 
frequency  of  a  system  below  the  excitation  frequencies.  The  couplings  can 
eaqiloy  rigid  and  elastic  linkages,  as  well  as  a  variety  of  torsional  and 
linear  springs.  Centrifugal  force  caused  by  the  rotating  links  plays  an 
important  role  in  determining  the  stiffness  characteristics  of  many  of  the 
couplings. 

The  geometry  of  the  coupling  (arrangement  of  components)  is  the  dominant 
factor  in  determining  the  torque  resulting  from  either  the  centrifugal 
force  of  the  links  or  the  force  due  to  the  elastic  deformation  of  the  links 
or  springs. 

This  section  describes  various  coupling  designs  which,  for  certain  operat¬ 
ing  conditions,  have  low  torsional  stiffness  and  the  capability  of  trans¬ 
mitting  high  torque. 

Elastic  Link  Designs 


The  elastic  materials  used  in  coupling  designs  range  from  rubber  to  steel. 
For  systesm  with  low  torque  requirements,  a  rubber  or  elastomeric  material 
may  be  used  successfully  to  achieve  low  stiffness.  For  the  more  common 
high  torque  systems,  steel  links  have  been  employed  to  achieve  low  tor¬ 
sional  stiffness. 

The  most  notable  of  the  high  torque  elastic  link  designs  is  the  half-moon 
(HU)  coupling  developed  by  C.  V.  Chapman.2  Figure  *  is  a  diagram  of  the 
HU  coupling  showing  the  elastic  (steel)  links.  Pins  X}  and  X2  are  equally 
spaced  about  0  and  connected  to  the  driving  flange.  Pins  Yj  and  Yj  are 
also  equally  spaced  about  0  but  are  con ec ted  to  the  driven  flange.  When 
one  flange  is  rotated  relative  to  the  other  about  point  0,  torque  can  be 
transmitted  across  the  coupling  through  the  mechanism  of  elastic  deforma¬ 
tion  of  the  steel  links.  The  zero-speed  torque-deflection  curve  of  the 
HM  cc  v  ling  is  shown  in  Figure  9.  At  other  than  zero  speed,  the  centrif¬ 
ugal  j.-ce  of  the  links  enables  additional  torque  to  be  transmitted  across 
the  ^.upllng.  The  additional  torque  is  proportional  to  the  square  of  the 
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coupling  speed.  Figure  8  shows  that  the  center  of  gravity  of  each  link 
aoves  toward  0  as  the  twist  of  the  coupling  Is  increased.  As  a  result  of 
the  path  of  the  center  of  gravity,  the  ajaount  of  additional  torque  trans¬ 
acted  across  the  coupling  resulting  from  centrifugal  force  is  reduced 
with  increasing  twist  angle.  If,  at  a  given  speed,  the  decrease  in  torque 
resulting  from  centrifugal  force  is  larger  in  magnitude  than  the  increase 
in  torque  resulting  from  elastic  link  deformation,  a  negative  torsional 
stiffness  will  result.  Chapman  overcomes  the  negative  stiffness  problem 
and  the  resulting  Instability  for  operation  up  to  a  maximum  design  speed, 
by  including  some  torsional  rigidity  in  the  bushings  at  pins  Xj ,  X2 ,  Yj  , 
and  Y2-  The  HM  coupling  has  been  quite  successful  in  low  speed/high  power 
applications.  Chapman  has  indicated  through  correspondence  that  he  has 
been  unable  to  design  a  steel  elastic  link  coupling  that  can  operate  at 
the  high  speeds  found  in  turbine  engine  helicopter  drive  trains. 

A  modified  arrangement  of  the  elastic  links  in  the  coupling  is  shown  in 
Figure  10,  where  the  pins  Xj ,  X2.  Y*  ,  and  Y2  are  not  necessarily  col  inear 
in  the  unstrained  position.  The  X  and  Y  pins  are  the  same  distance  from 
0,  as  opposed  to  Chapman's  staggered  pin  design.  The  next  section  presents 
an  analysis  technique  for  the  design  of  Figure  10  which  will  hopefully 
produce  elastic  link  couplings  with  improved  stability  characteristics  for 
high-speed  applications. 

Pinned  Link  Designs 


Two  pinned  link  coupling  designs  are  known  to  have  been  developed.  One  is 
the  toggle  link  (TL)  coupling  designed  by  C.  W.  Chapman.2  The  other  is 
a  design  of  Russian  origin.3  An  illustration  of  the  TL  coupling  is  shown 
in  Figure  11.  The  torsional  stiffness  of  this  coupling  is  effectively 
zero  at  zero  speed  but  increases  as  the  speed  increases.  Chapman  has  been 
able  to  produce  TL  couplings  with  torsional  stif fness-to-mean  torque  ratios 
of  4  or  less  for  low  speed/high  power  applications.  This  ratio  is  a  measure 
of  relative  stiffness  which  Chapman  uses  as  criteria  for  effectively  elim¬ 
inating  torsional  vibrations.  The  TL  coupling,  like  the  HM  coupling,  has 
some  torsional  rigidity  in  the  bushings  to  overcome  negative  stiffness  at 
operating  speed.  The  Russian  pinned  link  design,  shown  in  Figure  12,  con¬ 
sists  of  a  mass  connected  by  two  links,  one  to  the  driving  flange  and  the 
other  to  the  driven  flange.  The  stiffness  properties  of  this  coupling  are 
dependent  upon  the  mass,  weights  of  the  links,  speed  of  rotation,  link 
length,  and  pin  location  on  the  flanges.  Reference  3  gives  very  little 
analytical  background  which  is  of  much  use  in  designing  a  coupling  of  this 
type.  There  are  no  graphs  given  in  this  reference  and,  in  addition,  an 
equation  intended  to  represent  the  torque-deflection  characteristic  of  the 
coupling  is  dimensionally  Incorrect. 

Low  stiffness  couplings  up  to  now  have  been  designed  for  low-speed  oper¬ 
ation  (500-3000  rpm)  as  would  be  found  in  diesel  engine  systems.  For  tur¬ 
bine  engine  applications,  such  ns  found  in  many  helicopter  drive  trains, 
drive  shaft  speeds  of  10,000  to  20,000  rpm  are  not  uncommon. 
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A  new  pinned  link  design  is  developed  in  the  next  section,  "Design  Analysis 
for  licoptcr  Applications."  It  is  shown  that  a  simple  coupling  can  be 
synthesized  which  has  zero-torsional  stiffness  at  high  speeds  without  the 
aid  of  bushing  stiffness  or  other  added  torsional  stiffness  correctors 
which  are  used  to  avoid  negative  stiffness. 

DESIGN  AXALYSIS  FOR  HELICOPTER  APPLICATIONS 


Elastic  Link  Design  Analysis 


A  procedure  for  analyzing  the  elastic  link  coupling  shown  in  Figure  10  is 
developed  in  this  section.  The  parameters  involved  in  determining  the 
stiffness  properties  and  strength  characteristics  in  this  analysis  (see 
Figure  13)  are 

2  4 

Density  of  link  a.  lb-sec  in. 

Dimensions  of  link  u.b.d.t,  in. 

Initial  position  of  link  r.a  ,  rad. 

o 

Modulus  of  elasticity  of  link  E,  lb/in. 2 

Speed  of  rotation  uu,  rad/sec 

2 

Yield  strength  a  ,  lb  in. 

IQ8X 

Number  of  links  N  . 

The  location  of  the  center  of  mass  for  the  half-cylindrical  link  is 

2  2 
4  n  +  nb  +  b 

C  "  3r  a  +  b 

Due  to  the  symmetrical  distribution  of  pins  A  and  B  about  the  link  center 
line,  the  force  on  the  pins  due  to  angular  rotation,  u,  of  the  coupling 
for  a  link  of  mass  m  is 


_  cm2,  . .  .  .  . 

Fy  =  T  =  2  X  (C  +  r  COS<^0O>1 

(assuming  c  *  c(u ,F  ) 
c 

where 

1  2  2 
m  =  -  rr  d  o  (b  -  a) 


(1) 


(2) 
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Figure  13.  Elastic  Link  Coupling  Parameters. 
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The  force  in  the  x  direction  on  pins  A  and  B  is  a  function  of  the  deflec¬ 
tion,  Ax.  It  is  not  a  simple  relationship  and  requires  a  Theory  of  Elas¬ 
ticity  analysis  for  accurate  results.  An  expression  for  the  spring  con¬ 
stant,  Kp.  due  to  elastic  deformation  of  the  link  is  developed  from  the 
Theory  of  Elasticity  in  Appendix  I,  the  results  of  which  are  given  below. 
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For  a  rotation  of  radians  of  each  pin  the  change  in 
between  A  and  B  can  be  related  to  Ax  by 

0+8  0 
AX  =  2r[sin( - g-— *  ”  sin  “J"! 


relative  displacement 
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and.  thus,  the  force  directed  between  A  and  B  is 
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The  torque  about  0  resulting  from  forces  acting  at  A  is  equal  and  opposite 
to  the  torque  about  0  resulting  from  forces  acting  at  B.  Since  B  is  on 
the  driving  flange  and  A  is  on  the  driven  flange,  the  torque  transmitted 
by  the  coupling  is 

0 


0 


T  =  ~FBx  r  COS(* 


0  +  eo 

)  +  Fey  r  sin  ( - y-—) 


(6) 


Equations  1  and  5  arc  substituted  into  equation  6.  The  resulting  equation 
is  divided  by  d  and  the  right  side  is  multiplied  by  the  number  of  links,  S, 
giving 


T  N  E  r  ,,2 
_  =  - - - —  ( (a 

n(n  +  b  ) 


b2) 


.  2  J2.  ,  b. 

(a  +  b  )  log  — 1 
a 


•  [sin  (. 


Ntt  2,.2  2. 

-  sin  °  1  +  —r  c  x  (b  -  a  ) 
~T  4 


+  ab 


_4_  U> _ 

1 3n  a  + 


a  ) 


0  +  e 


+  r  cos(- 


■)]r  sin(- 


(7) 


It  is  necessary  to  determine  the  coupling  twist  angle,  *hich  produces 
the  maximum  allowable  stress  in  the  link.  At  zero  speed  the  stress  at 
point  D  in  Figure  13  will  be 


Mt  FAx 
CD  21  +  td 
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where 


•  •l-.,  »  -  V  *  I1'  *"d  ’  *  72 


giving 


(3a  +  2t) 

tZd 


(8) 


At  speeds  other  than  zero,  the  stress  at  points  A  and  B  oust  be  considered. 
However,  the  pins  at  these  points  can  be  made  large  to  account  for  stress 
without  affecting  the  stiffness  properties  of  the  »  >upling. 

Finding  a  coupling  with  zero  or  low  stiffness  which  can  transmit  the  speci¬ 
fied  torque  requirements  at  design  speed  is  a  formidable  problem.  Choosing 
»  material  such  as  steel  for  the  links  removes  E,  cx,  and  p  from  the  list 
of  free  design  parameters,  leaving 

0o  b 

r  d 

a 

Since  d  is  factored  out  of  the  right  side  of  equation  7,  the  ratio 

T 

“  =  d 

can  be  plotted  against  8  for  various  values  of  60,  r,  a,  and  b. 

The  determination  of  a  low  stiffness  coupling  requires  finding  a  *'<-9  curve, 
whose  slope,  (1/d)  (dT/d0) ,  is  zero  or  near  zero  but  never  negative.  This 
can  also  be  achieved  by  directly  plotting  the  slope 


versus  6  for  various  values  of  the  design  parameters.  Using  the  value  of 
D  where  low  stiffness  occurs,  and  the  amount  of  torque  which  the  coupling 
must  transmit  at  design  speed,  the  value  of  d  can  be  evaluated  from 


d  = 


T 

a 


Equation  8  must  be  applied  to  see  if  the  maximum  allowable  stress  In  the 
link  Is  exceeded  for  the  specified  set  of  design  parameters. 

A  slight  change  In  the  design  speed  or  design  torque  would  require  one  to 
start  the  search  technique  over  again.  The  practicality  of  designing  this 
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type  of  coupling  for  helicopter  use  by  trial  and  error,  where  the  design 
parameters  are  constantly  changing,  seems  very  low.  On  the  other  hand, 
if  design  charts  were  generated  showing  stiffness  propertie"  for  various 
speeds  and  torques,  this  type  of  coupling  might  prove  quite  useful. 

Under  a  set  of  simplifying  assumptions,  and  with  the  design  torque  and 
speed  known  and  invariant,  a  simplified  set  of  design  equations  can  be 
developed  which  renders  a  trial  and  error  solution  more  tractable,  as 
follows: 

The  assumptions  arc  (see  Figure  13) 

1.  9  =  0  at  design  torque  and  speed.  (All  of  the  torque  transmitted 
at  design  conditions  will  be  due  to  centrifugal  force.) 

2.  c  a  =  constant  =  0.636  for  small  9  (the  link  is  semicircular  and 
t  a  is  small). 

3.  Elastic  and  centrifugal  effects  are  independent. 

4.  Elastic  effects  are  approximated  by  the  bending  solution,  using 
the  method  of  Cnstigl inno. 17 

5.  The  pins  are  symmetrically  placed  at  equal  radii. 

6.  The  material  is  steel. 

The  design  speed  and  torque  are  calculated  from  the  requirements  of  the 
Swlck  and  Skarvan  model  as 


13.820  rpm  X  —  =  1447  — 

60  sec 

(9) 

(12) (550) 

~ X  5000  hp  =  22,850  in. -lb 
(1447 ) 

(10) 

Thus  the  torque  transmitted  per  link  in  a  two-link  coupling  is  11 ,425  in. -lb. 

The  weight  of  a  semicircular  steel  link  is  given  by  (see  Figure  14  for  link 
dimensions) 

W  =  0.89  abt  ♦  0.445  bt2  (11) 

The  centrifugal  load,  concentrated  at  the  eg  of  the  link,  is 

2  9 
Fc  =  raU)  (h  +  c)  ,  where  h  =  r  cos  (-^)  (12) 

or 

Fc  =  (h  ♦  .636a)  (4810  abt  ♦  1405  bt2)  (13) 
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The  resulting  centrifugal  tensile  stress  near  the  pin  locations  is  set 
equal  to  the  yield  stress  (with  an  appropriate  factor  of  safety). 

F 

0  =  Tirr  =  (h  +  .636a)  (2405a  +  1202t)  (14) 

y  2bt 

The  transmitted  torque  due  to  the  centrifugal  force  is 
F 

T  =  -£  a  =  (h  +  .636a)  (1202  abt2  +  2405  a2bt) 
c  2 

=  11,425  in. -lb  (15) 


By  the  method  of  Castigliano,  the  load-deflection  relationship  for  the 
elastic  link  (assuming  a  dinmetrinl  load  and  deflection)  is 


P  =  -  ~  ■■  g  5  =  (1.59)(106)  -5-  5 
3n(2a)  a 

where  P  is  the  value  of  the  bending  load  for  values  of  j  :  0. 
The  deflection  5  is  expressible  in  terms  of  6  as 


=  2(n-a) 


<6  -  9) 

2 1 r  sin  — — 


r  sin 


which  for  small  9  reduces  to 


e 

r  9  cos 


(16) 


(17) 


(18) 


The  torque  transmitted  by  elastic  bending  of  the  link  is  then 

T  =  Ph 
e 

2  3 

=  (1.59)(10G)li-4-  ^ 
a 


and  the  elastic  stiffness  is 


k  - 
e 


(1.59X106)  i— 


(19) 


(20) 


When  S  *  0,  the  torque  from  centrifugal  force  is 

T  =  2405  a2bt  (h  t  .636  a)  (21) 

c 
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where 


fi  =  r  cos 


a  -  r  sin 


(0  +0) 

o 

2 

(0+6) 

o 


(22) 


(23) 


The  centrifugal  torque  can  now  be  written  in  terms  of  the  twist  angle  8  as 

2  2  3/2  2  (V8)  'V9' 

T  =  2405  bt  (a  +n  )  sin  —  cos  — — — 

C  2  J 


+  0.63  sin 


(0  +0) 
o 


(24) 


The  centrifugal  stiffness  is 

*c  2  2  3  2  2  <V9»  (V8) 

k  =  =  2405  bt  (a  +h  )  1.9  sin  — —  cos  — St; - 

coo  22 

(0  +  0)  _  (0  *0)  ,  (6  +0) 
o  2  o  .  3  o 

+  sin — - -  cos  - - - jin  - - - 

2  2  2 

At  0  =  0,  the  centrifugal  stiffness  reduces  to 


(25) 


k  (6=0)  =  2405  bt  (2nh2  ♦  1.9  n2h-a3) 
c 


(26) 


The  overall  torsional  stiffness  of  the  coupling  is  the  sum  of  the  elastic 
and  centrifugal  stiffnesses,  which  is  set  equal  to  zero  for  6  =  0  to  give 


6  h^t3  2  9  3 

k  =  (1.59) (10  )  +  2405  bt  (2nh  +  1.9  a~h-a  )  =  0 

o  3 

a 


(27) 


A  coupling  can  now  be  designed  by  trial  and  error,  satisfying  equations  27, 
15,  and  14,  simultaneously. 

In  addition,  it  may  be  desired  to  introduce  n  size  constraint  equation. 

For  the  example  being  presented  here  with  5000  hp,  and  13,800  rpm,  a 
reasonable  constraint  on  the  diameter  of  the  coupling  might  be  taken  as 

a  +  h  =  I  in.  (28) 

A  solution  which  satisfies  equations  27,  15,  11.  and  28,  simultaneously, 
is  found  to  be 

a  =  3  in.  b  =  0.57  in. 

h  r  1  in.  t  =  0.302  in. 
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A  similar  procedure  would  be  required  In  order  to  design  a  coupling  for 
another  set  of  horsepower,  speed,  and  size  requirements. 


As  will  be  shown  lmawdlately  below,  an  advantage  of  the  pinned  link  cou¬ 
pling  over  the  elastic  link  coupling  described  above  is  that  a  very  simple 
design  equation  can  be  developed  for  the  pinned  link  coupling,  which  elim¬ 
inates  the  necessity  for  the  method  of  solution  just  described. 


Pinned  Link  Design  Analysis 

The  investigation  of  pinned  link  couplii^s  which  have  zero-  or  low- 
torsional  stiffness  begins  with  the  analysis  of  the  Russian  coupling  shown 
in  Figure  12.  Its  reduction  to  a  basic  kinematic  -  ->del  with  one  degree  of 
freedom,  6,  is  shown  in  Figure  14.  From  the  princ  pies  of  virtual  work 
and  conservation  of  energy  (assuming  no  friction),  the  following  equations 
result: 


8wf  =  F  6r 
6 Wfi  =  T  60 


T  = 


(29) 


Since  the  system  lr  Figure  14  has  one  degree  of  freedom,  then  R  =  f(0) 
where 


f(0)  =  r  cos  -  ♦ 


2 

r 


(30) 


Taking  the  variation  of  f(6)  with  respect  to  6  gives 


6R 


i  ,  .  e 

--  (r  sin  -♦ 


r  sin  9 


rt 2  2  ■  r 

•/  I  -  r  sin  ^ 


(31) 


The  force  F  caused  by  the  rotutlon  of  the  coupling  about  0  is 

2 

F  =  m  R  x 

Substituting  equations  30  through  32  into  equation  29  and  letting 


(32) 


(33) 


gives 


T 

.,2  2 

Nmr  u 


=  ![ 


sin  9  +  Z(sln  9/2)  + 


sin  9(cos  9/2) 
2Z 


■] 


(34) 
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where  N  is  the  number  of  linkages. 

Valves  must  be  determined  for  r,  t,  and  m  which  will  give  the  coupling 
zero  stiffness  at  operating  torque  and  speed.  If  the  ratio  i/r  is 
defined  as  a  new  parameter  equation  34  becomes 

— Tj— 2  =  (35) 

Nmr  x 

which  results  in  n  one-parameter  family  of  curves. 

The  function,  f(?,0),  plotted  against  8  for  various  values  of  r.  is  shown 
in  Figure  15.  As  can  be  seen  from  the  graph,  values  of  5  -  .  87  and 
9  =  100°  give  nearly  zero  stiffness  over  several  degrees  of  couplii>g  rota¬ 
tion. 

Since  the  operating  torque  and  speed  are  given  for  the  helicopter  problem, 
the  value  of  either  m  or  r  can  be  chosen.  The  other  value  can  be  calcu¬ 
lated  from  equation  35  where  5  -  .  87  and  8  ^  100°. 

Design  Equation  1 


- =  f  ( .  87 , 100)  ,  2.  115  >  1U~6  (36) 

Nmr  l 


where 


(37) 


m  lb  t  in. 

m 

Equations  36  and  37  are  now  genernl  design  equations  for  a  zero -torsional- 
stiffness  coupling  for  any  speed  and  torque.  It  is  necessary  to  see  if 
values  for  m.  r,  and  l  can  be  found  which  are  physically  realizable  for 
the  helicopter  coupling  design. 


For  high-speed  applications,  it  is  very  likely  that  the  mass  of  the  links 
which  has  been  ignored  up  to  now  will  have  a  significant  effect  on  the 
stiffness  characteristics  of  the  coupling.  Fortunately  the  mass,  mj,  of 
each  link,  and  the  mass,  mn,  at  point  I)  in  Figure  16  can  be  combined  for 
steady-state  analysis  (x  = constant) .  Figure  17a  shows  the  centrifugal 
force  F  resulting  from  the  distributed  mass  mj  of  a  link.  If  mass  mj  is 
equally  redistributed  to  points  A  and  B  (Figure  17b) ,  the  center  of 
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Figure  16.  Coupling  Mass. 
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gravity  of  the  link  does  not  change;  thus  the  force  F  does  not  change. 
The  force  vector  can,  however,  be  resolved  into  the  two  components  shown 
in  Figure  17b. 


fa  fai  '  FA2  =  I  "j  r  l2 


B  1  „  2 

T  r  2  mL  R  1 


Figure  17c  shows  the  two  links  combined,  giving 


F„  =  m .  R  i 
n  t 


'a  2  r  ‘ 


Including  m^  with  m^  gives 


fb  Kr>  •  ■nRii  ^  •  VR'2 


The  forces  F^i  and  F,\2  cannot  contribute  to  the  torque  T  since  they  are 
directed  radially  outward  from  point  0  and  have  no  moment  about  0. 

This  leaves  only  force  F«  to  produce  torque  T.  Thus,  the  relationship 
between  torque  and  twist  angle  can  be  determined  from  design  ««qun t  ion  1  by 
setting  n  a  m^  •  mn.  The  proof  of  this  statement,  given  below,  comes  from 
the  principle  of  virtual  work. 


where 


thus 


•  u. 

tHA 

6w 

T*  v 

*  *A 

Vr 

•  h 

n 

VR 

lr  ‘It 

1  fa  r  •  fb  ^ 


T  r T-r  ( m  .  •  m  )  Ri 
B  t  B 


2  -  It 
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Substituting  equations  31  through  35  into  the  above  equation  gives 


2  2 

X(mf+mD)r  i 

*  D 


=  f(3.f) 


which  is  the  same  as  design  equation  1,  with  m  =  m^.  ♦  m^. 

Therefore,  to  include  the  mass  of  the  links.  m£,  with  mass  mg,  simply  add 
to  mg  It  should  be  noted  that  adding  the  mass  of  one  link,  m  ,  at 
point  B  will  represent  the  total  effect  of  both  links."” 

Due  to  strength  considerations,  it  may  be  required  to  make  the  links  more 
massive  than  called  for  in  the  design  equations  36  and  37.  This  may  very 
likely  be  the  case  in  high-speed  turbine  applications.  Through  the  use 
of  a  form  of  counterbal ancing ,  the  links  can  be  made  as  massive  as  desired, 
without  altering  the  prescribed  stiffness  properties  of  the  coupling. 

After  mass  m  has  been  determined  from  design  equation  I,  it  is  necessary  to 
counterbalance  any  extra  mass,  me.  where 

Design  Equation  III 

m  =  m ,  -  m_  —  m  ( 38 ) 

e  l  B 

A  negative  value  for  me  means  the  link  is  too  light  and  +me  must  be  added 
at  point  B  or  to  the  eg  of  each  link.  If  me  is  positive  in  equntion  38. 
then  mass  mQ  must  be  ndded  at  point  E  to  each  link  at  a  distance  rp  as 
shown  in  Figure  18.  such  that 


Design  Equation  IV 


m  r  =  -  m  * 
i  e  2  e 


a  2r 


Coupling  Parameters  for  Helicopter  Analysis 

e. 

Representative  values  of  turbine  power  and  speed  from  Swick  and  Skarvan" 
are 

Normal  operating  power  =  5270  hp 
Normal  operating  speed  -  13.820  rpn 


Equation  *10  shows  that  a  zero-torsional-stiffness  coupling  must  be  designed 
to  transmit  n  torque  of  2004  fr-lb  at  13,820  rp:a  (1447  rad  sec). 


hp(33,000)  5270(33.000) 


(ft-lb)  2  ii 


(rpm) 


2^(13,820) 


=  2004 
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The  else  end  weight  of  the  coupling  can  diminish  the  helicopter  payload; 
therefore  the  value  of  r,  which  has  a  direct  bearing  on  the  coupling  flange 
•ize,  will  be  kept  to  a  minimum.  A  reasonable  choice  for  r  which  should 
allow  sufficient  clearance  of  various  existing  helicopter  hardware  is 
3  inches. 

For  design  equation  I  (36)  with 

r  =  3  in. 

X  -  13,820  rpm 
T  =  2004  ft-lb 
N  =  2(symmetric  coupling) 

8  =  . 87(zero  stiffness) 

n  value  of  m  =  .2385  lb«  is  obtained.  The  value  of  l  from  design  equation 
II  (37)  is  2.61  in.  All  kinematic  design  parameters  for  a  zero- torsional 
stiffness  have  been  determined  and  are  shown  in  Figure  19. 

Helicopter  Drive  Train  Parasieters 

The  helicopter  drive  train  shown  in  Figure  20  can  be  modeled  with  the  fol¬ 
lowing  parameters; 

Rotor  inertia  I 

R 

Engine  Inertia  1^ 

Speed  reducer  Inertia  1^ 

Mast  shaft  stiffness  K 

s 

Drive  shaft  stiffness 
Rotor  dnmplng  C 

La 

Engine  dnmplng  C 

e 

Mast  shaft  damping 
(see  Figure  20) 

For  systems  with  speed  reducers,  It  is  common  practice  to  reference  all 
springs  and  inertias  to  one  side  of  the  gearbox  (see  Appendix  II).  For 
this  helicopter  analysis  the  engine  input  side  will  be  used. 
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Numerical  values  for  the  above  Inertias  and  spring  constants  taken  from 
Swic)  and  Skarvan^  are 


xr  • 

3.888 

slug— ft2 

I  = 
e 

0.500 

slug-ft2 

XG  = 

0.144 

slug-ft2 

K  = 

212.1 

ft-lb/rad 

s*e 

t-  « 

ii 

21,417 

ft-lb/rad 

(all  values  referenced  to  engine  input  si>  e  of  gearbox) 

It  will  be  assumed  that  the  aerodynaaic  rotor  drag  and  engine  drag  are 
proportional  to  angular  velocity  squared.  Nearly  all  of  the  shaft  horse¬ 
power  of  the  engine  is  dissipated  in  the  aerodynamic  drag  of  the  rotor, 
where 

■  2 

T  =  C  0 
DRAG  L  R 

For  steady-state  operation  with  T  =  2004  ft-lb  and  6  =  6  =  1447  rad/sec, 

ft  ® 

c  =  9.564  X  lo”4  ft-lb-sec2 

La 

A  value  for  the  engine  damping  coefficient  Ce  is  determined  by  assuming 
a  25-30  percent  total  power  loss  due  to  turbine  drag;  thus,  total  hp  sj 
7400. 

From  equation  40  an  engine  drag  torque  of  about  800  ft-lb  is  obtained, 
giving 

C  8  =  800  ft-lb 

e  e 

or 

Ce  =  3.82  x  10-4  ft-lb-sec2 

The  speed  reducer  Inertia  is  small  compared  to  the  engine  inertia,  and  the 
drive  shaft  stiffness  is  much  larger  than  the  mast  shaft  stiffness  which 
is  effectively  reduced  by  the  speed  reduction  ratio. 

Since  the  investigation  of  torsional  vibrations  is  not  directly  concerned 
with  vibrations  of  very  high  frequency,  the  drive  shaft  can  be  assumed  to 
be  rigid,  allowing  the  speed  reducer  inertia  to  be  combined  with  the  engine 
inertia  giving  a  no*’  value  1^: 
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The  Model  for  the  helicopter  drive  train  with  Kt  -  ®  and  =  Ie+  lc 
reduces  to  the  two-degree-of -freedom  system  shown  in  Figure  21  (System  1). 

The  mast  shaft  damping  is  assumed  to  be  10  percent  of  critical  damping. 


where 


C  =  ( VK  I  >  =  2.24 

s  10  s  eff 


r  ID 
e  R 


I  +Ir, 
e  R 

The  drive  train  model  in  Figure  21.  with  the  zero  torsional  stiffness 

coupling  developed  in  the  previous  section  Inserted  along  the  drive  shaft, 

has  an  additional  degree  of  freedom.  The  inertia  of  the  driving  flange  of 

the  coupling  is  combined  with  the  engine  inertia,  since  Kt  is  assumed  to 

be  infinite,  giving  Ifc  =  !e*  !fe-  The  driven  flange  of  the  coupling  is 

combined  with  the  speed  reducer  inertia  to  give  i'  =  I_  +  I  . 

zr  u  ir 

Figure  22  shows  the  three-degree-of -freedom  system  incorporating  the  cou- 
pl lng  (System  2) . 

System  Simulation 


The  open-loop  helicopter  drive  train  model  without  a  zero-torsional-stiff¬ 
ness  coupling  can  be  dynamically  described  by  the  following  two  differen¬ 
tial  equations  of  motion: 


Ft  =  i'? 
e  e  e 


Fr„  =  10 

R  R  R 


or  by  inspection 


I  6  =  -K  (0  -6_)  -  C  0  -  C  (0  -0 „)  +  T(t) 

ee  seR  ee  seR 


Vr  *  WV  -  CL<t)9R  ‘  WV  « 

The  most  common  procedure  for  numerically  solving  these  equations  is  to 
first  use  two  dummy  variables  and  where 


V  ®R 


and  transform  equations  41  into  the  following  set  of  first-order  differ¬ 
ential  equations. 
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! 


K  -  7  «.<vv  -  cA  -  v  w  *  T<,» 

e 

S,  -  t;  ‘VW  -  Vt>xi!  *  ".‘VS.’1 


\  -  \  <«> 

The  equations  of  motion  for  the  helicopter  drive  t'ain  with  a  zero  tor¬ 
sional  stiffness  coupling  are  much  more  complex  th  n  the  above.  The  most 
effective  form  of  solution  is  through  the  use  of  analytical  dynamics,  more 
specifically  the  method  of  Vance  and  Sitchin*  for  directly  obtaining  first- 
order  difference  equations  from  Hamilton's  principle. 


The  method  will  first  be  used  with  the  two-degree-of-freedom  system.  The 
resulting  equations  will  be  used  for  the  actual  numerical  integration  due 
to  its  generally  improved  accuracy  over  other  methods.^ 

Hamilton's  principle  for  nonconservative  systems  is 


where 


2  K 


8l  =  6  /  L  dt  ♦  /  (  Z  F  6 

t.  t.  1=1  qi 


)dt 


0 


L  =  T-V  =  the  Lagrangl an 
T  =  the  kinetic  energy 
V  =  the  potential  energy 

=  the  nonconservntive  forces 


Hamilton's  principle  is  satisfied  by  the  set  of  k  second-order  Lagrange's 
equations 


d  &L  dL 

dt  (7T)  T“  =  Fi 


(43) 


’i 


'1 


The  generalized  coordinates  for  the  two-degree-of-freedom  system 
(Figure  21)  are 


K  =  2 
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For  the  three-degree-of-freedom  ayStea  (Figure  22),  the  generalized  coor¬ 
dinates  are 


V  e-  8r 


K  =  3 

Equations  of  Motion  for  System  1 


The  Logranginn  is 


L  =  i  i'  e2  +  I  I  e*  -  i  K  (9  -0  )2 
2  e  e  2RR  2seR 

The  elements  of  Lagrange's  equations  43  are 


cL 

e 

Jl.  ij 

de. 


=  i  e 

e  e 


R  R 


(44) 


it  -  V VV 

e 
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Momentum  is  defined  as 

_  ^ 

P‘  ' 

For  the  two  generalized  coordinates  8^  and  0^/  the  momenta  and  their 
respective  time  derivaties  are 

p  =  _ z.  -  i  e 


(45) 


(45) 


dS 


e  e 


P  =  =  I  0 

R  dOR  R  H 

p  =  ±  A 
*  dt  ae 


p  =  ±  <iL) 
’  dt 


(46) 


(47) 
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The  generalized  forces  can  be  determined  from  the  principle  of  virtual 
work: 


for  e 

e 


3q 


1 


iW  =  F  5  6 
e  e  e 


•  2 

f  =  -c  e 

e  e  e 


C  (0  -8  )  4  T  ( 
s  e  R  e 


t) 


and  for  0 

H 


R 


.  o 

Vr  * 


C 

s 


<0  -6 
e 


R 


) 


(•18) 


Using  equations  45,  47,  and  48,  Lagrange's  equations  43  can  be  written 
In  the  canonical  form: 


a 

P 

e 


K  (6  -0  )  -  C  0  -  C  (0  -e„> 
s  e  R  e  e  s  e  H 


T  (t) 
e 


(49) 


a  a 


P„  =  -K  (0  -0n)  -  c,  9,  4  C  (0  -9  ) 

R  scR  I,  R  s  e  R 

For  dynamical  systems  the  P^'s  are  always  linear  in  the  q^'s;  therefore, 
from  linear  algebra 


rPt}  =  [  A  | (  q  i } 


(50) 


or 


-1, 


fqi }  n  [A  ]  »’P  j  (assuming  A  is  nonsingulnr)  (51) 

From  equations  44  and  50  the  remaining  canonical  equations  are  formed. 
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0 
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1*R 

(52) 


Equations  49  and  52  form  a  set  of  four  first-order  differential  equations 
in  coordinates: 
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e  e 

P»  6n 

R  R 

These  equations  can  be  numerically  solved  by 

1.  Specifying  initial  conditions  for  P  ,  P  ,  0  ,  and  0 

°0  K0  °0  R() 

•  • 

2.  Finding  values  for  8  and  8„  from  equation  52 

e  R 

•  • 

3.  Finding  values  for  P^  and  PR  from  equntion  -19 

•  •  •  • 

4.  Numerically  integrating  P  ,  P  ,  9  and  8  using  routines  such  ns 

Euler  or  Runge-Kuttn  c  c 

Using  Euler  integration, 


8  =  8  4  8  .It 

B 
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en+l 

e 

n 

e 

] 

• 

8 

8 

4  e 

H  4 
n+1 

R 
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R 
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• 

P 

P 

4  P 

°n+l 

e 
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e 

1 

• 

P„ 

IV 

+ 

n+1 

R 

n 

R 

I 

‘n+l  = 

t  4 

n 

it 

(53) 


5.  Using  the  n4l  values,  return  to  step  2  nnd  repeat  procedure  for 
desired  length  of  time. 

Equntion  of  Motion  for  System  2 

In  order  to  find  an  expression  for  the  Lagrunglnn  of  the  three- degree-of- 
frecdom  system,  the  linear  nnd  angular  velocities  of  each  mass  nnd  inertia 
must  be  determined,  t  For  the  engine  and  rotor  inertias  nnd  i  ,  the  angu¬ 
lar  velocities  ore  8e  nnd  respectively.  The  velocities  of  the  coupling 
links  are  functions  of  $e  nnd  8  ns  well  ns  @e  and  6.  Referring  to  Figures 
22  nnd  23,  the  kinetic  and  potential  energies  nrc 


1  •  2  1  •  2  1  /  •  •  2  1  2  1  2  2 
T  -  2  V2  Vr*2  IfR(VC)  *  2  NVb  *  2  "V'Y'V 


v  =  4  k  <e  -e-eD)2 

2  s  e  R 


1  2  2 

♦  -  NI  .(u; ,  .  ) 

2  L  ^  l2 


(51) 
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whore 


v .  s  linear  velocity  of  eg  of  link  1 
ll 

w.  =  angular  velocity  of  link  1 

11 

v  .  =  linear  velocity  of  eg  of  link  2 

12 

id  |  =  angular  velocity  of  link  2 

l2 

(aee  Figure  23) 

For  the  static  analysis  of  the  coupling,  it  was  sh>wn  that  and  mg  can  be 
combined,  but  when  6  is  not  constant  the  mass  distribution  of  the  links 
will  have  some  effect  on  the  system  response.  A  significant  effect  would 
only  be  expected  when  the  links  are  massive  (low  speed/high  power)  but  will 
be  left  in  the  analysis  for  generality. 

The  expressions  for  the  v's  and  it's  in  the  kinetic  energy  equation  53  ore 
developed  in  Appendix  III  and  are  given  below. 

2  *2  2  •  •  2 
=  R  +  R  (6  -e/2) 

B  6 

vi  =  vi  +  VA  c  I  fr2(fl^-iefl  ♦  82/2)  -  r  R  0(sin  8/2) 

♦  R2  +  2rR(8  -8/2)2  cos  8/2  +  R2(8  -8/2)2}  (55) 

©  0 

2  2  2  9.2  2  •  •  *2  •  • 

<d.  =  uj.  +  <d  .  «  -=■  fr  (8  -8  8*8/2)  +  r  R  8(sin  8/2) 

1  i  l2  r  e  * 

+  R2  -  2rR(8  -8/2)2  cos  8/2  +  R2(8  -§/2)2} 

e  e  1 

s 

where  R  is  the  time  derivative  of  R  in  equation  30 

r(sin  8/2)  +  —  r  s,in  6. .  ..  J  (56) 

2  J /2-r2(sin  8/2) 

The  Lagrangiun  can  now  be  written  in  terms  of  the  generalized  coordinates 
as 


46 
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*•  =  !  V»* 

1  *o  2  •  •  2 

+  jn  +  r  <ee-e  2 ri 

-is  M  [ r  R  0(sin  0/2)  -2rR(0  -  0  2)2  cos  0/2] 
z  z  c 

+  7  N  M_[r  (9*  -  e  e  ^  0^  2)  ] 

Z  %S  ©  <? 


-  k  <e  -  e  -  9  ) ' 

SC  R 


where 

2 

Ml  =  raQ  +  2  -r  21  £  l 

‘<2  “  mX  2  -  21 X  S 


»3  =  "X  2  -  2'xl 


The  potential  forces,  dL'*q  .  in  Lagrange's  equations  are 


>L 


dR 

1  ^R 


^R 


2R^(6e"S2)  1 


-  \  NSIjIre  (sin  0-2)  4  1  r  R  ?(cos  0  2) 

-  2r  ^  (0  -  4  2) 2  cos  6  2 


r  R(6  -0  2)  sin  fc  2]  ♦  K„(0  -  0  -  6„)  for  0 


2  e 


|b-  =  -  K  (9  -  0  -  »  )  for  0 
dc  s  e  R  e 


■£-  =  K  (0  -  0  -  SD) 
30,.  s  e  R 

K 


for  0 


R 


The  momenta,  P  ,  in  Lagrange's  equations  are 

dl  •  • 

PT  =  —  =  -  1,(6  -0) 

T  aI  fR  ° 


1  •  2  •  • 

*  -  VM1|2R  -  H  (P^  -  0/2)) 


aC. 


(57) 


(58) 


48 


-  i  NMgfre  ^  (sin  0/2)  +  r  R(sin  0/2) 
d0 

+  2  r  R(0  -  0  '2)  cos  0/2]  +  i  NM„r2(0-  6  )  lor  0 
e  2  3  e 


P  =  =  I  +  1.  (0  -  0) 

•  se  e  fR  e 
e 


+  |  NM1[2R2(0e  -  0/2) ] 


+  s  NM_[4rR(6  -  G  2)  cos  0  2] 
+  \  NM3[2r2(Ge-0/2)] 


P  =  =  i  o 

R  d-  RR 

R 


for  9_ 


The  expressions  for  the  partial  derivatives  of  R  and  R  used  In  equations  58 
In  terns  of  the  generalized  coordinate  0  are 


AR  1  f”  /  .  a  r  sin  6  j 

*s  ’  'Ur  ‘*‘n  9  21  ♦ 

2 v  l  -  r  (sin  0  2) 


dn  i  .  r 

35  =  "4  8  [_r<co»  e/2)  ♦ 


r  cos  8 _ 

V  - r2(sin20  2) 


4  2 

r  sin  0 

.2  2  2  3  / 2_ 

4( l  -  r  (sin  9  2) ] 


dR  *R 


The  tine  derivatives  of  the  monentn,  P^'s,  need  not  be  taken, 


d  /& L  ' 


1  =  It 


for  9 


•  d  *1,  ' 

P  =  -£■  l !  for  e 
e  dt  t 


•  d  (dl  ,  . 

PR  •  3i  fe-J  ,or  \ 
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The  generalized  forces  F,  determined  from  the  principle  of  virtual  work 
are 


Av  =  F  A0 

e  e  e 

f  =  t  (t) -c  <0  -e-e  >  -c  e  |e  |  for  e 

e  e  s  e  R  e  e  e  e 

where  T  (t)  is  any  tlaie  function  of  the  input  torque. 

%  ■  VS 

PR  *  vs-^-v-v-'W  ,o*  s 

where  C  (t)  is  any  time  function  of  the  rotor  lift  coefficient. 
L 

Aw_  =  F  60 
T  T 


f  =  c  ( 9  -e-e>cJ-KJ0) 

T  s  e  R  D  fX 


for  0 


<61) 


where  Cp  is  viscous  hinge  friction  in  the  coupling  snd  K^O)  represents 

any  additional  springs  that  might  be  added  to  the  coupling.  and 
are  the  shaft  and  engine  damping  coefficients  respectively. 


Using  equations  38,  60,  and  61 ,  'Lagrange' s  equation^  can  be  written  in 
canonical  form: 

•  Al 

P  =  -rr-  +  F 

e  00  e 
e 


•  dL 
PT  =  35  +  FT 


(62) 


'R  =  W  '  'R 


From  equations  50  and  58, 


=  [A] 


(63) 


P„  =  I  »n 
R  R  R 


where  the  elements  of  the  A  matrix  are 


50 


■n  =  5NrV2*|f+5R2> 


•  a 

♦  Mgf-r  ^  (sin  0/2)  -rR(sin  0/2)  >  rR(cos  0/2)+  Mgr2)  +  1^ 

S0 


=  Nfi^R  +  2M2rR(cos  0/2)  +  Mgf  )  -  I, 


a21  =  °12 


°22  =  *e  “  *fR  "  2012 


R  =  R  6  *  'SJS 

Upon  inversion  of  equations  63,  the  following  expressions  for  the  three 
generalized  velocities  result: 

®=  (n22PR  "  al2Pe)  (alln22  “  n21a12) 

K  =  (ailPR  "  a21Pe)  (0lin22  "n21°l2) 

*R  *  PR  *R 

Equations  62  and  64  form  a  set  of  six  first-order  differential  equations 
in  coordinates: 


These  equations  can  be  numerically  solved  by 

1.  Specifying  Initial  conditions  for  P  ,  P  ,  P  ,  0  ,0  and  0 

T0  *0  R0  *0  °  R0 
•  •  • 

(initial  conditions  for  0  ,  0  ,  and  0  can  be  used  with  equn- 

®0  °  Ro 

tions  63  to  produce  P  ,  P  ,  and  P  ) 

e0  T0  R0 


2.  Finding  values  for  0  ,  0  and  0  from  equations  6-1 

e  R 
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3.  Finding  values  for  P  ,  P  and  P  from  equations  62 

®  T  R 

4.  Numerically  integrating  F>e,  PT>  PR,  8^  6,  8R 
Using  Euler  integration, 


e 

en+l 

=  e 

en 

+  e  it 

en 
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®R 

ml 

-  eR 
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=  P 
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t  P  fit 
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=  PTn 

♦  PT  At 
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PR  t 
ml 

"  PR 

n 

+  PR*1 
n 

t  i  = 
ml 

t  + 
n 

At 

5.  Using  the  ml  values  from  step  4,  return  to  step  2  and  repeat 
procedure  for  desired  length  of  time. 

THE  SYNTHESIZED  COUPLING  DESIGNS 
Description  of  Pinned  Link  Coupl lng 

Three  views  of  the  synthesized  zero  torsional  stiffness  (/TS)  pinned  link 
coupling  for  helicopter  application  are  shown  in  Figures  24,  25,  and  26. 


The  various  components  of  the  coupling 

are 

A 

driving  flange 

F 

connecting  tubes 

B 

driven  flange 

C 

deformable  stops 

C 

1  ink 

11 

link  pin  and  fastener 

D 

connecting  pin 

I 

torsion  bar  (not  shown) 

E 

alignment  shaft 
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Figure  25.  ZTS  Coupling  (Side). 
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a> 

I 

I 

I 


Figure  26.  ZTS  Coupling  (Dock). 
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The  connecting  tubes  nre  part  of  the  flanges  and  are  designed  to  fit  inside 
the  drive  shaft  of  the  helicopter.  A  froe  running  alignment  shaft  may  bo 

used  to  prevent  lateral  motion  of  the  flanges.  A  soft  torsion  bnr  can  be 

Inserted  through  the  center  of  the  alignment  shaft  which  will  close  the 
coupling  when  the  system  is  at  rest.  This  will  prevent  bnnging  of  the 
coupling  against  the  stops  during  initial  start-up.  The  ends  of  the  tor¬ 
sion  bnr  can  be  attached  at  any  point  inside  the  drive  shnft. 

Figure  27  shows  the  type  of  link  used  in  the  /.IS  coupling.  The  dimensions 

of  the  link  nre 


L  =  2.61  in. 

r  =  .375  in. 

P 

r  -  .625  in. 
r2  =  .6875  in. 

The  areas  shown  in  Figure  27  are 

Ax  =  2(11  -  h)  t  =  .375  in.  “ 

A  -  projected  bearing  surface  2tr 
r.  .6875  in.2 

A„  =  2(r  -r,  ) t  s  .25  in.“  A.  *t2  .1*118  in.“ 

3  2  1  -1  p 

The  bending  moment  caused  by  the  links  being  in  different  planes  can  be 
eliminated  by  replacing  each  link  connected  to  flnnge  n  with  two  similar 
links  of  thickness  t  2.  The  new  links  go  on  either  side  of  the  link  con¬ 
nected  to  flange  A. 

Determining  Counterweight  nu. 


t  =  .5  in. 

H  1.375  in. 

h  .625  in. 

r  1.3  in. 
e 


The  link,  before  counterweight  ing ,  is  divided  into  three  sections.  .Sec¬ 
tion  1  is  all  mass  symmetrically  distributed  about  point  H,  which  for  each 
link  is  mn  2.  Assuming  that  Jjhe  pin,  bearing,  aim  link  all  have  the  den¬ 
sity  of  steel  (,  =  .28  lb  in.  )  gives 


r  t  -  0 . 208  1 b 
2  m 


Section  2  is  the  portion  of  the  link  which  has  the  cross  section  shown  in 
Figure  27c.  The  center  of  mass  for  this  section  is  located  at  L  2;  thus 
it  can  be  treated  simply  as  m^. 

m„  m,  -  o (H-h) (f-2r„) (t)  0.130  lb 

2  t  2  m 
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Section  3  le  all  mass  symmetrically  distributed  about  point  A. 

m,  *pm-:?<2t)  =  .416  lb, 

J  z  in 

This  mass  does  not  affect  the  torque  characteristics  of  the  coupling  and 
can  be  neglected  when  determining  ma>  From  design  equation  (38), 

in  s  Bj  +  in  “  in 
e  l  B 

The  value  of  m  has  previously  been  determined  on  page  35  to  be  .2385  lb  . 
The  value  for  m  from  equation  38  Is  therefore 

m  =  .307  lb 
e  m 

Design  equation  IV  (39)  with  r^  =  i/2  gives 


Stress  Analysis 

The  coupling  geometry  used  for  the  stress  analysis  is  shown  In  Figure  2Hn 
where 


R  = 
h  = 


Of  = 

Figure  28b  shows 


r(cos  0/2)  ♦  J &  -  r2(sin20/2) 
r(sln  0/2) 
sin  1  h/i 

the  forces  acting  on  each  link  and  pin. 


Mass  m^  was  determined  so  that  me  would  not  contribute  to  the  torque  char¬ 
acteristics  of  the  coupling.  This  was  done  by  making  (me  2)  *  my  for  each 
link  have  a  resulting  center  of  gravity  at  point  A.  Figure  17  shows  that 
the  link  mass,  mj,  can  be  equally  distributed  to  points  A  and  n.  The 
force  acting  at  point  A  caused  by  mf,  m  ,  and  m  is 

_  ,  "•  "X  2 

Fi  =  (m&  *  T  +  T]  ™ 

The  force  at  point  A  caused  by  all  mass  located  in  section  3  is 


where 


m3  tw  Som^t 
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The  noncounterbnlnnced  mass  m  produces  a  force  at  point  B  with  magnitude 

2 

F,  =  milt 

u 

which  causes  each  link  to  have  n  tension  force 
1  2 

P.  --  -  m  »;<x  cos  a) 
a  2 

The  total  force  F^  acting  on  the  bearing  at  point  A  is  the  vector  sum  of 
,  Fq ,  and  F4  : 

Fa  =  f|(F  tF  )(cos  9  21+F^  cos  o]2 
♦  r(F1+F,,)  (sin  0  2)  -F  sin  a]2)1  2 

The  shenr  force  on  the  pin  at  B  caused  by  the  force  acting  between  the 
links  is 


(due  to  symmetry,  only  normal  force  between  the  links  exists). 


Equntion  39  shows  that  the  moment  caused  by  on  one  side  of  Point  A  is 
equnl  to  the  moment  caused  by  mc  2  on  the  other  side  of  A.  The  value  for 
the  moment  where  the  open  section  of  the  link  begins  (point  N)  is 


M 


h  I 


lit 


2 


I 


projected  mo-  force 
ment  arm 


The  mass 


m,  causes  a  tension  force  F_  in  each  link  where 
l  5 


1  2 
F_  -  7T  (m  *  ml)H(i  cos  a) 
O  £  II  * 


The  normal  stress  in  the  link  at  point  N  is 


c 


N 


F 

A 


5 

1 


The  stress  on  the 


bearing  at  pin  A  is 


60 


The  normal  stress  In  the  link  at  the  section  shown  In  Figure  27d  is 


and  the  sheur  stross  acting  on  the  pin  at  B  is 


The  value  of  cnch  stress  and  forco  is  a  function  of  6.  The  maximum  values 
are  computed  by  the  stress  program  in  Appendix  IV  and  are  given  below. 


Maximum  Stress 


c  ?  60,000  psl 

a.  —  25,000  psi 
A 

on  *  23,000  psl 
r  13,000  psi 


F.  17,200  lb 

Fb  5,700  lb 

F.  w  5,800  lb 
4 

F_  *  13,200  lb 
o 


M  «  4,800  in. -lb 

The  stress  distribution  about  pin  A  is  very  complex  and  depends  on  the 
rigidity  of  the  bearing,  flange,  and  link  as  well  as  the  moment  caused  by 
and  the  normal  stress  caused  by  m.  The  pin  has  been  made  large  to  help 
account  for  these  effects  and  to  increase  rigidity.  An  experimental  anal¬ 
ysis  is  required  in  order  to  determine  the  exact  mode  of  failure.  This 
region  is  not  critical  in  the  design  other  than  for  stress,  since  the  mass 
in  the  region  contributes  nothing  to  the  torque  characteristics  of  the 
coupl ing. 

The  steel  links  and  pins  should  be  ns  strong  or  rigid  ns  possible  to  avoid 
link  deflection  and  prevent  dnmnge  due  to  extreme  overloading  if  system 
failure  occurs.  AISI  2330  steel  drawn  at  400°F  hss  a  yield  strength  of 
195  ksl  and  a  Brlnell  hnrdness  of  425.  The  bearing  surfaces  are  subjected 
to  high  stress  and  are  difficult  to  lubricate.  Self-lubricnting  benrliurs 
made  with  woven  Teflon  have  been  designed  to  withstand  loads  to  80  ksi. 
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Design  of  Deformnble  Stops 


Stops  have  been  Included  in  the  coupling  design  tn  limit  6  between  1  and 
2  radians.  The  stops  can  come  in  contact  during 

1.  Start-up 

2.  Shutdown 

3.  High  speed/ low  torque  conditions 

4.  Low  speed/high  torque  conditions 

5.  Low  frequency /high  magnitude  excitation 

An  elastomeric  material  is  used  for  the  stops  to  reduce  Impact.  The  total 
contact  aron  of  each  set  of  stops  is  2  in.2.  The  spring  constant  of  the 
stops  is  designed  so  thnt  the  operating  torque  (2004  ft-lb)  at  zero  speed 
can  be  transmitted  with  a  .03  radian  coupling  rotation  from  initial  stop 
contnct. 

The  center  of  the  contact  region  is  3.25  in.  from  the  center  of  the  cou¬ 
pling;  thus  the  force  on  the  stop  is 

F  =  (2004  ft-lb) (12  in.  ft)  3.25  in.  -  7,100  lb 

The  linear  spring  constant  is 

K|  =  (7,400  lb)  (.03  rad  •  3.25  in.  rad) 

=  75,900  lb  in. 

The  torsional  spring  constant  1b 

Kt  =  (2004  ft-lb)/(.03  rad)  =  66,800  ft-lb  rad 

An  expression  for  the  approximate  value  for  the  deflection  of  the  elasto¬ 
meric  material  is 


where 


6 


PL 

AK 


AK 

I. 


Letting  I,  (the  initial  length  of  the  material)  be  equal  to  ,75  in.  gives 


(75,900  lb  in.  H.75  in.  )  ....  ...... 

h  - - -  1  ■  ■  ■  —  "  . .  - .  28,500  psi 

2  in. 


Low  density  polyethylene  has  an  elastic  modulus  range  of  11,000-38,000 
ns  i . 


Description  of  llMtlc  Link  Oouplli 

Two  views  of  the  synthesized  zero  torsionsl  stiffness  elastic  link  cou¬ 
pling  for  helicopter  application  ara  shown  in  Figures  29  and  30.  The  only 
two  aajor  differences  from  the  pinned  link  coupling  are: 

1.  Substitution  of  sesiicircular  elastic  links  for  the  straight 
pinned  links. 

2.  Reduction  of  the  angle  of  twist  (fron  stop  to  stop)  to  10  degrees 

As  shown  in  the  section  on  design  analysis,  the  determination  of  elastic 
link  dimensions  ia  an  Indirect  and  time-consuming  procees.  Link  dimensions 
are  given  on  page  25  for  the  example  taken  fron  Swlck  and  Skarvan. 

The  different  angles  of  twist  between  stops  for  the  pinned  link  design  and 
for  the  elastic  link  design  represent  two  different  philosophies  of  design, 
rather  than  special  requirements  of  each  type  of  coupling,  although  the 
small  twist  angle  in  the  elastic  link  coupling  eliminates  the  necessity  to 
avoid  negative  stiffness  at  large  angles  of  twist.  The  small  twist  angle 
also  eliminates  the  possibility  of  a  bending  failure  load  in  tha  elastic 
link. 


With  its  small  (10°)  angle  of  twist  between  stops,  the  elastic  link  cou¬ 
pling  acts  as  if  it  were  solid  (no  compliance)  during  acceleration  or 
deceleration  and  exhibits  zero  stiffness  only  at  the  steady-state  design 
conditions.  Thus  the  elastic  link  coupling,  as  designed  here,  is  best 
suited  to  providing  zero  stiffness  to  Improve  closed-loop  stability, 
without  the  capability  to  absorb  large  excltationa  or  dissipate  vibration 
energy  through  damping. 

The  two  designs  presented  here  (pinned  link  and  elastic  link)  can  be 
regarded  as  representing  the  two  extremes  of  design  philosophy  which  con¬ 
tain  the  zero  stiffness  concr^t.  The  questions  of  whether  pinned  links 
or  elastic  links  are  superior,  or  how  much  angle  of  twist  between  stops  Is 
optimum,  must  finally  be  answered  by  testing  for  a  specific  application. 


Provision  for  Da 


The  large  angle  of  twist  between  stops  Incorporated  in  the  pinned  link 
design,  along  with  the  available  space  represented  by  the  alignment  shaft, 
provides  an  excellent  opportunity  to  Introduce  significant  amounts  of 
shaft  damping  into  the  helicopter  system. 


Figure  31  shows  a  design  concept  developed  at  the  University  of  Florida 
to  accosq>llsh  this.  The  alignment  shaft  consists  of  a  hollow  steel  sleeve 
containing  an  annulus  of  elastomeric  material.  A  steel  rod  is  a  sliding 
fit  into  the  elastomeric  siaterlal.  The  rod  is  fixed  to  one  flange  of  the 
coupling  and  the  steel  sleeve  Is  attached  to  the  other  flange.  The  only 
bonds  between  steel  and  elastomer  are  at  areas  "A"  and  "B."  The  system 
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Flgur«  29.  Elastic  Link  Coupling,  Front  View. 
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Figure  31.  View  of  Proposed  ZTS  Coupling  Dampener  Configuration. 
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is  pretwisted  on  assembly  to  provide  an  ol  attic  clot  inf  torqua  which  kaapa 
the  coupling  agalnat  lta  atopa  until  centrifugal  fore*  froai  rotation  pulla 
it  opan.  Thua ,  "banging"  or  Impact  of  atopa  during  atart-up  ia  eliminated. 
Figure  32  shows  a  aaction  viaw. 

It  can  ba  aaan  that  toralonal  oaclllatlon  of  tha  coupling  producaa  tha 
following  af facta,  all  of  which  dlaaipata  anargy: 

1.  Twisting  (sliding)  of  tha  ataal  rod  with  raapact  to  tha  elasto¬ 
meric  annulus. 

2.  Twisting  (dafomation)  of  tha  elastomeric  material. 

3.  Twitting  (sliding)  of  tha  elastomeric  annulus  with  raapact  to 
tha  hollow  ataal  slaava. 

It  is  ah.yn  in  Part  1 1  that  shaft  damping  la  an  extremely  sansitlva  and 
influantlal  paraaatar  in  improving  cloaed-loop  toralonal  stability. 


RESULTS  OF  SYSTEM  SIMULATION  WITH  A  SYNTHESIZED  DESIGN 

6HH-LMP  VfbfcATi6M  I^lAYion - 


Tha  rasulta  of  tha  computar  simulations  daacribad  pravlously  ara 
praaantad  in  this  aaction.  Appendix  IV  contains  tha  computer  programs 
and  numerical  data  uaad  for  tha  simulation  of  tha  drive  train.  All  of 
tha  simulations  reported  in  Part  I  ware  dona  with  the  pinned  link  coupling 
without  damping  in  the  system  taken  from  Swick  and  Skarvan. 

Modes  of  Excitation 


Systea  1  (Figure  21)  and  System  2  (Figure  22)  ara  subjected  to  tha  fol low¬ 
ing  external  excitations: 

Teat  1— Variation  of  C. 

r  1  — —  L 

Ibis  test  models  tha  cyclic  change  of  tha  aerodynamic  rotor  drag 
caused  by  tha  rotor  blades  entering  and  leaving  tha  direction  of 
flight  as  wall  as  cyclic  variations  caused  by  pilot  excitation  ot 
collective  pitch.  Tha  lift  coefficient  of  tha  rotor,  C^,  la 
sinusoidally  varied  ±10  percent  about  its  average  value. 

C.  (t)  =  C,  (1  +  0.1  sin  uit) 

L  L 

uu  =  2.5  “*  10  Hz 
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Suction  View  of  Proponed  ZTS  Coupling 
D.impener  (  onf  igumt  ton  . 


f.H 


i  t  gore  :\2 


Tent  2— Variation  oi  T 

-  — —  —  ■  ■ .  i  ■—  — — «? 

Pilot  adjustment  o "  eng l no  speed  as  woll  as  possible  feedback 
Instability  is  represented  by  sinusoidally  varying  the  engine 
torque  110  percent  about  normal  operating  input  torque,  T 
(2804  ft-lb) .  * 

T  ( t)  -  T  (1+0.1  sin  at) 

e  e 

CL  =  2.  3  -  10H* 


A  helicopter  control  system  compares  the  engine  speed  with  s  desired  speed 
and  sltors  the  vnrioue  engine  torque  producing  parameters  accordingly. 

If  external  excitations  cause  large  variations  of  the  engine  speed,  con¬ 
trolling  system  stability  may  prove  to  bo  difficult.  A  beneficial  effect 
of  the  ZTS  coupling  la  to  effectively  isolate  the  engine  inertia  from  the 
rest  of  the  system,  resulting  in  a  reduction  of  the  engine  speed  variation 
due  to  external  excitation.  The  first  measure  of  system  response  is  A0e, 
the  engine  speed  variation  about  the  normal  operating  speed  (1447  rad/sec). 
The  second  meartre  of  system  response  la  ATS,  the  variation  of  the  torque 
transmitted  by  the  mast  shaft  about  the  normal  operating  torque  (2004 
ft-lb).  Hie  quantity  ATB  indicates  the  torque  overload  experienced  by  the 
mast  shaft,  gearbox,  engine  turbine,  and  rotor. 


The  speed  variation,  A0e,  has  units  of  rad/sec  and  the  torque  overload, 
AT^,  is  in  percentage  of  operating  torque. 


The  results  of  Test  1  given  in  Figure  33n  show  the  engine  speed  variation, 
A0O,  and  must  torque  variation,  ATg ,  ns  u  function  of  the  exciting  fre¬ 
quency  u).  Above  5  II/. ,  A0e  and  AT„  are  virtually  unaffected  by  the  exciting 
force.  At  slightly  above  3  Hz,  System  i  experiences  a  critical  frequency. 
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Fl*ure  33«.  Variation  of  C  (K  -  212.1) 


System  2  has  •  critical  region  below  2.5  Hz  where  the  coupling  hlta  against 
the  stops  during  each  cycle  of  excitation. 

e 

Due  to  the  nonlinearity  of  the  coupling  in  Systea  2,  the  values  of  A6#  and 

ATS  depend  not  only  on  the  exciting  frequency  but  on  the  state  (position 

and  velocity)  the  syaten  is  in  when  excitation  occurs.  If  either  Systym  1 

or  2,1s  excited  at  very  low  frequencies  u>  <  .  5  Hz,  then  the  values  of  8e 

and  will  be  nearly  the  sane.  For  this  type  of  excitation  the  coupling 

in  System  2  nay  close  against  the  atops,  causing  the  system  to  assume  the 

stiffness  of  the  mast  shaft,  K  . 

s 

Figure  33bshows  the  results  from  Test  2.  The  engine  inertia  is  directly 
excited  in  this  test,  causing  severe  oscillations  of  engine  speed  and 
mast  shaft  torque  when  excitation  occurs  at  or  near  the  critical  frequency. 

e  e 

A  time  response  .curve  9  f  0e  and  9p  during  start-up  and  shutdown  is  shown 
in  Figure  34.  9e  and  Cn  of  System  1  and  System  2  closely  follow  curve  1 
for  the  first  20  sefonds  (start-up  and  steady  state).  During  initial 
start-up,  9#  leads.9^  while  the  mast  shaft  "winds  up"  and  then  small  orcll- 
latlons  of  le  and  0R  about  curve  1  slowly  dampen  out.  After  20  seconds, 
when  T#  is  set  to  zero,  $e  and  eR  of  System  1  closely  follow  curve  2. 

Small  oscillations  again  occyr  as  the  mast  shaft  unwinds.  9r  in  System  2 
closely  follows  curve  2  but  9.  experiences  severe  oscillations  for  the 
first  1/2  aecond  wnen  the  system  is  shut  down.  This  is  due  to  the  coupling 
mass  traveling  radially  outward,  causing  the  coupling  to  quickly  close 
against  the  reverse  stops. 

Other  Systems 

Swlck  and  ;'ikarvan5  point  out  that  the  difference  in  the  natural  frequen¬ 
cies  of  helicopter  drive  traina  depends  primarily  on  mast  shaft  stiffness. 

A  common  range  for  the  natural  frequencies  is  2-6  Hz.  The  basic  helicopter 
design  presented  thus  far  has  a  critical  frequency  slightly  above  3  Hz. 
Increasing  the  mast  shaft  stiffness  of  the  helicopter  design  by  a  factor 
of  four  gives  the  system  a  natural  frequency  of  approximately  6  Hz. 

Figure  35  shows  the  response  of  Systems  1  and  2  with  increased  mast  shaft 
stiffness  to  the  excitation  of  Test  1.  At  the  critical  frequency  of 
System  1  (6  Hz)  the  torque  vnrintlon  is  over  22  percent  while  the  speed 
variation  la  over  11  rad/sec.  With  the  use  of  a  ZTS  coupling  (System  2), 
the  critical  frequency  of  the  system  is  shifted  below  4  Hz. 

•  t 

The  time  response  of  9e,  9p  and  mast  shaft  torque  to  a  6  Hz  excitation  of 
the  systems  in  Figure  35  is  shown  in  Figure  36.  The  zero  reference  speed 
in  Figure  36  is  1447  rad/sec  and  the  zero  reference  torque  is  2004  ft-lb. 

The  rotor  speed  varies  ±2.2  rad  sec  for  System  1  nnd  ±1.3  rad/sec  for 
System  2.  The  engine  varies  £12  rad  sec  for  System  1  while  the  engine 
variation  for  SyBtem  2  is  less  than  ±.5  rad/sec.  Variation  of  the  mast 
shaft  torque  is  ±23  percent  for  System  1  nnd  less  than  1  percent  for 
System  2. 


71 


Excitation  Frequency  u!Hz) 


2 


i 

V 


I 

£ 


Figure  34.  Tine  Response  of  Helicopter  Drive  Trein. 
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Figure  36.  Systea  Tise  Response  to  Excitation 
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Figure  37.  Variation  of  C  (K  -  1696.8) 


Systems  which  have  higher  shaft  stiffness  m iy  represent  other  turbine 
applications  such  as  pump  or  generator  systems.  The  final  investigation 
of  system  response  to  external  excitation  Is  with  mast  shaft  stiffness 
eight  times  that  of  the  original  helicopter  design.  Figure  1)7  shows  the 
system  response  from  Test  1.  The  critical  frequency  of  System  1  is 
slightly  under  7.5  Hz.  Using  a  ZTS  coupling  shifts  the  critical  frequency 
of  the  system  below  4  Hz. 


SUMMARY  OF  RESULTS 

1.  Torsional  vibrations  can  be  reducod  or  eliminated  with  a  ZTS  cou¬ 
pling.  The  factors  limiting  coupling  effectiveness  are 

a.  Very  low  frequency  excitation 

b.  Very  largo  torque  variations 

c.  New  critical  frequencies  in  range  of  excitation 

2.  Through  link  counterbalancing,  the  ZTS  pinned  link  coupling  can 
enter  the  domain  of  high-speed  turbo-powered  machinery. 

3.  For  systems  where  the  diag  torque  is  proportional  to  tho  angular 
velocity  squared  (helicopter  system),  the  ZTS  coupling  bocomes  speed  Inde¬ 
pendent  (see  design  equation  I  with  T  -  cu.2). 

4.  Low  stiffness  pinned  link  couplings  which  can  handle  largo  torque 
variations  can  easily  be  designed.  If  a  lower  value  of  6  or  0  is  chosen 
for  f(B,6)  in  equation  38,  thon  n  low  stiffness  coupling  will  result. 
Higher  torque  fluctuations  across  tho  coupling  will  be  possible  without 
causing  Impact  of  the  coupling's  stops  (see  Figure  15). 

5.  Zero  or  low  stiffness  pinned  link  couplings  enn  be  designed  with 
four  basic  equations,  36  -  39. 

0.  Tho  two  advantages  of  link  counterbalancing  nnd  ense  of  design 
analysis  make  the  pinned  link  design  apparently  superior  to  the  elastic 
link  design  for  very  high  speed  applications. 
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PART  II:  EFFECT  OF  ZTS  COUPLINGS 

ON 

CLOSED-LOOP  TORSIONAL  STABILITY 


PRELIMINARY  CONSIDERATIONS 


When  designing  n  dynamic  system,  It  Is  essential  to  be  certain  that  the 
system  will  not  become  unstable.  ThlB  is  especially  true  of  the  closed- 
loop  system,  In  which  a  small  perturbation  can  build  Into  a  destructive 
instability  through  self-excltatton.  A  closed-lot  p  system  is  one  In  which 
the  output  Is  compared  with  a  desired  output,  and  the  resulting  error  Is 
used  as  an  Input  to  the  system.  This  Is  sometimes  called  a  feedback  loop. 

The  complete  dynamics  of  a  modorn  helicopter  is  an  example  of  such  a  closed- 
loop  system.  The  elements  of  the  system  are  the  turbine  engine,  drive 
train,  and  speed  governor.  These  are  shown  In  block  dlsgram  form  in  Fig¬ 
ure  38.  It  has  been  shown  that  to  make  this  system  stable,  either  special 
features  must  be  provided  in  tho  design  of  the  governor  or  the  gain  of  the 
system  must  be  reduced  which  will  Increase  the  steady-state  droop. ®  The 
steady-state  droop  of  a  closed-loop  system  Is  the  difference  between  the 
desired  output  and  the  actual  output  obtained.  Equation  65  gives  the 
steady-state  speed  droop  for  the  closed-loop_hel icopter  system,  where  SSD 
is  the  steady-state  speed  droop  in  percent,  <2e  is  the  steady-state  engine 
speed,  and  Q  is  the  design  spoed  at  which  the  governor  is  set.  Hereafter, 
any  equation  nomenclature  not  defined  In  context  will  be  given  In  the 
List  of  Symbols  at  the  beginning  of  this  report. 

SSD  =  ^  X  100  (65) 

The  governor  modifications  mentioned  above,  introduced  to  ensure  stabil¬ 
ity,  will  decrease  efficiency  or  Increase  expense,  neither  of  which  is 
desirable. 

A  major  objective  of  Part  II  is  to  determine  if  system  stability  can  be 
achieved  by  means  of  effectively  modifying  the  drive  train  without  having 
to  complicate  the  design  of  the  governor.  The  proposed  method  of  modify¬ 
ing  the  drive  trnln  is  to  place  a  zero-torslonal-stif fness-coupl ing  in  the 
drive  shaft,  to  effectively  decouple  the  rotor  inertia  from  the  engine 
inertia.  A  particular  coupling  has  been  designed  in  Part  I  for  tho  oper¬ 
ating  conditions  of  a  helicopter.  It  will  be  attempted  in  Part  II  to 
determine  and  compnro  the  stability  of  the  closed-loop  helicopter  system 
with  and  without  a  zero  torsional  stiffness  coupling  in  the  drive  shaft. 

In  this  wny ,  it  will  be  determined  if  the  use  of  a  ZTS  coupling  is  an 
effective  way  of  maintaining  stability  while  retaining  the  dosirable 
properties  of  a  fast,  simple  governor. 
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In  the  discussion  of  speed  governors,  the  major  emphasis  will  be  plnced 
on  the  direct  fuel  flow  speed  governor.  It  1h  the  simplest  nnd  least 
expensive  speed  governor  available  for  use  in  helicopter  applications. 
Unfortunately,  the  ordinary  helicopter  system  is  unstable  with  a  direct 
fuel  flow  speed  governor.0  If  it  can  be  shown  that  the  helicopter  system 
with  the  ZTS  coupling  is  stable  with  a  direct  fuel  flow  speed  governor, 
the  effectiveness  of  the  coupling  will  have  been  demonstrated. 


TH1.  LINEARIZED  SYSTEM  WITH  AN  IDEAL  COVER  NO  It 

The  helicopter  system  is  inherently  nonlinear.  However,  the  govern¬ 
ing  equations  may  be  linearized,  which  grently  simplifies  the  problem  of 
evaluating  system  stability.  These  equations  are  mndc  linear  by  assuming 
that  the  physical  parameters  of  the  system  remain  within  some  small,  well 
defined,  operating  region.  Linearization  of  equations  will  be  discussed 
in  more  detail  Inter,  and  the  system  governing  equations  will  be  derived. 
The  results  obtained  for  the  linearized  equations  are  directly  applicable 
to  the  case  of  small  perturbations  about  an  operating  point.  They  can 
also  be  helpful  in  getting  nn  Intuitive  "feel"  for  the  more  realistic  case 
of  large  perturbations  nbout  an  operating  point.  More  sped  1  leal  1 y .  it 
will  be  seen  that  perturbations  of  the  gns  producer  and  speed  governor 
system  parameters  are  of  a  sufficiently  small  magnitude  for  realistic  pc 
turbntlons  of  the  drive  train  system  parameters,  that  the  linearized 
governing  equations  for  the  gas  producer  and  speed  governor  are  fairly 
accurate  even  for  the  "large  perturbation"  case.  However,  this  is  unfo. 
tunately  not  the  case  for  the  linearized  drive  train  equations,  especi¬ 
ally  when  the  7.TS  coupling  is  included. 

A  reasonably  accurate  open-loop  model  of  the  helicopter  drive  train  is 
shown  in  Figure  39.  This  system  can  be  represented  by  equations  fili  and  (57. 

..  •  2  ix 
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The  Inertia  lef,  is  the  sum  of  the  power  turbine  inertia,  lL>,  and  the  gear¬ 
box  inertia,  lKtj.  The  only  nonlinear  terms  in  the  two  equations  are  the 
turbine  drag  nnd  rotor  drng.  Actually,  the  engine  torque  is  also  non¬ 
linear  for  the  closed-loop  system  but  thut  will  be  dealt  with  as  one  of 
the  gas  producer  system  parameters.  Under  the  assumption  of  small  changes 
in  the  system  parameters  from  some  steady-state  operating  conditions, 
equations  fid  and  fi7  can  tie  linearized.  This  is  made  less  difficult  by 
first  applying  a  coordinate  transformation.  A  new  coordinate,  6S ,  is 
introduced  nnd  the  coordinate  0^  is  eliminated.  6S  is  defined  as  the 
angle  of  twist  in  the  shaft  and  is  given  by 


HO 


for  this  model.  Using  this  substitution  in  equations  66  and  67  yields 
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Since  Pe  is  now  a  cyclic  coordinate  (i.e.,  does  not  appear  explicitly  in 
the  equations),  these  equations  can  be  rewritten  with  the  following 
substitutions: 


9  =  x 

e  e 


0  =  ID 

c  e 


where  x  is  the  engine  speed,  to  give 
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(68) 
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These  equations  are  linearized  about  a  steady-state  operating  point,  by 
replacing  the  generalized  coordinates  as  follows: 

x  =  x  +  x  8=6+9  T  =  T  ♦  T 

eee  sss  eee 


where  xe  is  the  steady-state  engine  speed,  and  xc  is  the  difference  between 
the  engine  speed  and  its  steady-state  value,  or  the  engine  speed  perturba¬ 
tion.  Making  these  substitutions  gives 
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Since  all  system  variables  are  now  small,  any  above  first  degree  can  be 
ignored  as  being  negligible,  resulting  in 
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By  lettlrg  the  system  perturbations  and  their  time  derivatives  go  to  zero, 
the  above  equations  can  be  made  to  yield  the  steady-state  operating  condi¬ 
tions,  which  also  represent  the  particular  solutions  to  the  differential 
equations.  When  the  system  perturbations  and  their  time  derivatives  are 
set  equal  to  zero,  the  resulting  equations  are 
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These  equations  can  be  solved  for  the  steady-state  operating  conditions. 
Meanwhile,  the  particular  solution  can  now  be  removed  from  the  differen¬ 
tial  equations,  leaving 
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From  these  equations  it  can  be  seen  that  the  linearized  engine  and  rotor 
dampir.g  coefficients  are  represented  by  equations  70  and  71. 
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(70) 
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(71) 

Substituting  these  linearized  damping  coefficients  into  the  differential 
equations  yields  equations  72  and  73  which  represent  the  linearized  heli¬ 
copter  drive  train. 
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Figure  40  is  a  plot  f  the  torque  on  the  engine  due  to  turbine  damping 
as  a  function  of  the  engine  speed.  It  is  of  interest  to  note  that  the 
linearized  engine  damping  coef f icier t  given  by  equation  70  is  the  same  as 
the  slope  of  the  curve  at  the  steady-state  engine  speed.  This  is  reason¬ 
able  since,  in  considering  small  changes  of  torque  and  engine  speed,  it 
is  the  slope  of  this  curve  at  the  operating  point,  rather  than  the  average 
slope,  which  gives  the  linearized  damping  coefficient  of  interest.  It  is 
also  of  Interest  that  these  linearized  damping  coefficients  are  equal  to, 
precisely,  twice  the  average  linearized  damping  coefficient  (average  slope 
of  the  curve) . 


Figure  iO.  Turbine  Damping  Torque  Versus 
Engine  Speed. 
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Again  assuming  small  changes  In  the  system  variables,  the  use  of  a  ZTS 
coupling  would  reduce  the  shaft  stiffness  and  shaft  damping  to  zero, 
effectively  decoupling  the  rotor  Inertia  from  the  engine  inertia.  This  is 
represented  schematically  In  Figure  41.  The  differential  equations  repre¬ 
senting  this  system  are  given  by 

I  '6  +C  02  =  T  -  T 
e  o  e  e  e  c 
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CA 
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where  Tc  is  the  torque  that  the  ZTS  coupling  is  designed 
speed  xe.  Those  equations  are  cyclic  in  both  0  and  0 
substituting  in  6 

to  transmit  at 
Therefore,  by 
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they  can  be  rewritten  as 

•  2 

I  x  ♦  C  u.  =  T  -  T 
e  e  e  e  e  c 

(74) 

Vr  ^  CI.aR  =  Tc 

(75) 

Linearizing  these  equations  and  eliminating  the  particular  solutions, 
using  the  same  method  as  used  above,  yields  equations  76  and  77. 

I  i  *  C'l  =  T 
e  e  T  e  e 

(76) 

Vr  +  Vr  =  0 

(77) 

It  is  noted  that  the  equations  representing  the  system  without  a  ZTS 
coupling  inc’  ide  the  gearbox  inertia  in  with  the  engine  inertia,  while 
the  equations  representing  the  system  with  n  ZTS  coupling  do  not.  This 
is  because  it  is  assumed  that  the  shaft  connecting  the  engine  to  the 
gearbox  is  completely  rigid  (a  reasonable  assumption,  in  comparison  with 
the  mast  shaft  stiffness),  and  that  the  ZTS  coupling  will  be  located  in 
this  shaft.  Therefore,  the  ZTS  coupling  will  also  decouple  the  gearbox 
inertia  from  the  engine  inertin. 

To  obtain  a  feel  for  the  effect  of  the  vnrlous  components  on  the  stability 
of  the  drive  train  system,  the  preceding  equations  will  be  coupled  with 
an  ideal  3peed  governor,  while  the  gas  producer  dynamics  are  ignored. 

In  effect,  this  system  represents  a  single-spool  turbine  engine  with  an 
ideal  governor.  Figure  *42  shows  the  system  of  Figure  39  coupled  with  an 
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Figure  41.  Schematic  of  Helicopter  Drive  Trnin 
With  ZTS  Coupling. 
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ideal  governor.  An  ldenl  speed  governor  is  one  in  which  there  arc  no  lags 
or  other  dynamics  and  the  output  signal  is  directly  proportional  to  the 
input  signal. 

For  the  problem  presently  being  considered,  the  input  signal  is  the  engine 
speed  error  and  the  output  signal  is  the  engine  torque.  The  engine  speed 
error  is  the  difference  between  the  actual  engine  speed  and  the  desired 
engine  speed  to  which  the  speed  governor  lias  been  set.  The  speed  gover¬ 
nor  can  be  represented  by  equation  78. 
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Making  the  same  small  change  substitutions  as  before  yields 
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Solving  thes'j  equations  for  and  X^  produces  equations  83  through  85. 
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This  can  be  written  in  matrix  form  as  shown  in  equation  86. 
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In  this  form  the  system  equations  easily  lend  themselves  to  the  evaluation 
of  stability  by  the  first  method  of  Lyapunov.  When  using  this  method  of 
stability  analysis,  it  is  necessary  to  have  the  differential  equations  in 
the  form  of  equation  87,  where  x  and  x  are  vectors  and  A  is  a  matrix. 

*  --  |x_  (  87  ) 

This  is  the  form  of  equation  86.  The  procedure  is  to  subtract  X  from 
each  element  along  the  main  diagonal  of  the  matrix  A  from  equation  87. 
Following  this  procedure  for  the  matrix  A  of  equation  (86)  results  in  the 
matrix  giver,  by  equation  88, 
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The  result  will  be  u  polynomial  in  X  that  is  the  same  order  as  the  matrix 
.  If  this  polynomial  is  set  equal  to  zero,  the  roots  of  X,  or  eigen¬ 
values  of  the  dynamic  system,  can  be  solved  for.  If  all  of  the  roots  have 
negative  real  parts,  the  system  is  stnble,  since  tin  system  solution  will 
converge  to  a  steady-state  value.  If  any  of  the  roots  have  positive  real 
parts,  the  system  will  be  unstable,  since  the  system  solution  will  grow 
toward  infinity.  If  one,  or  more,  of  the  roots  has  no  real  part  and  the 
rest  of  the  roots  have  negative  real  parts,  the  stability  of  the  system  is 
uncertain.  If  the  above-mentioned  polynomial,  known  as  the  characteristic 
equation,  is  higher  than  second  order,  finding  its  roots  may  be  difficult. 
In  such  a  case  it  may  be  advantageous  to  use  the  "R<>uth-Hurwitz  Criteria,” 
a  method  for  testing  the  coefficients  of  the  characteristic  equation  for 
the  existence  of  roots  with  nonnegative  real  parts  ithout  having  to  actu¬ 
ally  solve  it. 


The  characteristic  equation  derived  from  equation  88  is  given  in  equation 
89. 
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It  appears  that  the  three  eigenvalues  would  be  very  difficult  to  soive 
for.  Since  the  actual  values  of  the  eigenvalues  are  unimportant  for  the 
stability  problem,  it  is  much  more  efficient  to  use  the  " Routh-Hurwi tz 
Criteria"  to  determine  the  signs  of  the  real  parts  of  the  eigenvalues. 


In  order  for  all  of  the  roots  of  a  real  polynomial,  of  the  form  of  equa¬ 
tion  90  to  have  negative  real  parts,  ail  of  the  coefficients  of  the  poly¬ 
nomial  must  be  positive. 
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In  addition, 


the  coefficients 


of  equation  90  must  satisfy  equation  91. 
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From  equation  89  the  coefficients  of  equation  90  for  this  problem  are 
given  by  equations  92  through  94. 
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Since  all  of  *he  physical  system  parameters,  including  the  gain  of  the 
speed  go v erne.  ,  must  be  positive,  then  all  of  the  coefficients  must  also 
be  positive.  Therefore,  the  first  test  for  roots  with  negative  real  parts 
is  satisfied.  Substituting  equations  92  through  91  into  equation  91 
results  in  equation  95,  which  can  be  rearranged  to  produce  equation  96. 
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Since  all  of  the  physical  system  parameters  are  positive,  then  equation  96 
is  satisfied.  Therefore,  the  previously  described  linearized  model  of  the 
helicopter  drive  system  with  an  ideal  governor  nnd  no  gas  producer  dynam¬ 
ics  is  stable. 

By  the  same  method  as  that  just  used,  the  stability  of  the  reduced  system 
represented  by  equations  76  and  77  is  evaluated.  It  can  be  seen  that 
equations  76  and  77  are  decoupled,  and,  since  it  is  the  dynamics  of  the 
engine  speed,  j.e,  that  is  of  primary  interest,  equation  77  can  be  ignored. 
As  in  the  unreduced  problem,  the  system  will  include  an  ideal  speed 
governor  represented  by  equation  78.  This  system  is  illustrated  in  Fig¬ 
ure  -13.  Substituting  equation  79  into  equation  76  results  in  equation  97. 

I  i-  ♦  C_i  =  -Ol  (97) 

e  e  T  e  c 

Using  the  same  method  of  substitution  shown  previously,  equation  97  can 
be  converted  into  the  matrix  form  of  equation  87  shown  in  equation  100, 
using  the  substitutions  of  equations  98  and  99. 


91 


Figure  43.  Schematic  of  Helicopter  Drive  Train 
With  ZTS  Coupling  and  Ideal  Governor 


92 


(99) 


U) 

e 


<C  +G) 

i  xi 

o 


(100) 


The  characteristic  equation  for  the  system  ran  be  found  easily  and  is 
given  by  equation  101. 

(C  +G) 

X  +  — i -  =  0  (101) 
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Since  the  engine  damping  nnd  speed  governor  gain  are  both  positiv. (  the 
one  eigenvalue  of  this  system  of  equotions  is  definitely  negative.  There¬ 
fore,  the  linenrized  model  of  a  helicopter  drive  system  with  an  ideal 
speed  governor  and  no  gas  producer  dynamics  is  stable  whether  or  not  it 
has  a  7.TS  coupling.  However,  in  the  design  and  operation  of  realistic 
helicopter  drive  systems,  stability  problems  are  often  encountered. 
Therefore,  it  must  be  assumed  that  these  instabilities  are  due  to,  or  at 
least  aided  by,  either  the  dynamics  of  realistic  speed  governors  or  real¬ 
istic  gas  producer  dynamics,  or  both.  Consequently,  it  is  necessary  to  go 
to  a  more  sophisticated  nnd  realistic  model  of  the  helicopter  drive  sys¬ 
tem,  including  speed  governor  and  gas  producer  dynnmics,  in  order  to 
determine  the  cause  of  these  system  instabilities  and  the  effectiveness 
of  a  ZT.C  coupling  in  eliminating  them. 


THE  LINEAR! ZED  SYSTEM  WITH  A  DIRECT  FUEL  FLOW  GOVERNOR 


The  equations  for  more  realistic  models  of  the  speed  governor  nnd  gas  pro¬ 
ducer  are  derived  in  this  section.  These  models  are  then  combined  with  the 
previously  described  linearized  model  of  the  drive  trnin  to  form  a  more 
realistic  model  for  the  helicopter  drive  system  undee  discussion. 


The  speed  governor  of  interest  is  a  direct  fuel  flow  speed  governor,  and 
the  conceptual  model  is  shown,  schematically  in  Figure  -H.1  Conceptually, 
the  speed  governor  operates  in  the  following  manner.  The  speed  stick  is 
set  to  a  desired  engine  speed,  f),  thereby  determining  the  value  of  the 
coordinate  S.  The  prr  sure,  p,  on  the  piston  of  mass.  Mg,  is  produced  by 
a  mechanical  or  electrical  transducer  from  the  actual  engine  speed,  xe. 

This  pressure  is  assumed  to  be  directly  proportional  to  the  actual  engine 
speed.  This  is  represented,  analytically,  in  equation  102. 


p  K  x 
j  e 


(102) 


When  operating  at  steady  state,  the  piston  will  be  in  the  position  to 
allow  the  correct  amount  of  fuel  flow,  such  that  the  engine  speed  resulting 
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from  that  fuel  flow  will  create  a  pressure,  p,  which  will  result  in  a 
force  on  the  piston  equal  in  magnitude  to  the  force  exerted  on  the  piston 
by  the  spring  of  stiffness,  Kg.  If  the  engine  speed  suddenly  decreases, 
possibly  due  to  an  increased  load,  then  the  force  on  the  piston  due  to  the 
pressure,  p,  will  become  less  than  the  spring  force  on  the  piston.  The 
piston  will  move  to  the  right,  decreasing  the  spring  force,  until  the 
forces  on  the  piston  are  balanced.  This  will  result  in  an  inc*  jased  fuel 
flow,  increasing  the  engine  speed,  and  subsequently  forcing  the  piston 
to  move  back  to  the  left.  This  damped  oscillation  of  the  piston  will  con¬ 
tinue  for  a  short  time  until  the  piston  settles  down  into  a  new  steady- 
state  position.  The  damping  of  the  piston  is  assumed  to  be  purely  viscous. 
This  appears  to  be  an  accurate  assumption.*”  Also,  an  increase  in  the 
engine  speed  or  a  pilot-induced  change  in  the  set  speed  will  produce 
dynamic*  of  the  piston  similar  to  those  described  above. 

The  differential  equation  of  motion  for  the  piston  is  derived  by  means  of 
Newton's  law,  represented  by  equation  103. 

U  MX  (103) 

x 

Under  the  assumptions  that  the  damping  is  viscous  and  the  spring  is  l'~*»nr, 
the  equation  of  motion.  104,  will  be  linear. 

M  X  •  C  X  ••  K  X  =  K  A  x  -  K  iO  +  K  D  (104) 
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The  speed  governor  can  be  physically  designed  so  that  equation  105  is 
satisf ied. 

K  a  K  A 
g  i  P 

If  this  is  done,  then  equation  (104)  will  reduce  to  equation  106, 

M  X  *  C  X  *  K  X  K  faU  -fi)  ♦  D] 
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When  relating  parts  of  the  helicopter  drive  system,  it  is  much  more  con¬ 
venient  to  use  the  fuel  flow,  Wf,  as  a  system  parameter,  than  to  use  the 
piston  displacement.  Xg.  Assuming  that  the  fuel  flow  orifice  is  rectangu¬ 
lar.  that  the  fuel  flow  reduction  near  the  boundary  of  the  orifice  is 
negligible,  and  that  the  fuel  supply  pressure  is  constant,  the  fuel  flow 
can  be  related  to  the  piston  displacement  by  equation  107.  Under  the 
above-listed  assumptions,  the  partial  derivative  of  the  fuel  flow  with 
respect  to  the  piston  displacement  is  a  constant. 

dWf 
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(105) 
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This  constant  is  negative,  since  an  increase  in  Xg  will  cause  a  decrease 
in  the  fuel  flow.  Making  the  substitution  of  equation  108  into  equa¬ 
tion  107,  noting  that  Q  is  positive,  and  solving  the  resulting  equation 
for  Xg  will  give  equation  109. 
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The  first  two  time  derivatives  of 
110  and  111. 


Equations  109  through  111  are  subs 
equation  112.  This  equation  is  in 
piston  displacement. 


are  taken  as  shown  in  equations 
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ituted  into  equation  106,  resulting  in 
terms  of  the  fuel  flow  rather  than  the 
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A  schematic  representation  of  the  gas  producer  and  power  turbine  is  given 
in  Figure  45.  3  This  is  known  as  a  two-spool  or  free  turbine  engine. 

The  fuel  flow  is  delivered  from  the  speed  governor  to  the  combustor.  The 
fuel  is  burned  in  the  combustor,  which  creates  a  driving  torque  on  the  gas 
producer,  as  well  as  on  the  power  turbine.  The  gas  producer  creates  a 
driving  torque  on  the  power  turbine  as  well  as  on  the  compressor.  There¬ 
fore,  when  the  rate  of  fuel  flow  is  suddenly  changed  it  has  a  d'reet, 
rapid  effect  on  the  engine  torque.  However,  it  also  has  an  indirect, 
delayed  effect  through  the  speed  of  the  gas  producer.  That  is,  a  change 
in  fuel  flow  causes  a  change  in  the  gas  producer  speed  that  is  not  imme¬ 
diate,  due  to  the  inertia  of  the  gas  producer.  This  change  in  gas  pro¬ 
ducer  speed  then  causes  a  change  in  the  engine  torque  created  by  the  gas 
producer.  This  delayed  reaction  of  engine  torque  to  fuel  flow  is  signif¬ 
icant  in  considerations  of  helicopter  drive  system  instability. 


The  differential  equation  of  motion  for  the  gas  producer  is  derived  by 
means  of  Newton's  law,  as  given  in  equation  113. 

ET  =  I  i  (113) 

g  g  g 
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Under  the  assumption  that  the  relative  changes  in  I  icl  flow,  gas  producer 
speed,  and  engine  torque  are  fairlj  small,  equatioi.  Ill  is  derived  from 
equation  113. 
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Each  of  the  partial  derivatives  actually  represents  a  set  of  curves  such 
as  those  shown  in  Figure  46."*  For  example,  equation  115  represents  the 
graph  of  Figure  *16.  If  the  region  of  operation  is  relatively  small,  the 
curves  can  be  assumed  tc  approximate  straight  line-  through  the  operating 
point,  thereby  making  the  partial  derivatives  cons  int  values  for  that 
operating  point. 
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The  same  argument  can  be  used  to  find  constants  to  replace  the  partial 
derivatives  in  equation  114.  It  is  likely  that  these  linearized  equations 
are  sufficiently  accurate  to  give  reliable  results  for  realistic  pertur¬ 
bations  about  the  steady-state  operating  conditions. 

Figure  47  is  an  illustration  of  a  single-spool  or  fixed  turbine  engine. 

The  two-spool  engine  has  the  advantage  of  allowing  the  gas  producer  spool 
to  spin  free  of  the  power  turbine  if  a  large  enough  load  is  suddenly 
applied  to  the  power  turbine  to  stop  it,  or  at  lenst  slow  it  down  signif¬ 
icantly  very  quickly.  If  such  a  load  were  applied  to  a  single  spool 
engine,  the  compressor  would  also  be  forced  to  slow  down  or  stop,  and  the 
potential  for  damage  to  the  turbine  engine  would  be  much  greater.  However, 
a  single-spool  engine  has  the  advantage  of  being  less  dynamically  complex 
and  more  inherently  stable.  A  curve  representing  the  relationship  between 
engine  torque  and  fuel  flow  for  a  fixed  turbine  engine  would  look  like  one 
member  of  the  family  of  curves  shown  in  Figure  16.  Under  the  small  pertur¬ 
bation  assumption,  the  engine  torque  for  a  one-spool  turbine  engine  would 
be  given  by 


dT 

fc  ^  *f 

where  the  partial  derivative  can  be  considered  a  constant. 

In  order  to  carry  out  a  Bode  plot  stability  analysis  of  a  helicopter 
drive  system,  it  is  necessary  to  convert  the  system  differential  equa¬ 
tions  into  Laplace  transfer  functions.  The  transfer  function  for  the 
speed  governor  is  derived  from  equation  112.  Equation  117  is  substi¬ 
tuted  into  equation  112  and  the  Laplace  transform  is  taken  of  the 
resulting  equation,  giving  equation  118. 
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Figure  17.  One-Spool  Turbine  Kim i no 
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The  speed  error,  C,  is  a  variable  rather  than  a  constant,  and  therefore 
its  Laplace  transform  is  the  Lnplace  function  E(s).  Equation  118  is 
rearranged  to  form  equation  (119).  Also,  since  the  system  model  is  based 
on  the  assumption  of  small  changes,  the  transfer  function  is  written  in 
terms  of  small  changes  in  the  system  parameters. 
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Dividing  the  numerator  and  denominator  of  the  transfer  function  by  Kg 
results  in  equation  120. 
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(120) 


The  denominator  of  this  transfer  function  can  be  fnc  .red,  and  with  the 
substitutions  of  equations  121  through  123,  equation  12-1  will  be  the  result¬ 
ing  speed  governor  transfer  function. 
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(121) 


(122) 


(123) 


(121) 


Equation  124  is  the  same  transfer  function  as  given  in  Reference  5. 

The  transfer  function  for  the  complete  turbine  engine  is  derived 
from  equations  114  and  115.  The  Laplace  transform  of  those  two  equations 
is  taken  after  the  substitutions  of  equations  125  through  128  are  made. 
Equations  129  and  130  result  from  equations  114  and  115  respectively. 
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( 126) 
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It  Is  advantageous  to  eliminate  the  gas  producer  speed  Laplace  function. 
Wg(s),  from  the  equations.  This  is  accomplished  by  solving  equation  129 
for  Wg(s)  and  substituting  the  result  into  equation  130  to  give  equation 
131. 


ATe(s)  =  AWf(s)K2  4  AWf(s)  .  •  (131) 

g'  1 

Both  sides  of  the  equation  are  divided  by  AWf(s)  and  the  terms  on  the 
right  side  of  the  equation  are  regrouped  to  form  equation  132.  Also,  since 
the  equation  is  based  on  the  assumption  of  small  changes,  the  system  var¬ 
iables  are  written  as  small  changes  in  system  parameters. 
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The  transfer  function  is  further  rearranged  and  the  substitutions  of 
equations  133  through  135  are  made  in  order  to  derive  equation  136. 
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This  equation  is  the  Laplace  transfer  function  representing  the  two-spool 
turbine  engine  and  gas  producer.  It  is  the  same  as  the  gas  producer  trans¬ 
fer  function  given  in  References  5  and  4. 

The  transfer  function  for  the  one-spool  turbine  engine  is  derived  from 
equation  116.  The  Laplace  transform  is  taken  of  this  equation,  and  equa¬ 
tion  127  is  substituted  in  to  give  equation  137. 

AT  (s)  =  K-AW.(s)  (137) 

e  2  f 


This  leads  to  the  transfer  function  given  by  equation  138,  where  K2  is  not 
numerically  equal  to  the  found  in  equation  133),  since  they  are  found 
from  different  curves. 
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K2  from  equation  137  is  assumed  to  be  of  the  order  of  magnitude  of  Kg  from 
equation  136.  That  is,  the  steady-state  gain  of  a  one-spool  turbine  engine 
is  assumed  to  be  about  the  same  as  that  for  a  two-spool  turbine  engine. 

Equations  72  and  73  are  used  to  derive  the  transfer  function  for  the  linear¬ 
ized  drive  train  without  a  ZTS  coupling.  The  Laplace  transform  of  both 
equations  is  taken,  resulting  in  equations  139  and  140  with  the  A’s  repre¬ 
senting  small  changes  in  the  system  parameters,  from  their  steady-state 
values,  the  same  as  the  's  in  time  domain. 

I  AW  (s)S  +  OIW  (s)  ♦  C  A?  ( s) S  +  K^e  (s)  =  AT  (s)  (139) 
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re  s  D  e  s 

-  C  A6  (s)S  -  K  -A 9  (s)  0  (140) 
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Equation  140  is  solved  for  A9<?(s)  to  give  equation  141. 
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This  is  substituted  into  equation  139  in  order  to  eliminate  A0  (s)  from 
the  transfer  function,  resulting  in  equation  112.  This  equation  is  con¬ 
solidated  to  give  equation  1-13. 
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This  equation  is  divided  by  AWe(s)  and  inverted.  Then,  the  substitutions 
of  equations  144  through  119  are  made  in  order  to  get  the  transfer  func¬ 
tion  into  the  more  tractable  form  of  equation  150. 
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Equation  150  is  the  Laplace  transfer  function  which  represents  a  linear¬ 
ized  helicopter  drive  train  without  a  ZTS  coupling. 


The  only  Laplace  transfer  function  left  to  derive  for  the  Bode  plot 
stability  analysis  is  the  one  representing  the  linearized  helicopter 
drive  train  with  a  ZTS  coupling;  that  is,  the  system  with  the  decoupled 
rotor  represented  by  equations  76  and  77.  Again,  equation  77  has  no 
effect  on  the  engine  speed  and  is  therefore  of  no  interest.  This  leaves 
Just  equation  76,  which  is  reproduced  as  equation  151. 
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The  Laplace  transform  of  this  equation  is  given  by  equation  152 

I  SAW  (s)  ♦  C  AW  (s)  =  AT  (s)  (152) 
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Solving  this  equation  for  gives  equation  153. 
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The  substitutions  of  equations  151  and  155  into  equation  153  result  in 
equation  156,  which  is  the  Laplace  transfer  function  which  represents 
n  linearized  helicopter  drive  train  with  a  ZTS  coupling. 
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Equation  153  nay  also  be  derived  by  allowing  the  torsional  •  tiffness,  F  , 
and  torsional  damping,  C  ,  of  the  nast  shaft  to  go  to  z  o  1  replacing 
Igg  with  Ie  in  equation  143,  which  is.  in  effect,  what  the  ZTS  coupling 
is  supposed  to  do. 

The  total  open-loop  transfer  function  for  a  helicopter  drive  system 
without  a  ZTS  coupling  is  given  by  the  product  of  equations  124,  136,  and 
150,  as  shown  in  equation  157. 
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The  total  open-loop  transfer  function  for  the  system  with  a  ZTS  coupling, 
the  product  of  equations  124,  136,  and  156,  is  given  in  equation  158. 
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A  flow  chart  of  the  closed-loop  helicopter  drive  system  is  shown  in 
Figure  48.  Since  it  has  unity  feedback,  the  closed-loop  transfer  func¬ 
tion  can  be  found  from  the  open-loop  transfer  function  using  equation  15J. 
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The  open-loop  transfer  function  represents  the  ratio  of  the  engine  speed 
to  the  speed  error.  This  is  substituted  into  equation  159  to  form  equa¬ 
tion  160. 
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The  Laplace  transform  of  equation  117  is  given  in  equation  161,  which  is 
substituted  into  equation  160,  resulting  in  equation  162, 
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To  find  the  steady-state  speed  droop  of  the  system,  it  is  necessary  to 
find  the  steady-state  value  for  the  engine  speed.  The  transfer  functions 
Yy(s)  V2(s)  ,  and  Y2(^)  cm  be  used  in  their  present  form.  However,  the 
engine  and  rotor  damping  coefficients  of  Y.^(s)  and  Y^(s)  were  derived  for 
the  purpose  of  giving  a  change  in  damping  torques.  In  order  to  give  the 
actual  damping  torques  at  steady-state  operating  conditions,  it  is  neces¬ 
sary  to  use  the  average  damping  coefficients,  which  are  equal  to  exactly 
half  of  the  linearized  damping  coefficients.  With  this  substitution  in 
Y3(s)  ,  the  A's  are  removed  from  equation  162.  To  find  the  steady-state 
droop,  the  set  speed  is  considered  to  be  a  step  input,  and  then  the  engine 
speed  can  be  represented  by  equation  163. 

W  (s)  =  Y  (s)  -  (163) 

e  c  s 

In  this  equation,  Yc(s)  is  the  closed-loop  transfer  function  which  will 
give  the  actual  system  parameters,  with  the  intent  of  finding  the  steady- 
state  system  parameters. 

The  steady-state  value  of  a  Laplace  transfer  function  can  he  found  by 
using  the  final  value  theorem,  given  in  equation  161. 
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Using  the  final  value  theorem,  the  steady-state  value  of  the  engine  speed 
is  found  in  equation  165. 
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Substituting  this  equation  into  equation  65  for  the  steady-state  speed 
droop  results  in  equation  166. 

SSD  =  (1  -  Y  (0 ) )  y  100  (166) 
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For  the  helicopter  drive  system  without  the  ZTS  coupling,  Y  (0)  is  given 
by  equation  167. 
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Therefore,  the  steady-state  speed  droop  of  the  system  without  the  ZTS 
coupling  is  given  by  equation  168. 

SSD°  (rrahnr) x  100  <168) 

p  B  L 

For  the  system  with  a  ZTS  coupling,  the  transfer  function  Y3<s)  was  deter¬ 
mined  for  small  changes  such  that  a  constant  torque,  Tc ,  was  dropped  from 
the  governing  equations.  Consequently,  that  transfer  function  cannot  be 
used  to  determine  the  steady-state  system  parameters,  since  they  would  all 
be  off  by  some  factor  of  that  constant  torque.  Therefore,  the  above  pro¬ 
cedure  cannot  be  used  to  find  the  steady-state  speed  droop  for  a  helicop¬ 
ter  drive  system  with  a  ZTS  coupling. 

The  steady-state  speed  droop  can  also  be  derived  directly  from  the  differ¬ 
ential  equations  of  motion.  This  method  of  derivation  is  applicable 
whether  the  system  has  a  ZTS  coupling  or  not.  An  attempt  is  made  to  find 
the  stationary  values,  or  steady-state  solutions,  of  the  system  parameters. 
This  is  done  by  setting  the  highest  time  derivatives  in  the  system  differ¬ 
ential  equations  equal  to  zero,  and  solving  for  the  system  parameters. 

For  a  helicopter  drive  system  without  a  ZTS  coupling,  equations  68,  69, 

112,  lid,  and  115,  at  steady  state,  become  equations  169  through  173 
respectively.  The  engine  and  rotor  damping  terms  are  replaced  by  the 
average  linearized  damping  terms. 
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Also,  the  values  of  the  partial  derivatives  of  equations  172  and  173 

are  assumed  to  be  constant  all  the  way  from  zero  to  steady-state  engine 

speed,  for  the  purpose  of  determining  the  steady-state  system  parameters 

only.  Since  §  is  equal  to  zero  at  steady  state,  equations  169  and 
*  s 

170  reduce  to  equations  17*1  and  175. 
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Equations  171  through  173  r.re  manipulated  in  order  to  get  equations  for 
the  steady-state  values  for  the  system  parameters  as  functions  of  the 
steady-state  engine  speed.  This  results  in  equations  176  through  178. 
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Using  the  substitutions  of  equations  117,  121,  125  through  128,  ami 
133),  these  reduce  to  equations  179  through  181. 
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Equations  174  and  175  are  added  together  and  equation  181  is  used  to 

substitute  for  T  ,  giving  equation  182. 
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Using  equation  144  and  resubstituting  for  the  speed  error  from  equation 
117  results  in  equation  183. 
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Solving  for  the  steady-state  engine  speed  gives  equation  181,  and  insert¬ 
ing  this  into  equation  65  gives  equation  185. 
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Equations  168  and  185  are  identical. 

The  above  procedure  is  repeated  for  a  helicopter  drive  system  with  a  ZTS 
coupling.  Equations  7<1  and  75  with  the  average  linearized  engine  and 
rotor  damping  terms  become  equations  186  and  187  at  steady  state. 
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Equations  186  and  187  are  added  and  equation  181  is  used  to  substitute 
for  T  ,  giving  equation  188,  considering  that  the  steady-state  engine 
speed  is  equal  to  the  steady-state  rotor  speed. 
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Using  equation  111  and  resubstituting  for  the  speed  error  from  equa¬ 
tion  117  results  in  equation  189. 


x  2K  K  K  (."  -  i  )  (189) 
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Solving  this  equation  for  the  steady-state  engine  speed  gives  equation 
190,  and  inserting  this  into  equation  65  gives  equation  191. 
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It  can  be  seen  from  equations  185  and  191  that  the  steady-state  speed 
droop  for  this  helicopter  drive  system  is  the  same  whether  or  not  it  has 
a  ZTS  coupling. 
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BODE  ANALYSIS  OF  THE  LINEAR  SYSTEM 


A  Bode  plot  of  a  frequency  transfer  function  is  a  plot  of  the  gain  and 
phase  of  the  transfer  function  versus  frequency.  The  frequency  transfer 
functions  of  interest  are  derived  by  taking  the  open-loop  transfer  func¬ 
tions  of  equations  157  and  158  and  making  the  substitution  of  equation  192. 

s  =  (192) 


This  results  in  equations  193  and  191. 
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The  gain  of  the  frequency  transfer  func*  on  F(jx),  given  by  equation 
195,  is  expressed  in  decibels. 

GAIN  (x)  =  20  log  |f(,1x)|  (195) 


The  phase  of  the  same  frequency  transfer  function  is  given  by  equation 
196  and  is  expressed  in  degrees. 


PHASE  <x) 


/  F(Jx) 


where 

-180  degrees  ^  /  K(jx)  -  180  degrees  (196) 


A  Bode  plot  is  actually  comprised  of  two  plots.  One  plot  is  of  the  gain 
of  the  frequency  transfer  function  versus  x  and  the  other  plot  is  of  the 
phase  of  the  frequency  transfer  function,  versus  x .  The  gain  margin  for 
a  frequency  transfer  function  is  defined  ns  the  gain  increase  needed  to 
drive  the  system  unstable.  ^  The  gain  margin  is  calculated  ns  the  nega¬ 
tive  of  the  gain  at  the  point  of  phase  crossover.  The  phase  crossover 
is  where  the  phase  plot  crosses  over  180  degrees,  or  minus  180  degrees. 
The  phase  margin  is  equal  to  the  phase  shifted  180  degrees  at  the  point 
of  gain  crossover.  The  gain  crossover  is  where  the  gain  plot  crosses 
over  zero  decibels.  [f  the  gain  or  phase  margin  of  a  system  is  negative, 
then  that  system  is  unstable.  If  the  gain  and  phase  margins  of  a  system 
are  both  positive,  then  that  system  will  converge  to  some  steady-state 
operating  condition  following  a  perturbation.  However,  the  lower  the 
gain  or  phase  margin,  the  longer  the  system  will  take  to  stabilize. 
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For  the  type  of  helicopter  drive  system  under  discussion,  a  minimum  gain 
margin  of  4  to  6  decibels  and  phase  margin  of  20  tc  30  degrees  is  neces¬ 
sary  to  achieve  an  acceptable  transient  response. ^ 


The  helicopter  drive  system  of  interest  is  the  helicopter  analyzed  in 
Reference  5  with  a  normal  operating  engine  power  output  of  5270  horsepower 
at  a  speed  of  13,820  revolutions  per  minute  or  1447  radians  per  second. 

It  is  assumed  that  at  steady-state  operating  conditions,  all  of  this  power 
is  dissipated  as  aerodynamic  rotor  damping.  Torsional  aerodynamic  damping 
is  generally  assumed  to  be  of  the  form  of  equation  197. 
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The  relationship  of  torque  to  Horsepower  and  speed  is  given  in  equation  198, 
with  speed  in  radians  per  set  and  ar.d  torque  in  foot-pounds. 


T 


L 


550  hp 


•a 

e 


(198) 


The  value  for  C,  can  be  found  by  substituting  the  torque  from  equation 
198  into  equation  197  and  solving  for  the  damping  coefficient.  The  result 
is  given  as  equation  199. 

Cj  =  9.564  >  Hf1  ft-lb  (rad  per  sec)2  (199) 

From  equation  71  the  linearized  rotor  damping  coefficient  is  calculated 
and  given  in  eqiation  200. 

CR  =  2.768  ft-lb  rad  sec  (200) 

Since  the  model  for  the  helicopter  power  plant  is  borrowed  from  Refer¬ 
ence  5,  then  it  is  reasonable  to  use  the  helicopter  drive  train  and 
speed  governor  models  from  the  same  reference.  Tlx.*  system  parameters 
for  the  direct  fuel  flow  speed  governor  of  Reference  5  are  given  as 
follows : 


T  =  0.02  sec 

T  =  0.03  sec 
2 

The  gas  producer  and  turbine  engine  from  the  same  reference  have  the  fol¬ 
lowing  system  parameters: 


K 


b 


0. 914 


ft-lb 
lb  hr 


T  =  0 . 1869  sec 

n 
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0 . 3  sec 


T 

e 


K  =  0.57 
2  lb/hr 


K_  =  3.2563 
3  rnd  sec 


The  drive  train  system  parameters,  as  referenced  to  the  engine  design 
speed,  from  Reference  5  are  as  follows: 


I  =0.5 
e 


ft-lb 

r 

rad  sec 


I 


0.1944 


Kb 

I  =  3.888 


ft-lb 

r 

rad/ sec 
ft-lb 


rad  'sec 


k  ,  =  21,417  — 
TM  rad 


Kg  =  93.0  to  837.3 


ft-lb 

rad 


The  torquemeter  stiffness,  Kj.^,  is  at  least  25  times  the  mast  stiffness, 
Kq.  Therefore,  the  torquemeter  can  be  assumed  to  be  rigid  as  compared  to 
the  mast.  The  range  of  mast  stiffness  listed  will  give  the  required  crit¬ 
ical  frequencies  of  the  model  of  2  to  6  cycles  per  second.  In  order  to 
reproduce  the  stability  results  reported  in  Reference  5  for  the  direct 
fuel  flow  speed  governor,  it  is  necessary  to  use  a  mast  stiffness  of 
about  800  foot-pounds  per  radian,  corresponding  to  a  critical  frequency 
of  about  5.86  cycles  per  second. 


The  shaft  damping,  C  ,  is  assumed  to  be  about  5  percent  of  the  critical 
shaft  damping  for  the  torsional  two-degree-of -freedom  system.  The  crit¬ 
ical  shaft  damping  for  such  a  system  is  represented  by  equation  201. 


C 


SCR 


(201) 


Using  the  system  parameter  values  given  above,  the  critical  shaft  damping 
is  found.  Five  percent  of  that  value  Rives  the  shaft  damping  used;  see 
equation  202. 
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2.24  f t-lb/rad 'sec 


(202) 


C  = 
s 


The  torque  necessary  to  overcome  the  turbine  damping  is  assumed  to  be 
about  40  percent  of  the  load  torque  of  equation  198.  The  turbine  damping 
is  assumed  to  be  aerodynamic,  as  shown  in  equation  203. 


T  =  C  t2  (203) 

T  e  e 

Using  T  ,  equal  to  40  percent  of  T  ,  and  equation  203  ,  is  calculated 
and  given  in  equation  204. 


C  =  3.82  x  10-4  ft-lb  (rad  sec)2  (204) 

e 

The  linearized  turbine  damping  coefficient  is  calculated  from  equation  70 
and  is  given  in  equation  205. 

C  =  1.106  ft-lb  rad  sec  (205) 

T 


The  gain  of  the  drive  train,  K  ,  is  a  function  of  the  linearized  turbine 
and  rotor  damping,  and  is  founa  in  equation  206. 

K  =  - - - =  0.258  rad  sec  ft-lb  (206) 

L  CT  *  CD 


The  only  system  parameter  of  this  helicopter  drive  system  model  which 
is  not  yet  specified  is  the  governor  gain,  K  .  Equation  168,  for  the 
steady-state  speed  droop,  can  be  solved  for  the  governor  gain  resulting 
in  equation  207. 


100-SSD 
p  "  2K  K  ( SSD) 


(207) 


With  a  5-percent  steady-state 
helicopter  drive  system5)  and 
culated  and  given  in  equation 


speed  droop  (standard  for  this  type  of 

values  of  K  and  K.  from  above,  K  is  cal- 
D  L  P 


K  =  40.27  lb  hr  (rad  sec)  (208) 

P 


The  transfer  functions  used  in  Reference  5  for  the  direct  fuel  flow  specu 
governor  and  the  gas  producer  and  turbine  engine  are  identical  to  those 
derived  in  this  report.  However,  the  transfer  function  for  the  drive 
train  used  in  Reference  5,  given  by  equation  209,  is  a  simplification  of 
the  one  derived  in  this  report,  given  by  equation  150. 
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This  transfer  function  is  derived  in  Appendix  V.  To  obtain  equation  209 
from  equation  150,  first  shaft  damping  is  set  equal  to  zero.  Then,  the 
denominator  is  approximately  factored  and  the  following  substitutions  are 

made: 
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2  2 
Cl  +C  I 
D  eg  T  R 


2vK.InI  (I  +In) 
°  R  eg  eg  R 


3  2 


The  error  introduced  by  the  approximate  factoring  is  negligible  However, 
the  error  introduced  by  ignoring  the  shaft  damping  can  be  significant. 

In  fact,  for  the  value  of  damping  given  by  equation  202,  the  error  of 
stability  margin  for  the  helicopter  drive  system  of  Reference  5  is  con¬ 
siderable.  This  will  be  demonstrated  by  a  comparison  of  the  gain  and 
phase  margins  of  the  helicopter  drive  system  model  of  Reference  5  with 
zero  shaft  damping  and  with  the  value  of  shaft  damping  calculated  above. 
The  gain  and  phase  margins  are  found  from  Bode  plots  which  are  made  by 
inserting  incremental  values  of  0.01  to  10.0  cycles  per  second  into  the 
frequency  transfer  functions  representing  the  helicopter  drive  system 
model  with 
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1.  Zero  shaft  damping — from  equation  193 

2.  Nonzero  shaft  damping — from  equation  193 

3.  A  ZTS  coupling  in  the  drive  shaft--from  equation  191 

The  resulting  Bode  plots  are  given  in  Figures  49,  50,  and  51  respectively. 
The  gain  and  phase  margins  arc  given  in  Table  I.  It  can  be  seen  that 
the  addition  of  shaft  damping  makes  the  unstable  system  of  Reference  5 
stable.  However,  it  still  does  not  meet  the  minimum  phase  margin  criteria, 
in  order  to  give  an  acceptable  transient  response.  It  can  also  be  seen 
that  the  addition  of  a  ZTS  coupling  gives  a  system  stability  which  is 
better  than  the  zero  shaft  damping  case,  but  is  not  significantly  better 
than  the  case  with  shaft  damping.  Furthermore,  the-  system  with  the  ZTS 
coupling  does  not  meet  the  minimum  phase  margin  criteria.  Therefore,  ii 
appears,  at  least  for  the  linearized  model,  that  the  addition  of  a  ZTS 
coupling  to  the  helicopter  drive  system  of  Reference  5  does  not  increase 
the  stability  of  that  system  enough  to  be  acceptable. 

The  points  for  the  Bode  plots  are  calculated  and  plotted  by  the  computer 
program  listed  in  Appendix  VI.  The  same  computer  program  also  makes  Bode 
plots  for  the  two  helicopter  drive  systems  using  two,  other,  more  sophisti¬ 
cated,  speed  governors.  These  two  speed  governors  are  typical  examples  of 
the  results  of  adding  additional  dynnmics  to  a  direct  fuel  flow  speed 
governor  in  order  to  achieve  stability.  Unfortunately,  they  also  add 
additional  cost  and  complexity  to  the  direct  fuel  flow  speed  governor. 

One  of  these  other  speed  governors  is  a  direct  fuel  flow  speed  governor 
with  lagged  gain  reset.  It  is  basically  similar  to  a  direct  fuel  flow 
speed  governor  with  the  addition  of  a  feedback  loop  from  the  speed  gover¬ 
nor  output  to  vary  its  speed  setting.  This  feedback  loop  acts  as  a  reset 
function  through  a  time  lag.°  Block  diagrams  for  a  direct  fuel  flow 
speed  governor  and  a  direct  fuel  flow  speed  governor  with  lagged  gain 
reset  are  reproduced  from  Reference  5  in  Figures  52  and  53,  respectively. 

A  transfer  function  for  the  second  type  of  speed  governor  is  found  from 
Figure  53  and  is  given  in  equation  210. 


£Wf(s) 

TeTsT 


TTdTT  (Tr 


(Tis+1)1nr  s+ii(t2s+1) 

r 


Y^s) 


(210) 


Bode  plots  of  the  helicopter  system  without  a  ZTS  coupling,  using  a  direct 
fuel  flow  speed  governor  with  lagged  gain  reset,  are  given  in  Figure  54 
for  zero  shaft  damping  and  in  Figure  55  for  nonzero  shaft  damping.  To  get 
these  Bode  plots,  values  of  0.2  second  for  the  reset  time  constant,  t  ,  and 
0.5  (dimensionless)  for  the  reset  gain,  ,  are  used.0  Also,  a  value  of 
20.14  pounds  per  hour,  per  radian  per  second,  is  used  for  the  speed 
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Figure  50.  Bode  Plots — Direct  Fuel  Flow  Governor, 
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Figure  51.  Bode  Plots--Direct  Fuel  Flow  Governor  Kith  ZTS  Coupling  (Decoupled  Rotor  Model). 
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TABLE  I  -  Continued 
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Figure*  53.  IJirect  Fuel  Flow  S|>«*e*d  Governor  With  I. Jigged 
Gain  He*sc*t  Block  Di;igrnm. 
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Figure  54.  Bode  Plots--Direct  Fuel  Flow/Lagged  Gain  Reset,  C  =0.0 


(SIOqtDOQj  UIBJ 


126 


Figure  55.  Bod?  Plots--Direct  Fuel  Flow/Lagged  Gain  Reset,  C  *  2.24 


governor  proportional  gain,  Kt>  to  retain  a  5-percent  steady-state  speed 
droop.  Then  equation  210  is  used  in  place  of  equation  124  in  finding  the 
open-loop  transfer  function  of  the  helicopter  drive  system,  used  to  get 
the  frequency  transfer  function,  which  subsequently  yields  the  Bode  plots 
for  the  system,  using  the  computer  program  in  Appendix  VI.  The  gain  and 
phase  margins  determined  from  all  of  these  Bode  plots  are  given  in  Table  II. 
It  can  be  seen  that  this  system  is  stable  even  with  zero  shaft  damping, 
although  it  does  not  meet  the  minimum  gain  or  phase  margin  criteria  for 
a  practical  system.  However,  with  '.he  nonzero  value  of  shaft  damping 
used,  this  system  is  not  only  stabi  .  but  it  easily  meets  and  exceeds 
those  minimum  gain  and  phase  margin  criteria.  This  demonstrates  the  value 
of  adding  complexity  to  the  speed  governor. 

The  other,  even  more  complex,  type  of  speed  governor  presented  is  a  direct 
fuel  flow  compressor  discharge  pressure  speed  governor  with  lagged  gain 
reset.  This  type  of  speed  governor  is  very  similar  to  a  direct  fuel  flow 
speed  governor  with  lagged  gain  reset,  with  the  further  addition  of  a 
feedback  signal  from  the  compressor  discharge  pressure  to  be  used  as  an 
additional  criterion  for  determining  the  fuel  flow.  ® 

A  block  diagram  for  a  direct  fuel  flow  compressor  discharge  pressure  speed 
governor  with  lagged  gain  reset  is  reprr  (uced  from  Reference  5  in  Figure  56. 
A  transfer  function  for  this  type  of  go.,  .‘nor  is  found  in  the  same  refer¬ 
ence  and  is  given  by  equation  211. 
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Bode  plots  are  made  of  the  helicopter  drive  system  without  a  ZTS  coupling, 
with  a  direct  fuel  flow  compressor  discharge  pressure  speed  governor  with 
lagged  gain  reset,  with  zero  and  nonzero  shaft  damping,  using  the  same 
method  as  outlined  above,  and  are  given  in  Figires  57  and  58  respectively. 
To  get  this  Bode  plot,  in  addition  to  the  values  given  above,  a  value  of 
0.773  second  is  used  for  *j.J  The  gain  and  phase  margins  from  these  Bode 
plots  are  given  in  Table  II.  It  can  be  seen  from  the  positive  gain  and 
phase  margins  that  this  system  is  not  only  stable,  but  it  meets  the  min¬ 
imum  gain  and  phase  margin  criteria,  with  either  zero  or  nonzero  shaft 
damping.  Also,  it  is  more  stable  than  either  of  the  speed  governors  pre¬ 
viously  presented,  but  is  also  much  more  complex. 

It  was  previously  shown  that  the  helicopter  drive  system  with  an  ideal 
speed  governor  and  no  gas  producer  dynamics  (one-spool  turbine  engine)  is 
stable,  whether  or  not  it  includes  a  ZTS  coupling.  This  is  substantiated 
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by  a  Bode  analysis.  Replacing  the  direct  fuel  flow  speed  governor  with  an 
ideal  governor  and  using  a  one-spool  turbine  engine,  whose  transfer  func¬ 
tion  is  given  by  equation  138,  equations  193  and  194  become  equations  212 
and  213. 
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Assuming  that  the  value  of  K0  for  a  one-spool  turbine  engine  is  the  same 
as  the  value  of  K.  for  a  two-spool  turbine  engine,  these  equations  can  be 
derived  from  equations  193  and  19  1  by  setting  all  governor  and  gas  pro¬ 
ducer  time  constants  equal  to  zero.  Bode  plots  are  constructed  for  this 
system  for  the  following  cases: 


1.  Zero  shaft  damping — from  equation  212 

2.  Nonzero  shaft  damping — from  equation  212 

3.  ZTS  coupling  in  drive  shaft — from  equation  213. 


The  resulting  Bode  plots  are  shown  in  Figures  59,  60,  and  61,  respectively. 
It  can  be  seen  that  there  is  no  phase  crossover.  Therefore,  this  system 
is  stable  for  any  positive  value  of  speed  governor  gain.  This  agrees  with 
the  previous  characteristic  equation  analysis  results.  Therefore,  it 
appears  that  the  principal  cause  of  the  closed-loop  instability  is  not  in 
the  drive  train,  but  is  rather  in  the  governor  dynamics  or  in  the  gas  pro¬ 
ducer  dynamics  or  in  a  combination  of  both.  However,  if  the  drive  train 
is  not  the  principal  cause  of  the  system  instability,  it  does  not  necessar¬ 
ily  follow  that  it  is  not  a  contributing  factor. 


In  an  attempt  to  determine  the  relative  effect  of  the  speed  governor  and 
gas  producer  dynamics  on  the  system  stability,  a  Bode  analysis  is  done  on 
the  helicopter  drive  system  with  an  ideal  speed  governor  and  a  two-spool 
turbine  engine.  That  is.  the  speed  governor  dynamics  are  eliminated,  but 
the  gas  producer  dynamics  are  included.  The  frequency  transfer  functions 
for  this  system,  corresponding  to  equations  193  and  194,  are  given  by 
equations  213  and  215  respectively . 

K.>KnKi,Tn(Jjl)4l,|Ai(jA,2+Ao(Ji)*1l 

^(.U)  '  - 5 - ^ -  (214) 

I  '  ,<j  Mill!  (ji)  4B0(ji  )”+B  (Jr)4l) 


F2(Jl  > 


KPKnWJi,  +  11 

[T  <J  .  )4ll[B.(.)l>4ll 
e  1 


(215) 
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Figure  60.  Bode  Plots--Ideal  Governor,  Single-Spool  Turbine,  Engine, 
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Figure  61.  Bode  Plots--Ideal  Governor,  Single  Spoo 


Bode  plots  are  made  for  the  same  three  cases  as  listed  above,  except  that 
the  frequency  transfer  functions  given  by  equations  214  and  215  are  used 
instead  of  those  given  by  equations  212  and  213.  The  resulting  Bode  plots 
are  given  in  Figures  62,  63,  and  64.  Once  again  there  in  no  phase  cross¬ 
over;  therefore,  this  system  is  stable  for  any  positive  gain  of  the  ideal 
speed  governor. 

For  comparison,  the  gas  producer  dynamics  can  be  eliminated  by  using  a  one- 
spool  turbine  engine  in  the  model  ,  whose  transfer  function  is  given  by 
equation  138,  while  the  dynamics  of  the  speed  governor  are  included.  Bode 
plots  are  made  for  the  helicopter  drive  system  model  with  a  direct  fuel 
flow  speed  governor  and  a  one-spool  turbine  engine  for  the  same  three  cases 
as  were  considered  above.  The  necessary  frequency  transfer  functions  are 
derived  by  substituting  equation  138  for  equation  136  in  the  derivation  of 
equations  193  and  19-1  in  order  to  obtain  equations  216  and  217. 

K  K, ,K  (A  (Jr)2  ♦  A0jr+U 

F  (Jr.)  =  - P  - - ^ - 5 -  (216) 

r1<j*)nir9(ji.)+:il  B  <ji>\b0(Ja)  +b  (.i«')+i| 


Wlc 

F  (Jr)  =  - ------ -  (217) 

[T1(Jx)*l]['2(jx)UHB|(ji)+lJ 

Using  tlie  value  for  K9  as  the  same  'is  the  value  for  KQ  of  the  two-spool 
turbine  engine,  the  Bode  plots  resulting  for  those  three  cases  are  given 
in  Figures  65,  66,  and  67  and  the  gain  and  phase  margins  are  given  in 
Table  II.  It  appears  that  this  system  is  even  less  stable  then  the  system 
with  both  speed  governor  and  gas  producer  dynamics.  Therefore,  it  seems 
that  the  speed  governor  lags  are  the  principal  causes  of  the  closed-loop 
system  instability. 

In  general,  the  addition  of  a  7.TS  coupling  to  all  of  these  systems  greatly 
simplifies  the  Bode  plots  of  the  systems.  That  is,  it  reduces  the  number 
of  gain  crossovers  and  eliminates  the  peak  in  the  gain  plot,  which  repre¬ 
sents  a  critical  frequency  of  the  system.  These  are  definite  advantages 
to  the  use  of  a  7.TS  coupling.  However,  from  the  Bode  analysis,  it  also 
appears  to  be  generally  true  that  the  addition  of  an  undamped  7TS  coupling 
to  an  unstable,  or  insufficiently  stable,  system  does  not,  in  itself, 
improve  the  stability  of  the  system  enough  to  be  practically  useful. 


TRANSIKNT  S0UT10N  FOB  HIE  I,  I  NT  AH  SYSTEM 

The  transient  solution  of  the  engine  speed,  in  time  domain,  is  calculated 
from  the  linear  Laplace  transfer  functions.  Equation  162  is  solved  for 
the  variation  in  the  engine  speed  from  a  steady-state  value  to  give  equa¬ 
tion  218. 
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Figure  63.  Bode  Plots-- Ideal  Governor,  Two-Spool  Turbine  Engine  C  ■  2  24  ft~lb  • 
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1  igure  (i4 .  Bode  Plots--Idcal  C.ovcrnor,  Two-Spool  Turbine  lingine  Kith  ITS  Coupling  (recoupled 
Itotcr  Mode  1 )  . 
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(218) 


AW  (s)  =  Y  (s)  A0(«) 
e  c 


The  set  speed  is  assumed  to  be  a  step  function  of  magnitude  Cl,  making  the 
variation  of  the  set  speed  an  impulse.  Therefore,  AD(s)  can  be  replaced 
by  AfJ.  giving  the  ratio  of  the  variation  of  the  engine  speed  from  its 
steady-state  value  to  the  change  in  governor  set  speed  ns  equation  219. 


AW  (s) 

IT 


Y  (s) 
c 


(219) 


Since  it  is  the  variation  of  the  engine  speed  from  the  steady-state  value 
that  is  being  solved  for,  rather  than  the  engine  speed  itself,  the  solu¬ 
tion  will  be  harmonic  about  zero.  In  order  to  fim'  the  solution  of  Ai 
in  the  time  domain,  it  is  necessary  to  take  the  inverse  Laplace  transform 
of  the  closed-loop  transfer  function.  Equations  157  and  159  are  combined 
to  find  the  closed-loop  transfer  function  for  the  system  without  a  ZTS  cou¬ 
pling,  resulting  in  equation  220. 

K  K  K  (T  S+l )  (A,  S^+A-Sf  1 ) 

Y  (s)  =  ,B - L_2 - 

C  K  K  K  (T  S+1)(A  S  +A,S*1) 

p  B  L  d  1  i 


+ 


(T  Sfl)(T  Sfl)(T  sa)(B-S3vB0S2+B_S+l) 
1  £,  C  1  <5  ^ 


(220) 


Using  the  substitutions  in  equations  219  and  221  through  227,  this  equa¬ 
tion  cnn  be  rewritten  ns  equation  228. 


D1  = 


TlVeBl 


'221 ) 


°2  =  V2W(TlVTlVYe,Bl 


(222) 


°3  =  TlT2W(TlYTlVT2VY(YYVBl 

D4  =  TlT2Te+(TlT2+TlTetT2Te)V(T1<T24Te)Vni+KpKBK.,TnAl 

(T^yAj) 


5  ~  1  2  1  e  2  e  i  2  e'  3  2 


°6  ^  (YYV+B3+KpKBKl/YV 


D  =  K  K  K,  +  1 
7  p  B  L 


(223) 

(221) 

(225) 

(226) 

(227) 


AW  (s)  K  K  K  (1  S»1)(A.S  »A  Stl) 

e  _ pBL  1) _ 1  2 _ 

L‘  (D,  S6+D_S^+  D„S^4  D.S^f  U  S2*  IJ  St  D  ) 

1  2  3  4  5  6  7 


( 22K) 
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In  order  to  find  the  inverse  Laplace  transform  of  this  function,  the  denom¬ 
inator  must  first  f>e  factored  and  the  function  must  be  broken  up  by  partial 
fraction  expansion.  This  has  been  done,  with  the  same  system  parameters 
used  previously  for  the  Bode  plots,  by  the  computer  program  in  Appendix  VII 
This  is  done  for  both  zero  shaft  damping  and  for  the  same  nonzero  value  of 
shnft  damping  that  was  used  for  the  Bode  analyses.  For  zero  shaft  damping, 
equation  229  results  and  for  nonzero  shaft  damping,  equation  230  results. 


:w  (S) 
e 

ir. 


|0. 138  (0.02  4  .j  0.082)  (0.02  -  j  0.082) 

1 )  “  (S-0.055-J  13.9)  ”  (S-O.OSS+j  13.9) 


0.180  (0.011  -  J  0.062)  (0.0‘11  +  j  0.062)1  ,  -2 

(S+ 11 . 2)  (S  4  1. 87  -  j  1. 66)  (s+'l.  87  +  j  4. 66)  J 


(229) 


“"e(S)  0.139  (0. 022  -4  j  0.083)  (0.  022  -  j  0. 083) 

Lt  L(S»68.  5)  "  (Sil.73  -  J  -M.  1)  "(8+1.73  4  J  11 . 1 ) 

0.175  (0.01  .  J  0.062)  (0.01  +  J  0.062)  1  ,n-2 

(S+ll.l)  4  (8+  1.85  +  j  1.66) 4  (S+1.85  -  J  1.66)  J  Xl 


(230) 


The  inverse  Laplace  transform  of  equation  229  is  given  by  equation  231, 
while  the  inverse  Laplace  transform  of  equation  230  is  given  by  equation 
232. 


x  (t) 
e 


=  0. 001 38e 


-68.lt 


0  05  5 1 

0 . 00168e  ‘  sin  (13. 9t  -  0.21) 


-  0.0018e 


-11. 2t 


..j  ^7 

0.001l9e  sin  (1.66t  +  0.58) 


(231) 


— -  =  0.00139e”68  5t 

* 

- 1 7  3 1 

4  0.00172c  sin  (ll.lt  -  0.26) 


0.00175c 


-ll.lt 


-•I  85t 

+  0.00117c  '  sin  (  1 . 66 1  (  0.57) 


(232) 


These  functions  are  plotted  in  Figures  68  and  69  For  the  system  with 
a  /.TS  coupling,  the  closed-loop  transfer  function  is  found  to  be  repre¬ 
sented  by  equat ion  233. 
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Figure  68.  Linear  Time  Response,  C.  =  0.0 
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Figure  69.  Linear  Time  Response,  Cs  =  2.24 
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Y  (s)  _ WLC1^11 _ 

C  '  KpKBKLC(TBS4-l)+  (T^l)  (T2Sfl)  (TeSf  1)  (B^Sfl) 


The  substitutions  of  equations  219  and  234  through  238  reduce  the  above 
equation  down  to  equation  239. 


E1  * 

TlVeB4 

(234) 

E2  3 

TlT2Te+T1T2B4+TlTeVT2TeB4 

(235) 

E3 

TlT2+TlTe+TlB1+T2Te+T2B4+TeB4 

(236) 

E4  3 

Vb^cVVYV^i 

(237) 

E5  3 

kpkbklc+1 

(238) 

4  Vs’  Wlc<tbs>1) 

(239) 

on 

4  3  2 

E,  S  +  E_S  +  E.S  +  E .S+  E_ 

1  2  3  4  5 

Similar  to  the  above  situation,  the  denominator  of  this  function  must  be 
factored,  and  the  resulting  equation  must  be  broken  up  by  partial  fraction 
expansion.  This  has  also  been  done  by  the  computer  program  in  Appendix  VII, 
resulting  in  equation  240. 


AWe(S)  !"  0.115  (-0.137Sf8.25)  0.0011  “1  -2 

00  =L(Sf74.7)+(^4i34)2+(31_8)2-(Sf5.46)J  * 


(240) 


The  inverse  Laplace  transform  of  the  previous  function  is  given  by  equa¬ 
tion  241. 


U)  (t) 
e 

n 


[0.145e-‘M'7t  +  0.28e“4'3t  sin 


(31. 8t  -0.421) 


-0.00110e“5-'16tl  x  10~2 


(241) 


This  function  is  plotted  in  Figure  70. 

The  zero  ordinate  lines  in  Figures  68  through  70  represent  the  new  values 
of  steady-state  engine  speed  after  the  change  in  set  speed.  Therefore,  the 
initial  peak  of  each  curve  represents  the  overshoot  of  the  system.  The 
entire  curve  shows  the  transient  response  of  each  system  as  i  t  approaches 


146 


0015 


cr> 


(GSOTUOTSUOUITa)  - 5“ 

U)  « 

V 


1 17 


Figure  70.  Linear  Time  Response  With  ZTS  Coupling 


its  steady-state  values.  From  Figure  68  it  can  be  seen  that  the  engine 
speed  for  the  case  with  zero  shaft  damping  diverges,  and  is,  therefore, 
completely  unstable.  However,  from  Figures  69  and  70,  the  other  two  cases 
are  seen  to  have  engine  speeds  which  converge,  and  so  arc  at  leust  margin¬ 
ally  stable.  Comparing  Figures  69  and  70,  it  is  apparent  that  the  helicop¬ 
ter  system  with  a  7TS  coupling  has  a  greater  overshoot,  but  also  a  shorter 
settling-down  period,  than  the  system  without  a  ZTS  coupling.  The  greater 
overshoot  is  reasonable  since  the  inertia  of  the  linearized  system  is 
reduced  by  the  use  of  n  ZTS  coupling.  That  is,  when  the  set  speed  is  first 
changed,  the  initial  change  in  the  fuel  flow  will  be  the  same,  and  the 
Initial  change  in  engine  torque  will  be  the  same,  whether  or  not  the  system 
has  a  ZTS  coupling.  Therefore,  the  system  with  the  smaller  inertia  will 
have  the  greater  initial  acceleration,  or  deceleration,  and  will  likely 
have  the  greater  overshoot.  The  shorter  settling-down  period  is  also  rea¬ 
sonable.  This  is  becnusc,  for  the  linearized  system  with  a  ZTS  coupling, 
only  the  oscillation  of  the  engine  inertia  must  be  damped  out.  However, 
for  the  linearized  system  without  a  ZTS  coupling,  oscillation  must  be 
damped  out  of  both  the  engine  and  rotor  inertins,  which  are.  in  addition, 
out-of-phase  oscillations.  It  would  seem  that  a  shorter  settling  time  is 
a  more  important  advantage,  in  considering  transient  response,  than  a 
smaller  overshoot  is.  Therefore,  according  to  the  linearized  analysis, 
adding  a  ZTS  coupling  to  a  helicopter  drive  system  will  improve  the  time 
response  of  the  system. 

For  comparison,  as  was  done  with  the  Bode  analysis,  the  time  response  of 
the  engine  speed  can  also  be  found  for  the  helicopter  drive  system  without 
speed  governor  dynamics,  without  gas  generator  dynamics,  and  without  both. 

First,  for  the  helicopter  drive  system  with  an  ideal  governor  and  a  one- 
spool  turbine  engine,  the  frequency  transfer  functions  of  equations  212  and 
213  are  converted  to  Laplace  transfer  funct'ons  by  replacing  each  jx  with 
an  S.  Then  the  closed-loop  transfer  functions  are  derived,  with  the  aid  of 
equation  159,  to  become  equations  2-12  and  2-13,  respectively. 


Y  (s) 
c 


K  K  K  (A  S~* A>  1  ) 
|)  D  I.  1  2 


3  2 

B,  S  ♦(B.M.K  K  K  )  S  +  ( !)„♦  A  K  K  K  )S+K  K  K  *  1 
1  2  1  p  B  L  3  2  p  B  1.  p  B  I. 


Y  (s)  = 
c 


K  K  K 
•>  B  LC 


Wd'S-c1  1 


(212) 


(213) 


The  inverse  Laplace  transforms  for  these  transfer  functions  for  the  same 
three  cases  considered  above,  again  us‘ng  the  computer  program  in  Appen¬ 
dix  VII,  are  given  by  equations  241,  245,  and  246. 
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(244) 


-  (1.05  y  10_8)e“15-4t 


_q  _2l  2t 

— -  =  (8.76  xlO  )e  sin  (16. 5t  +  1.54) 


,1 


„  ori  -8  -16, 7t 

-  (1.32  x  10  )e 


(245) 


—75  8t 

— -  =  -  (1.08  x  10  )  e  75'8t 


(246) 


These  functions  are  plotted  in  Figures  71,  72,  nnd  73.  All  three  plots 
illustrate  a  very  quick  response.  It  is  much  quicker  than  the  response 
of  the  system  with  speed  governor  and  gns  producer  dynamics.  However, 
there  is  a  marked  similarity  between  the  response  of  the  two  systems. 
That  is,  both  of  these  systems,  when  including  a  ZTS  coupling,  have  a 
larger  Initial  variation  from  steady  state,  but  converge  to  steady  state 
much  quicker. 


The  next  model  used  Includes  gns  producer  dynamics,  in  the  form  of  n  two- 
spool  turbine  engine,  nnd  an  ideal  speed  governor.  The  closed-loop  Laplace 
transfer  functions  for  this  model  are  derived  from  the  open-loop  frequency 
transfer  functions  of  equations  214  and  215  in  the  same  manner  as  described 
above.  The  resulting  equations,  247  and  248,  are  Inverse  transformed  by 
the  computer  progrnm  in  Appendix  VII  to  give  equations  249,  250,  and  251 
for  the  same  three  cases  used  previously. 


Y  <S> 
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K  K„K,  (T0S*l)(AtS  +A  S>-1 ) 
p  B  L  B  1  2 


K)}KbK  j  ( T  0S+ 1 )  ( Al  S2+  A2S+ 1 )  +  ( T  1 )  ( S34  B2S2+  B.jS+1 ) 


(247) 


Y  (S)  = 
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K  K_K  (T  ♦ 1 ) 
p  D  LC  B 


T  B.S  +(T  +B.+K  KJ,  T  )  S+K  K„K,04l 
e4  e  4  p  BT  LC  B  p  B  LC 
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4  O. 041 He 


-4 , 62t 
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Figure  73.  Linear  Time  Response — Ideal  Governor,  One-Spool  Turbine  Eneine 
With  ZTS  CouDlin*. 
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00381c"5  5t 
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These  three  functions  are  plotted  in  Figures  71,  75,  and  76,  respectively. 
It  can  be  seen  from  the  first  two  graphs  that  this  system  without  a  ZTS 
coupling  has  a  rather  large  initial  variation  from  steady  state,  but  it 
converges  very  quickly.  However,  the  time  response  of  this  system  with 
a  ZTS  coupling,  ns  shown  in  Figure  76,  is  much  improved.  Therefore,  it 
appears  that  a  ZTS  coupling  does  help  to  slgnif icently  reduce  the  tran¬ 
sient  oscillations  in  a  helicopter  drive  system  with  gas  producer  dynam¬ 
ics,  but  with  no  speed  governor  dynamics. 


The  next  model  used  Includes  speed  governor  dynamics,  in  the  form  of  n 
direct  fuel  flow  speed  governor,  but  no  gas  producer  dynamics,  modeled  ns 
n  one-spool  turbine  engine.  The  closed-loop  l.nplace  transfer  functions 
for  this  model  are  derived  from  the  open-loop  frequency  transfer  functions 
of  equations  216  and  217  in  the  same  manner  ns  described  above.  The  re¬ 
sulting  equations,  252  and  253,  arc  inverse  transformed  by  the  computer 
progrnm  in  Appendix  VII  to  give  equations  254,  255,  and  256. 
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This  is  for  the  same  three  cases  used  previously. 
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These  three  functions  are  plotted  in  Figures  77,  78,  and  79.  It  can  be 
seen  from  Figure  77  that  this  system  with  zero  shaft  damping  diverges, 
and  is,  therefore,  completely  unstable.  This  system  with  nonzero  shaft 
damping,  from  Figure  78,  is  stable,  but  it  converges  rather  slowly.  With 
a  ZTS  coupling,  this  system  has  both  a  smaller  initial  variation  from  the 
steady  state  and  a  more  rapid  convergence  to  the  steady  state.  Therefore, 
it  appears  that  a  ZTS  coupling  certainly  improves  the  time  response  of 
this  system).  As  compared  to  the  previously  analyzed  ideal  governor,  two- 
spool  turbine  engine  drive  system  model,  this  system  has  a  smaller  Initial 
variation  from  the  steady  state,  but  it  also  has  a  much  slower  convergence 
to  the  steady  state. 


In  general,  it  appears  that  the  results  and  comparisons  of  the  several 
systems  anil  cases,  deduced  from  the  time  response  nnalyses,  agree  very 
well  with  the  stability  information  found  from  the  Bode  analyses. 
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In  Part  I,  a  ZTS  coupling  is  designed  for  high  speed,  high  power  applica¬ 
tions,  such  as  a  helicopter  drive  system.  This  ZTS  coupling  is  illustrated 
in  Figures  80  through  82.  The  drive  train  with  the  ZTS  coupling  is  shown 
in  Figure  83.  The  governing  equations  for  this  drive  train  are  given  by 
equations  257  through  262. 
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A  coordinate  transformation  is  made,  once  again  employing  9., 
of  twist  in  the  shaft.  For  the  helicopter  drive  train  with’  a 
8  is  given  by 

S  6  9  -  9,  -  6. 

s  e  II 

Substituting  6  Into  equations  257,  2G3  through  265,  261, and 
ing  8|{  and  1  ts‘  der ivnt i  ves ,  and  replacing  9  by  and  8  by 
in  equations  266  through  271. 
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Taking  the  time  derivatives  of  equations  267,  269,  and  271  and  inserting 
them  into  equations  266,  268,  and  270,  respectively,  results  in  equations 
272  through  27  1. 
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In  order  (o  linear l/e  these  equations,  a  Taylor's  series  expansion  c  . 
each  equation  is  taken  about  the  steady-state  values  of  engine  speed, 
coupling  angle,  and  shaft  twist  angle.  Assuming  that  the  variations  of 
these  parameters  from  their  steady-state  variables  are  relatively  small, 
all  terms  of  the  Taylor's  series  above  the  first  degree  can  be  Ignored. 
An  ordinary  Taylor's  series  expansion  is  given  by  equation  27S.  s 
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f(x)  f(u)  +  (x-«)f2(a)  +  .  .  .4  f(n) 


Expanding  this  to  cover  more  thnn  tine  variable  and  eliminating  all  terms 
above  first  order  results  in  equation  27G, 
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Using  equation  276  with  j  ,  8  and  6  and  their  time  derivatives  as  the 
variables  and  substituting 
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where  i  Ik  the  steady -state  engine  speed  and  x  is  the  variation  of  tlu* 
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engine  speed  from  the  steady  state,  the  result  Is  equation  277. 
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Using  equation  27 i  *o  llnearl/e  equat ions  272  through  271  results 
in  equations  278  through  280. 
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Kquntlons  70  and  71  arc  substituted  into  equations  278  through  280. 

In  addition,  the  following  substitutions  are  made  for  constant  coeffi¬ 
cients  ; 
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Ths  resulting  differential  equations  are  given  by  equations  281  through 

283. 
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Setting  ell  system  variables  and  their  time  derivatives  equal  to  zero 
yields  equations  284  through  286. 
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These  equations  represent  the  steady-state  solutions  of  equations  281 
through  283,  It  Is  the  transient  solutions  that  are  of  primnry  Interest, 
since  they  determine  stability.  The  stoady-st ate  solutions  enn  be  removed 
from  equations  281  through  283  by  subtracting  equations  28-1  through  281 
respectively,  yielding  equations  287  through  289. 
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These  are  the  linearized  equations  for  a  helicopter  drive  train  with  the 
/.IS  coupling  designed  in  Part  I. 


A  frequency  transfer  function  is  derived  from  these  equations  to  be  used 
in  a  llode  analysis.  The  l.iplnce  transform  is  taken  of  equations  287 
through  289  to  give  equations  290  through  292. 
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Equation  292  is  solved  for  A@g(S)  giving 
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which  is  substituted  into  equations  290  and  291  in  order  to  eliminate 
A“s(s)  from  the  equntions.  The  resulting  equations  are  given  by  equa¬ 
tions  293  and  29-1. 
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which  is  substituted  into  equation  293  in  order  to  eliminate  t  (S)  from 
the  equation.  In  addition,  the  following  substitutions  are  made: 
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The  resulting  equation  Is  given  by  equation  295 
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This  e(|uat  ion  is  converted  into  a  I.aplncc  transfer  function  with  the  aid 
of  the  following  substitutions: 
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The  transfer  function  is  given  hy  equation  296. 
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The  correspond  In#  frequency  transfer  function  is  found  hy  substiMitlng 
,jw  for  s  in  equntion  296. 

In  order  to  determine  the  values  of  the  coefficients  of  equations  296,  it 
Is  necessary  to  find  the  steady-state  values  for  the  system  parameters. 
These  are  found  by  solving  the  differential  equations  of  the  system  for 
the  stationary  solutions  of  the  system  parameters.  This  is  done  by  setting 
the  highest  time  derivatives  in  the  system  differential  equations  equal 
to  zero,  and  solving  for  the  stationary  solutions.  Kquntlons  179  through 
181  for  stendy-state  fuel  flow,  gas  producer  spied  and  engine  torque, 
respectively,  can  be  npplicd  to  the  nonlinear  case,  because  the  differ¬ 
ential  equations  for  the  speed  governor  and  gas  producer  arc  approximated 
by  assuming  that  the  average  values  of  the  coefficients  in  those  differ¬ 
ential  equations  are  approximately  equal  to  the  values  of  those  coeffi¬ 
cients  used  for  small  changes  about  steady  state.  This  does  not  compromise 
the  major  objective  of  determining  the  effect  of  the  ZTS  coupling  on  the 
closed-loop  system.  To  find  the  stationary  solutions  for  the  engine  and 
rotor  speeds,  equations  272  and  271  are  reduced  as  described  above  to 
produce  equations  297  and  29H 
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Adding  these  equations  and  substituting  equation  181 
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Substituting  into  this  equation  for  the  speed  error  and  rearranging  the 
terms  gives 
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This  equntion  is  solved  for  u  using  the  quadratic  equation  to  obtain 
equation  299,  the  steady-state  engine  speed. 
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Using  the  value  of  Kj,  determined  for  a  steady-state  speed  droop  of 
5  percent  for  the  linearized  system,  and  using  the  vnlues  given  previ¬ 
ously  for  the  rest  of  the  parameters  in  equation  299,  -c  is  solved  for. 
The  steady-state  engine  speed  for  the  nonlinear  case  is  equal  to  1378 
radians  per  second  for  a  set  speed  of  1H7  radians  per  second.  Substi¬ 
tuting  these  values  into  equntion  63  gives  a  steady-state  speed  droop  of 
about  *1.8  percent  for  the  nonlinear  system  with  the  same  gain  that  gives 
the  linear  system  a  5-percent  steady-state  speed  droop.  To  find  the 
stationary  solution  for  the  angle  of  the  coupling,  9,  equation  273  is 
reduced  as  described  above  to  yield 


Equation  298  is  subtracted  from  this  equation,  and  the  result  is  divided 
-2 

by  a  to  give  equntion  300. 
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The  ZTS  coupling  is  designed  for  this  helicopter  drive  system,  to  give 
the  following  coupling  parameters  (see  Pinned  Link  Design  Analysis,  Part  I): 

N  2  m  0  2385  lh-mass 

r  3  in  4.  2 .61  in. 
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Using  the  previously  introduced  value  of  rotor  damping,  the  stationary 
solution  for  the  coupling  angle  Is  found  from  equation  300  to  be  1.61  radi¬ 
ans.  The  steady-state  speed  error  is  found  from  equation  117  nnd  the 
values  for  C.  and  w  calculated  above.  The  steady-stato  fuel  flow  for  this 
set  speed  is  found  from  equation  179  to  be  approximately  2760  pounds  per 
hour.  Kquatlon  180  gives  the  steady-state  value  of  the  gas  producer  speed 
as  about  294  radians  per  second.  From  equation  181,  the  steady-state 
value  for  the  torque  applied  to  the  power  turbine  is  found  to  be  about 
2540  foot-pounds. 

To  find  the  steady-state  values  of  engine  and  rotor  speed  for  the  same 
helicopter  drive  system  model ,  but  without  a  ZTS  coupling,  equations  68 
and  69  are  reduced  as  above  to  give 
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These  equations  are  identical  to  equations  297  and  298.  Therefore,  the 
stationary  solutions  of  engine  and  rotor  speed  for  this  helicopter  drive 
system  nre  the  same  whether  or  not  it  has  a  ZTS  coupling.  It  follows  that 
the  steady-state  values  for  engine  torque,  fuel  flow,  and  gas  producer 
speed  are  also  independent  of  the  use  of  this  ZTS  coupling. 

Without  a  special  damper,  the  viscous  damping  coefficient  C,.  would  l>  • 
negligible.  However,  in  the  hope  of  increasing  the  stability  ol  th 
tern,  this  damping  coefficient  may  he  artificially  increased,  througl 
use  of  an  elastomeric  damper.  Therefore,  in  further  stability  analyse, 
of  this  ZTS  coupling,  Ct,  will  he  used  as  zero  and  as  50  percent  of  crit 
leal  damping,  for  comparison.  Since  the  helicopter  system  with  the  ZTS 
coupling  is  no  longer  a  simple  two-degree-of -f reedom  torsional  spring- 
mass  system,  the  determination  of  th?  critical  damping  is  greatly  compl i 
ented.  However,  it  is  assumed  that  the  critical  damping  of  this  system 
is  of  the  same  order  of  magnitude  as  the  critical  damping  of  the  heli¬ 
copter  system  without  a  ZTS  coupling.  Previously,  this  was  found  to  be 
nbout  45  foot  oounds  per  radian  per  second.  Therefore,  it  is  safe  to  use 
a  value  of  10  foot-pounds  per  radian  per  second  for  coupling  damping 
without  fear  of  exceeding  50  percent  of  critical  damping.  Since  the  cou¬ 
pling  damping  in  produced  by  an  elastomeric  dnmper,  the  exact  value  of 
damping  is  arbitrary,  because  the  damper  can  be  designed  to  create  what¬ 
ever  damping  is  desired.  Part  of  the  intention  of  this  investigation  is 
to  determine  the  effect  of  coupling  damping  on  the  stability  of  the  system, 
rather  than  to  determine  the  effect  of  a  specific  amount  of  coupling 
damping  on  the  stability  of  the  system. 

There  »s  now  sufficient  information  to  determine  the  value  of  substitutions 
Fj  through  Fg  and  all  of  the  coefficients  of  equation  296.  The  values  of 


Ei  nnd  E2  turn  out  to  bo  much  smaller  than  the  values  of  the  Inertias  It,, 
Ije,  nnd  Ijp.  Rased  on  this,  equations  287  through  2H9  enn  be  reduced  to 

equations  301  through  303. 
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It  enn  be  seen  from  equation  302  that  tills  coupling  docs  not  truly  have 
zero  torsional  stiffness  when  dovinted  from  the  equilibrium  solution,  but 
rather  has  a  low-torsional  stiffness  with  a  spring  constant  given  by  F“1  j 

The  frequency  transfer  function  for  the  helicopter  drive  train  with  this 
ZTS  coupling,  derived  from  equation  296,  is  combined  with  the  frequency 
transfer  functions  for  the  direct  fuel  flow  speed  governor  nnd  two-spool 
turbine  engine,  to  give  equation  301. 
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The  coefficients  A,  and  n  can  be  factored  out  of  the  numerator  and 
5  6 

denominator,  respectively,  resulting  in  equation  30.r>,  where  K(  is  the  same 
ns  defined  earlier. 


(305) 


This  transfer  function  results  in  the  same  steady-state  speed  droop  as 
the  similar  transfer  functions  in  previous  sections.  Therefor1*  to  achieve 
n  5-percent  steady-state  speed  droop,  tin*  saint*  value  ol  speed  governor  gain 
is  used  as  was  used  before 

A  Rode  analysis  of  the  frequency  transfer  function  ol  equation  305  was  con¬ 
ducted.  This  was  done  for  values  of  zero  and  2.21  for  shaft  damping,  and 
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zero  and  10.0  for  coupling  damping,  all  in  units  of  foot-pounds  per  radian 
per  second.  All  four  combinations  were  analyzed  and  the  resulting  Bode  plots 
are  given  in  Figures  84  through  87.  The  gain  and  phase  margins  found  from 
these  Bode  plots  are  given  In  Table  II.  It  appears  that  all  four  of  these 
cases  are  not  only  stable,  but  they  even  meet  the  minimum  margin  require¬ 
ments  for  practical  stability  as  given  in  Reference  5.  Furthermore,  it 
appears  that,  although  shaft  damping  seems  to  have  little  effect  on  system 
stability,  Introduction  of  coupling  damping  makes  the  system  even  more 
stable.  This  is  a  very  desirable  situation,  because  the  ZTS  coupling  can  be 
designed  with  whatever  amount  of  damping  is  necessary  to  meet  realistic  sta¬ 
bility  requirements.  Therefore,  from  this  linearized  stability  analysis, 
it  appears  that  the  ZTS  coupling  designed  in  Part  I  can  be  used  to  make  the 
helicopter  drive  system  stable  enough,  with  a  direct  fuel  flow  speed  gover¬ 
nor,  to  be  practically  usable  without  necessitating  redesigning  the  speed 
governor. 


TABLE  II. 

BODE  RESULTS  FOR  LINEAR 

COUPLING  DESIGN 

Shaft  Damping 
(ft-lb  rad  per  sec) 

Coupling  Damping 
(ft-lb  rad  per  sec) 

Gain  Margin 
(db) 

Phase  Margin 
(deg) 

0.0 

0.0 

7.1 

26.5 

2.24 

0.0 

7.1 

26.6 

0.0 

10.0 

7.2 

64.0 

2.  24 

10.0 

8.2 

63.0 
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Figure  87.  Bode  Plots--Direct  Fuel  Flow  Governor  Kith  ZTS  Coupling  (Linearized  Design  Model) 


NUMERICAL  SIMULATION  OF  THE  NONLINEAR  CLOSED-LOOP  SYSTEM 


In  general,  there  are  no  known  methods  to  solve  the  differential  equations 
of  motion  for  the  nonl lnearlzed  system  analytically.  It  is  possible  to 
solve  the  differential  equations  numerically  with  the  r.id  of  a  digital 
computer.  This  is  also  known  as  simulating  the  dynamic  system.  The  numer¬ 
ical  routine  involves  marching  the  differential  equations  through  small 
increments  in  time.  For  the  purposes  of  the  numerical  simulation,  the  dif¬ 
ferential  equations  of  the  drive  trnin  are  used  in  their  cnnonicnl  form. 


The  drive  train  differential  equations  for  the  system  without  a  ZTS  coupling 
are  converted  to  canonical  form  from  equations  66  and  67  to  form  equations 
306  through  309. 


(306) 


(307) 


(308) 


(309) 


For  the  drive  train  with  the  ZTS  coupling,  shown  in  Figure  83,  the  govern¬ 
ing  differential  equations  are  given  by  equations  257  through  262.  The 
differential  equation  for  a  direct  fuel  flow  speed  governor  is  given  by 
equation  112.  After  making  the  substitution  of  equation  121,  the  result 
is  equation 


M  *  *  C  W# 

g  1  g  f 


K  W, 
K  f 


K  K  Or  ) 

IP  e 


(310) 


The  differential  equations  for  a  two-spool  turbine  engine  are  derived 
from  equations  114  and  115,  under  the  assumption  that  the  average  values 
of  the  partial  derivatives  are  approximately  equal  to  the  values  used 
under  the  small  changes  assumptions.  The  substitutions  of  equations  125 
through  128  are  made,  resulting  in  equations  311  and  312. 
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A  computer  program  Is  developed  in  Part  I  to  numerically  Integrate  equa¬ 
tions  306  through  309  and  another  computer  program  is  also  developed  to 
numerically  integrate  equations  257  through  262.  These  computer  programs 
are  modified  here  to  include  equations  310  through  312,  thereby  effectively 
closing  the  loop  on  the  helicopter  drive  system.  These  modified  computer 
programs  are  listed  in  Appendix  VIII.  These  computer  simulations  are  used 
to  compare  the  stability  of  the  helicopter  drive  systems  with  respect  to 
both  time  response  and  frequency  response. 

The  time  response  of  each  system  is  determined  by  running  the  system  at 
steady  state  and  then  suddenly  stepping  the  set  speed  from  1447  radians 
per  second  to  1500  radians  per  second.  The  time  response  of  the  engine 
speed  is  then  observed  and  plotted.  It  is  nlso  determined  what  length  of 
time  is  required  for  the  system's  engine  speed  oscillations  to  settle  down 
to  within  1  rndian  per  second  of  its  final  value.  For  the  helicopter  drive 
system  without  a  ZTS  coupling,  anti  with  zero  shaft  damping,  the  time 
response  is  plotted  in  Figure  H8.  It  can  be  seen  that  this  system  is  diver¬ 
gent  and  is,  therefore,  completely  unstable.  Also,  when  operating  at 
steady  state,  any  slight  perturbation  in  the  system  will  cause  the  engine 
speed  to  diverge  from  its  steady-state  value.  This  is  in  complete  agree¬ 
ment  with  the  results  of  the  linearized  stability  analyses. 

Also,  iidr  the  helicopter  drive  system  with  a  ZTS  coupling,  and  with  zero 
shaft  damping  and  zero  coupling  damping,  any  slight  perturbation  from 
steady  state  causes  the  engine  speed  to  diverge.  In  addition,  in  attempt¬ 
ing  to  find  the  time  response  of  this  system,  the  coupling  numerically 
comes  apart.  Therefore,  this  system  is,  also,  completely  unstable.  This 
does  not  agree  with  the  linearized  stability  analyses,  which  shows  that 
some  accuracy  may  have  been  lost  in  the  linearization  of  the  ZTS  coupling. 
However,  if  a  coupling  damping  of  10  foot-pounds  per  radian  per  second  is 
introduced  to  this  system,  it  becomes  convergent  as  shown  in  Figure  88. 

Also,  from  Table  III,  it  takes  only  1.27  seconds  for  the  engine  spved  of 
this  system  to  settle  down  to  within  1  radian  per  second  of  its  final 
value.  Therefore,  this  system  appears  to  be  quite  stable. 

Next,  the  time  response  of  the  same  three  systems,  hut  with  a  shaft  damping 
of  2.24  foot-pounds  per  radian  per  second,  is  compared.  The  system  without 
a  ZTS  coupling,  whose  time  response  is  given  in  Figure  89,  converges,  but 
very  slowly.  From  Table  III,  it  takes  6.41  seconds  to  settle  down  to 
within  1  radian  per  second  of  its  final  value.  Therefore,  this  system  is 
stable,  but  it  is  not  quick  enough  to  be  practical.  The  time  response  of 
the  systen  with  a  ZTS  coupling,  but  with  zero  coupling  damping,  is  shown 
in  Figure  89.  Although  its  engine  speed  oscillations  are  rather  large  at 
first,  it  settles  down  quite  rapidly.  From  Table  III.  the  engine  speed 
takes  only  1.66  seconds  to  settle  down  to  within  1  radian  per  second  of 
its  final  value.  Therefore,  the  time  response  of  this  system  is  a  dis¬ 
tinct  improvement  over  that  of  the  system  without  a  ZTS  coupling,  despite 
the  large  initial  oscillations.  However,  the  time  response  of  this  system 
is  not  as  good  as  the  time  response  of  the  system  with  zero  shaft  damping 
and  with  a  ZTS  coupling  with  n  coupling  damping  of  10  foot-pounds  per 
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Figure  89.  Ti«e  Response  (C  =  2.24. 


TABLE  III.  ENGINE  SPEED  SETTLING  TIME 

ZTS  Coupling 

c  (  ft-lb  ) 

2  rad/sec 

c  (“;“■> 

c  rad/sec 

Time 

(sec) 

Wi thout 

2.  24 

0.0 

6.41 

With 

2.24 

0.0 

1.66 

With 

0.0 

10.0 

1.27 

With 

2.  24 

10.0 

1.05 

With 

2.  24 

20.0 

1.14 

With 

2.  24 

5.0 

1.21 

With 

2.  24 

6.0 

1.07 

With 

2.24 

16.0 

1.04 

With 

2.24 

17.0 

1.14 

With 

2.24 

25.0 

1.15 

radian  per  second.  Therefore,  It  appears  that  the  ZTS  coupling  acts  as 
a  means  for  artificially  Increasing  "shaft"  damping,  In  order  to  Improve 
the  closed-loop  system  stability. 


The  system  with  a  ZTS  coupling  and  shaft  damping,  and  with  a  coupling 
damping  of  10  foot-pounds  per  radian  per  second,  shows  an  even  more 
improved  time  response,  as  illustrated  In  Figure  90.  Also,  from  Table  III 
the  engine  speed  of  this  system  takes  only  1.05  seconds  to  settle  down  to 
within  1  radian  per  second  of  its  final  value.  Therefore,  this  system  Is 
much  more  stable  than  the  system  with  a  ZTS  coupling.  It  appears  that 
Increasing  the  coupling  damping  increases  the  stability  of  the  system. 
However,  this  Is  only  true  to  a  point.  This  is  demonstrated  by  determining 
the  time  response  of  the  above  system  with  an  increase  in  coupling  damping. 
A  coupling  damping  of  20  foot-pounds  per  radian  per  second  Is  used  and  the 
time  response  of  the  system  Is  shown  In  Figure  90.  From  Table  III,  the 
engine  speed  of  this  system  takes  1.14  seconds  to  settle  down  to  within 
1  radian  per  second  of  its  final  value.  Looking  at  Figure  90,  it  doesn't 
appear  that  doubling  the  coupling  damping  has  significantly  improved  the 
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tine  response,  and  it  has,  in  fact,  slightly  increased  the  settling  time. 
This  means  not  only  that  increasing  the  coupling  dnmping  beyond  a  certain 
value  docs  not  improve  stability,  but  also,  and  more  importantly,  as  long 
as  the  coupling  damping  is  some  minimum  value,  its  exact  value  is  rela¬ 
tively  unimportant.  This  fact  greatly  simplifies  the  problem  of  design¬ 
ing  and  constructing  the  coupling  damper.  In  order  to  determine  the  range 
of  coupling  damping  that  will  give  this  system  its  optimum  stability,  tin- 
time  response  of  the  system  is  found  for  several  different  values  of  cou¬ 
pling  damping.  It  is  found  that  the  time  response  of  the  system  is  vir¬ 
tually  unchanged  for  values  of  coupling  damping  from  fi  to  lf>  foot-pounds 
per  radian  per  second.  The  time  response  of  this  system  for  each  of  these 
values  for  coupling  dnmping  is  shown  in  Figure  91.  Also,  the  length  of 
the  time  required  for  the  engine  speed  of  this  system  to  settle  down  to 
within  1  radian  per  second  of  its  final  value,  for  several  different 
values  of  coupling  damping,  is  given  in  Table  III.  It  is  of  principal 
importance  that  the  coupling  damping  be  at  least  the  minimum  value  in  this 
range,  because  if  it  Is  below  this  range,  the  time  response  of  this  system 
gets  perceptively  worse.  However,  if  the  coupling  damping  exceeds  the 
maximum  value  of  this  range,  the  time  response  of  the  system  becomes  only 
slightly  worse.  Therefore,  it  is  necessary  to  design  the  ZTS  coupling 
with  a  coupling  damping  of  at  least  6  foot-pounds  per  radian  per  second, 
and  preferably,  not  more  than  1G  foot-pounds  per  radian  per  second,  in 
order  to  optimize  the  time  response  of  the  system. 

The  frequency  response  of  each  system  is  determined  by  varying  the  rotor 
damping  coefficient  sinusoidally,  as  would  occur  in  forward  flight,  at 
different  frequencies.  The  variation  in  engine  output  speed  is  plotted 
as  a  function  of  the  frequency  of  excitation.  The  frequency  response  is 
found  for  the  helicopter  drive  system  with  a  shaft  damping  of  2.2-1  foot¬ 
pounds  per  radian  per  second.  For  this  system  without  a  ZTS  coupling,  the 
frequency  response,  as  plotted  in  Figure  92,  displays  two  apparent  crit¬ 
ical  frequencies.  The  first  is  found  below  2  Hertz  which  is  below  the 
range  of  normal  system  excitation.  The  second  crlticnl  frequency  for 
this  system  is  at  about  7  Hertz.  System  excitation  at,  or  near,  this 
frequency  would  create  engine  speed  oscillations  that  could  be  destruc¬ 
tive.  Tills  critical  frequency  lies  dangerously  near  the  range  of  normal 
system  excitation,  which  is  about  8  to  10  Hertz. 

When  a  ZTS  coupling  is  inserted  in  this  system,  the  frequency  response 
changes  drastically,  as  il'ustratcd  In  Figure  92.  The  second,  and  more 
dangerous,  critical  frequency  is  completely  eliminated.  In  addition,  the 
engine  speed  variation  is  educed  for  all  excitation  frequencies  between 

2  and  10  Hertz.  There  Is  virtually  no  change  in  the  frequency  response  of 
this  system  with  a  ZTS  coupling,  if  coupling  damping  is  added,  for  excita¬ 
tion  frequencies  above  3  Hertz.  However,  for  excitation  frequencies  below 

3  Hertz,  the  existence  of  coupling  damping  tends  to  reduce  the  variation  of 
engine  speed.  Therefore,  it  appears  that  the  introduction  of  a  ZTS  cou¬ 
pling  greatly  Improves  the  frequency  response  of  this  helicopter  drive  train 
system,  and  that  adding  coupling  damping  improves  the  frequency  response 
even  a  little  more. 
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Figure  91.  Time  Response  (Acceptable  Extremes  of 


Kxcitotion  Frequency  au  (Hz) 

Figure  92.  Frequency  Response  (Variation  of 
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SUMMARY  OF  RESULTS 


The  helicopter  drive  system  analyzed  in  Reference  5,  with  a  direct  fuel 
flow  speed  governor,  is  not,  inherently,  a  very'  stable  system.  Although 
the  instability  is  primarily  due  to  the  speed  governor  lags,  the  amount 
of  shaft  damping  in  the  drive  train  has  a  critical  effect  on  the  stabil¬ 
ity  of  this  system. 

The  ZTS  coupling  designed  in  Part  I  acts  in  two  ways  to  improve  the  stabil¬ 
ity  of  this  system.  First,  it  provides  a  mode  for  artificially  increasing 
the  "shaft"  damping  of  the  drive  train,  in  order  to  improve  the  closed-loop 
stability  of  the  helicopter  drive  system.  This  is  done  by  using  an  elasto¬ 
meric  damper  as  a  torsional  damper  in  the  ZTS  coupling.  The  closed-loop 
stability  of  this  system  is  optimized  by  designing  the  coupling  damping  to 
be  in  the  range  of  6  to  16  foot-pounds  per  radian  per  second,  although 
letting  the  coupling  damping  be  greater  than  this  will  not  significantly 
hurt  the  stability  of  this  system. 

Second,  the  ZTS  coupling  acts  to  isolate  the  rotor  inertia  of  the  heli¬ 
copter  drive  system  from  the  engine  inertia.  That  is,  it  removes  rotor- 
induced  oscillations  from  the  system  before  they  reach  the  engine.  In 
this  respect,  the  ZTS  coupling  is  very  effective  in  that  it  completely 
eliminates  one  critical  frequency  from  the  system  and  reduces  the  ampli¬ 
tude  of  the  engine  speed  oscillations  over  the  entire  range  of  excitation 
frequencies  examined.  The  existence  of  coupling  damping  has  little  effect 
on  the  frequency  response  of  the  system. 

Therefore,  using  the  ZTS  coupling  designed  in  Part  I  with  the  coupling 
damping  in  the  range  discussed  above  is  an  effective  way  of  stabilizing 
the  helicopter  drive  system  presented  in  Reference  5,  with  a  direct  fuel 
flow  speed  governor.  In  this  way,  the  relative  simplicity  of  this  type 
of  speed  governor  does  not  have  to  be  sacrificed  in  order  to  stabilize 
the  system. 
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CONCLUSIONS 


The  principal  conclusions  of  this  investigation  and  analysis  are  as 
follows: 

1.  A  shaft  coupling  can  be  designed  to  meet  the  horsepower  and  speed 
requirements  of  modern  U.  S.  Army  helicopter  drive  trains  with 
zero  torsional  stiffness  at  the  design  speed  and  horsepower. 
Specifically,  two  different  types  of  zero  torsional  stiffness 
(ZTS)  couplings  were  designed  to  transmit  5270  horsepower  at 
13,820  rpm.  For  one  type  of  coupling,  simple  design  equations 
were  derived  to  size  a  coupling  for  any  helicopter  application. 

2.  Three  significant  design  innovations  for  ZTS  couplings  were 
developed  during  this  investigation: 

a.  The  counterbalanced  link  for  the  pinned-link  coupling  which 
allows  applications  at  much  higher  shaft  speeds  than 

were  heretofore  possible. 

b.  The  ZTS  coupling  damper  which  allows  the  introduction  of 
significant  amounts  of  shaft  damping  into  a  drive  train. 

c.  The  off-center  elastic  link  which  allows  operation  of  an 
elastic-link  ZTS  coupling  at  high  speed  without  the  negntive 
stiffness  exhibited  by  the  English  "llalf-iioon"  design. 

3.  A  ZTS  coupling  in  a  helicopter  drive  train  acts  as  a  torsional 
vibration  isolator  and  can  provide  significant  damping  which  is 
not  otherwise  available  to  the  designer.  Engine  and  control 
response  to  excitation  at  the  rotary-wing  can  be  reduced  or  com¬ 
pletely  eliminated  with  a  ZTS  coupling,  depending  on  the  ampli¬ 
tude  and  frequency  of  excitation.  As  an  example,  the  engine 
response  to  a  rotary-wing  excitation  of  ±  10r  full  load  torque 
was  found  to  be  almost  totally  eliminated  at  all  frequencies 
above  4  Hz  by  the  use  of  a  ZTS  coupling  in  the  drive  train  of 
Reference  5. 

4.  The  use  of  a  ZTS  coupling  with  damper  greatly  improves  the  tor¬ 
sional  stability  characteristics  of  a  helicopter  drive  train  with 
automatic  speed  control.  For  example,  the  drive  train  of  Refer¬ 
ence  5  with  a  simple  direct  fuel  flow  governor  and  sufficient 
gain  to  produce  5r  droop  was  found  to  have  negative  stability 
margins  (unstable)  with  no  ZTS  coupling,  and  large  positive 
stability  margins  (stable)  with  a  ZTS  coupling.  These  stability 
predictions  were  verified  by  numerical  simulations  of  closed- 
loop  drive  train  dynamics  on  a  digital  computer. 
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APPENDIX  I 

ELASTICITY  ANALYSIS  OF  CURVED  LIMC 


The  elasticity  analysis  of  the  elastic  link  shown  in  Figure  94  proceeds 
along  a  similar  vein  as  for  the  curved  bar  shown  in  Figure  93  which  has 
been  analyzed  by  Timoshenko  and  Goodier.*6  The  pin  force  P  at  A  and  B  is 
replaced  with  the  shear  force  shown  in  Figure  94  so  that  polar  coordinates 
can  be  used  to  specify  this  force  as  a  boundary  condition.  The  total  error 
caused  by  this  nssu^tion  is  small  as  a  consequence  of  the  principle  of 
Saint  Venant. 


Considering  the  external  force  p  only  and  treating  this  ns  a  plane  stress 
problem  (no  stress  in  Z  direction),  then 


ct=t  =t  -o 
7.  rz  z0 


(313) 


The  equations  of  equilibrium  and  compatibility  must  be  satisfied  in  con¬ 
junction  with  the  boundary  conditions  by  some  assumed  stress  distribution. 

It  has  been  shown  for  polar  coordinates  that  if  a  stress  function  $  is 
defined  such  that 


d20 

CTP  = 

Sr 


and 


A 

W 


Tre  -  5r 

the  equations  for  equilibrium 

Sa  ,  o  -o. 

r  1  r6  r  0  . 

■■y  ♦  —  ■  ■■  ♦  — — .  =  0 

r  °r  r 


1 


dT  2T 

r0  r0 


=  o 


are  satisfied  for  any  function  0(r,0). 


(314) 


(315) 


If 


0  =  f (r)  sin  0 


(316) 
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: ic  Link. 


? 


la  substituted  into  the  biharmonic  aquation  in  polar  coordinates 

then  compatibility  is  also  satisfied  by 
f(r)  =  Ar3  +  Bp+Cr+Drlogr 
It  aust  be  determined  if  all  the  boundary  conditions 


a  (r=a)  =  a  (r*b )  =  0 
r 


Tr0(r=a)  =  Tr9(r=b)  =  ° 


for  all  8 
for  all  6 


b  b 

f  Tpfl(9=TT)  dr  =  J*  -  rr0(0=O)dr  =  P 


can  be  satisfied  by  the  stress  function 
3  1 

0  =  (At  ♦  B  —  ♦  Cr  ♦  Dr  log  r)  sin  0 

T 

Substituting  0  into  the  stress  equations  (314)  and  applying  the  B.C 
<j(r=a)  =  0  gives 
2B  D 

2Aa  -  —ir  +  —  =  0 

7  ■ 

a(r=b)  s  0  gives 


2B  D 


Mb  '  ~S*  E=  0 


w 

f  Tr0<fts°)  =  -P  gives 


[ao.2-.2)  *  Olo*r]  * 


Solving  equations  320  through  322  for  A,  B,  and  D  and  letting 


(317) 


(318) 


(319) 


(320) 


(321) 


(322) 


199 


2  2  2  2 
N  =  a  -b  +  (a  +b)  log  b/a 

gives 

A  =  -P  2N 

B  =  (a2b2P  2N)  (323) 

2  2 
„  P(a  +b  ) 


The  term  coefficient  C  in  equation  (317)  is  not  needed  to  satisfy  the  U. C. 
The  stresses  are  now  found  from  equations  31  *1  to  be 


r-  2  2  2  2  -i 

■pi  r  *  2-|_  -  la  *b  il  si„  e 

r3  r  J 


r-  2  2  2  2  -f 

-P/K  [3r  -  2-5-  -  sl„  8 


2  2  2  .2,-, 

T  =  P  'N  j  r  + 

re  L  rJ  r  J 


'  N  Tr  +  - 


cos  8 


(324) 


which  are  the  same  as  the  Timoshenko  and  Ooodier’s  analysis  for  the  curved 
bar  in  Figure  93. 

To  find  the  displacements  u  and  v  in  the  radial  and  tangential  directions 
respectively,  the  strain-displacement  relationships  325  are  used. 

C  =  ^u/*r 

r 

c0  =  u/r  +  1  r(dv  d8)  (325) 

v  .  -  l/r(du  ’^0)  *  *v  *r  -  v  r 
rt 

In  conjunction  with  Hooke's  law, 

€  =  1/E (C  -  uo.) 

r  r  9 

€g  =  l/E(Cg  -  MCr>  (326) 

vr6  =  1  G(~  » ) 
r8 
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where  E  is  the  Modulus  of  elasticity,  G  is  the  modulus  of  rigidity,  and 
li  is  Poisson's  ratio  to  give 


du/dr  -  l/E[2Ar  -  2B/r3  +  D  r  -  uGAr  +  2B/r3  +  D  r)  1  sin  0  (327) 

Integration  of  327  gives 

u  =  sin  0/E[ (l-3u)Ar2  +  (l+u)B/r2  +  (l-u)D  r]  4  f(8)  (328) 

From  the  second  of  equations  325, 
dv'de  =  refl  -  u 

and  substituting  e  from  equation  326  and  u  from  327  gives 

0 

dv'30  =  sin  6  E[(5+u)Ar2  +  (l+p)R  r2  +  (l-u)D 

-  (l-u)D  i-  r]  -  f(0)  (329) 

Integration  of  the  above  gives 

2  2 

v  =  -  cos  0  f:((5+u)Ar  +  (l+u)B  r  +  (l-^i)D 

-  (l-u)D  1-  r]  -  /  f (9)d8  ♦  g(r) 

From  the  last  of  equations  325  and  326, 

\e  --  rllr  **  *r  -  v  r  *  TrC  0 

Replacing  G  with  E/2(lip),  differentiating  u  with  respect  to  8  and  v 
with  respect  to  r  and  then  substituting  into  equntion  331  gives,  after 
some  rearrangement  and  canceling  of  terms, 

-4D  cos  8  F  =  f(0)  +  rg(r)  -  g(r)  4  j  f(8)d6 

or 

f'  (8)  4  f  f  (0)d8  4  (41)  cos  5  E)  =  rg  (r)  -  g(r)  (332) 

The  left  side  of  equation  332  is  a  function  of  8  only  and  the  right  is 
a  function  of  r  only.  A  solution  to  this  equation  can  be  found  if  both 
sides  of  equation  332  are  set  equal  to  zero  which  gives 


(330) 


(331) 
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^  k  •« 

or 

g(r)  =  Hr  (333) 

Differentiating  the  left  side  of  equation  332  with  respect  to  6  gives 

f  (8)  +  f (8)  =  4D  sin  8/E  (334) 


with  homogeneous  solution 

f(8)  =  K  sin  8  L  cos  8  (335) 

H 

To  find  the  particular  solution  of  334,  let  f(8)p  =  M8  cos  8,  then 
f/(9)  =  M  cos  8-  M8  sin  8 
f  (8)  =  -M  sin  8  -  M  sin  8  -  M  cos 

and 

f (6)  +  t" (8)  =  -2M  sin  8  =  4D  sin  8/e 
or 

U  =  -2D/E 

f (8)  =  (-2D/E)  8  cos  8 

P 

f (8)  =  (-2D/E)8  cos  8  +  k  sin  9  +  L  cos  8  (336) 


The  Integral  of  f(8)  is  needed  in  equation  330  and  is  given  in  equa¬ 
tion  337  after  integrating  336  by  parts. 

/f(8)d8  =  J(-2D/E) 8  cos  8  -  k  cos  8  +  l  sin  8 

p  =  8  dp  =  d8 

dq  =  cos  8  d8  q  =  sin  8 

Jf(8)d8  =  (-2D/E) f 8  sin  8  -  J*  sin  0  d8]  -  K  cos  8  +  L  sin  8 

Jf(8)d8  =  ( -2D/E) ( 8  sin  8  +  cos  8)  -  K  cos  8  +  L  sin  8  (337) 

After  substitution  of  equations  333,  336,  and  337  into  equations  328  and 
330,  the  expressions  for  u  and  v  become 
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(338) 


u 


v 


a  *111  0/E[  (l-3p)Ar2  ♦  (l+p)H/r2  +  (l-n)D  On  r) 

-  (2D@  cos  0/E)  +  K  sin  0  +  L  cos  0 
=  -cos  0/E[(3+p)Ar2  +  1  +  pB/r2  +  (l-P)D-  (1-WD  in  r] 


-2D(0  sin  9  r  cos  9) 
E 


K  cos  0  +  L  sin  0  +  Hr 


(339) 


The  arbitrary  constants  Hf  K,  and  L  oust  be  evaluated  froa  the  boundary 
conditions  for  displacement,  which  are: 

u(r,o)  s  0 

v(R,0)  a  0  where  R  a  (afb)/2  (340) 

v(R,tt)  =  0 

The  first  E^Cq  when  applied  to  equation  338  gives  L  =  0,  the  second  B^Cq 
applied  to  equation  339  results  in  the  following  relationship  between 
K  and  H: 

K  a  HR  ♦  l/E(D(p-3)  -  (i»5)AR2  -  (l+p)B/R2  ♦  (l-p)D  Sn  R]  (341) 

Finally,  applying  the  last  BqCo  to  equation  339  gives 

K  =  -HR  -  l/E[-D(p-3)  +  (p*5)AR2  4  (l+u)B/R2  -  (l-p)D  6n  r]  (342) 

Adding  341  to  342  and  rearranging  terms  gives 
H  =  0 

K  =  (^"2)D  -  1/E(  (54-u)AR2  ♦  (l+n)B/R2-  d-u)D  U  R] 

Substituting  the  values  of  A,  B,  D,  H,  K,  and  L  into  xiuatlon  338  along 
with  the  coordinates  of  the  position  where  the  deflection  is  desired, 
u(R,tt),  yields 


u(R,tt)  =  2ttD/E 


-2TTP(a2+b2) 

3  2  2  3 

E[  (b  -a  )  +  (a  +b  )  log  a^b 


(343) 


Equation  343  gives  for  displacement  A  between  points  A  and  B,  due  to  F: 
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A  = 


2rrP(a^+b2) 


E[a2-b2)  +  (a2+b2)  log  b./a] 


the  spring  constant  K  is  P/A;  thus, 

P 


Kp  = 


I((a2-b2)  +  (a2+b2)  fc  b/a] 

2  2 
2tt(b  +b  ) 


(unit  thickness) 


The  deflection  A  is  twice  that  of  curved  bar  problem  solved  by  Timoshenko 
and  Goodier.  This  result  might  be  expected  if  one  considers  the  curved 
bar  problem  in  Figure  93  as  simply  half  of  the  symmetric  elastic  link. 

The  analysis  presented  here  with  the  details  of  each  step  should  suggest 
a  method  of  solution  for  other  curved  link  problems. 
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APPENDIX  II 

REPRESENTATION  OF  SPEED  REDUCERS 


To  simplify  the  dynamic  analysis  of  systems  with  speed  reducers,  it  is 
common  practice  to  reference  all  spring  constants  and  inertias  to  one 
side  oi  the  speed  reducer.  Figure  95a  shows  a  three-degree-of -freedom 
system  which  has  the  following  relationships  between  the  dependent  coor¬ 
dinates  ©2  and  03: 

03  =  N02  (314) 

03  =  N§2  (345) 


where  N  is  the  gear  ratio. 

If  the  coordinate  6^  is  referenced  to  the  left  side  of  the  speed  reducer, 
then  the  referenced  position  0,  and  the  referenced  velocity  6'  are 

4  *i 

0'  --  6^  S  (346) 

0'  =  0 .  /N  (347) 

4  4 


The  kinetic  energy  of  the  system  is 

KE  =  1/21,9*  .  1  2I29*  .  1  2I0  8*  .  1  2I0  9* 


or  upon  substitution  of  equations  345  and  347, 

KE  =  V'2^02  ♦  1'2I9N2(01)2  *  1  2 ( I c  +  lc  N2)02 


The  potential  energy  of  the  system  is 

PE  =  1,  2K1(0O-01)2  +  1  2K0N2(e.>-e'J) 

Letting 

/  2 
K2  =  K2N 


X2  &  X2N 
lG  "  lGl  + 


tne  three-degree-of -freedom  system  takes  the  form  shown  in  Figure  95b. 
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Figure  95.  Representation  of  Speed  Reducer. 
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APPENDIX  III 

FORMULATION  OF  LINK  VELOCITIES 

An  expression  for  the  angular  velocity  of  each  link  as  well  as  the  linear 
velocity  of  the  center  of  gravity  of  each  link  is  determined  by  using 
velocity  polygons. 

The  velocities  shown  in  Figure  96  are 

V .  =  velocity  of  eg  of  link  1 

X1 

V-  =  velocity  of  eg  of  link  2 
2 

and 

V. ,  V, /,  and  B  =  velocities  of  points  A,  a' ,  and  B. 

A  A 

A  velocity  in  the  form  V  is  the  velocity  of  point  i  relative  to  point  J. 

*  /  J 

Velocities  with  subscripts  x  and  y  are  component  velocities  in  the  x  and 
y  directions. 

The  equations  representing  the  vectors  in  Figure  96  are 


V 

VA  +  Vf/A 

(348) 

II 

K 

H 

> 

VA  *  Vi/A 

X  X 

(343) 

vV 

VA  *  VX/A 

y  y 

(350) 

Vi/A 

1/2vb/a 

(351) 

VB/A  = 

VB  -  VA 

(352) 

VB/A 

X 

=  VD  -  V. 

B  A 

X  X 

(353) 

VB/A 

=  V„  -  V. 

B  A 

(354) 

y  y  y 


The  magnitude  of  from  equations  349  and  350  is 


(355) 
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Using  equations  351,  352,  and  353,  equation  355  takes  the  following  form: 


4 


(356) 


The  expression  for  the  magnitude  of  the  velocity  of  link  2  Is  similar  to 
equation  356  with  A  replaced  with  a'  giving 


(v. 


+  v. 


)2  + 


(V, 


>' 


(357) 


For  point  B, 


vn  =  R( 8  -  0/2) 

B  © 

y 

thus 

9  •  2  2  •  •  ? 

v“  =  R  ♦  R  (0  -  0/2) 

B  e 


For  point  A, 

v.  -  -r0  sin  0/2 
A  e 

x 

v.  =  rfi  cos  9/2 
A  e 

y 

For  point  a' , 

v  /  =  r(8  -9)  sin  S  2 
A  e 

x 

v .  /  --  r (9  -9)  cos  9  2 

A  e 

y 

2  2  2 

Let  v i  ~  v i  v i  (see  equation  54) 


(358) 

(359) 

(360) 

(361) 

(362) 

(363) 

(364) 

(365) 


Substituting  equations  358  through  364  into  equations  356  and  357  gives 
the  following  equation  for  v^  in  terms  of  R,  R,  9,  and  0. 
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v?  =  l/4{(-r0  sin  0/2  ♦  R)2  +  r0  cos  0/2  +  R(0  -B/2)}2 
JL  e  e  e 

+  [r(0  -0)  sin  0/2  +  R]2 
e 

+  [r<0  -0)  cos  0/2  +  R(0  -0/2) ]2 
e  e 

Squaring  and  combining  terms  in  equation  366  gives 

v2  =  l/2f r2(02  -  0  6  +  02/2)  -  rR0  sin  0/2 
i  e  e 


♦  R2  +  2rR(0  -  0/2)2  cos  0/2 
e 


2  •  •  2 
♦  R  (0  -0'2)  } 
e  ‘ 


The  angular  velocity  of  link  1  is 
v. 


B/A 


or 

2  1  r,  . 2 

(JU.  =  -r  =  r  (v  -  V 

X-  B  A 

IX  xx 

and  for  link  2 


)2  .  (v,  -,A  ,*) 


ID , 


vb/a' 


where  V„,./  =  V  -  V./ 
B/A  B  A 


or 


<  =  i/i2r('B  -vA/)2  4  (V  -vAo2} 

2  xx  y  y 


2  2  2 
Let  ui .  =  ai  t  +  it, 

*  *1  2 


(see  equation  54) 


Substituting  equations  358  through  364  into  368  and  369  gives 

i)2.  =■  2/X2fr2(02  -  0  0+  92,  2)  +  rR0  sin  0/2 
X  e  e 

+  R2  -  2rR(0  -  0/  2) 2  cos  Q'2  +  R2(0  -  6/2)2} 
e  e  * 


(366) 


(367) 


(368) 


(369) 


(370) 
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appendix  IV 
COMPUTER  PROGRAMS 

WOK  COUPLING  DESIGN 

Design  Equation  I — *(P.6> 


Input  Data 

Quantity 

Value 

Divans ions 

N 

2 

none 

■ 

.2385 

lbm 

r 

3. 

in. 

8 

.87 

none 

O) 

13820 

rp* 
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Stress  and  Force  Analysis 


Input  Dats 


Quantity 

Value 

Dimensions 

■ 

.2385 

lb 

m 

"i 

.130 

lb 

m 

r 

e 

1.3 

in. 

L 

2.61 

in. 

X 

1447 

rad / sec 

r 

3.0 

in. 

A1 

.375 

in2 

A2 

.6875 

in2 

A3 

to 

cn 

in2 

A4 

.4418 
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I/C 

.143 

in3 

"b 

.416 

lb 

m 

*3 

.416 

lb 

m 

213 


• 

4* 

03 

04 

U. 

• 

• 

• 

» 

o 

** 

o 

< 

«-l 

X 

•  H- 

a. 

.  *. 

mi 

z  - 

< 

M  <t 

mm 

•  u. 

z 

»  • 

m* 

<o  * 

to 

fO  ■£> 

ft 

>-  O' 

>»• 

•>  h- 

LL 

•  « 

1 

<  " 

to 

</>  m 

ft 

f0 

•  u. 

«« 

X 

•>  • 

04 

» 

<0  • 

U. 

CO  fO 

<N  O 

♦ 

X  X 

»-  CO 

•-4 

»  *. 

•>  h- 

u. 

o  o 

C'  » 

«4 

H  X 

«>•  « 

mm 

«  » 

a  < 

♦ 

<4-  O 

<  LL 

04 

<  •— 

Of  - 

ft 

w  »> 

ft 

x  'O 

mm  mm 

«n  <  < 

►-  0- 

<  — 

X  -  • 

•  ►“ 

X  < 

»N« 

*  • 

a.  X 

C3  <  < 

<o  • 

-i  Cl. 

X  -  - 

—  X 

**  mi 

•  m*  Ol 

*-  o 

—  < 

o  <  < 

•  X 

X 

»■* 

to  — 

X  »  • 

m  m 

• 

O  to 

*■  of  -1 

mm  m 

3 

• 

u  o 

UJ  •  < 

o  o 

• 

to 

—  xo 

X  3  « 

LU  <© 

of 

• 

—  3  • 

•  «  Of 

a  *- 

• 

of 

<  *  <r 

X  J  *■ 

mm  m 

• 

X  3  U  U. 

O  -3 

X  - 

• 

of 

CL  ft  •  ♦ 

L  111  • 

*-  CD 

04 

1 

«■» 

mi  —  Of  — 

»  Of  -1 

• 

-4 

-4 

<  <  Of  04 

X  O  *  i 

m  ™ 

^  mi 

• 

>» 

—  X  •  u. 

<  •*  JUJ  - 

CO  • 

U4  _l 

a 

X 

to  &  *» 

of  -  x  a'  o  »-  'O  of  X  w  Ojj-- 

O  j  •  •  •  -  irv  r-  o-  r»  *  ♦  3  o  •  i  v-  z  3  o  <  X  a. 

o  c  r-i  •-*  x  <m  uj  •  co  n  m  a;  »-  •  •  ~  *  w 

.mf»  .  .  .  .  t  —  x  3  h  n  s  u  to  3  •  z  a  — 

(LZh  »u.  •»  »  h^nnnn  a  s  —  •  s  o  of  •  (M  —  X  -. 

•  —  oOxj-  •  co  m  -4  •  •  x  -J  ♦  of  k  ♦  <  of  *—  to  —  l- 

l/l  I  Ifl  ••-''V-l/lHfv/NSSSMI'J  ♦  •  u  •  •  0.  — '  —  l_>  K  Of  S  ■  •  Of 

*/>  —  ►-  <  oj  <  -  ««  jawojHrt  jujrt  tH«Zi/)«iM«itn^ir>u 

1U  JDZlhX  •nslSj;o£\S>.  X-i.O*h"U«  •  X  X  LL  U.  •  lO 

af<<«of— of<o  II  x  u  h  n  u  _j  af  «  n  n  n  /  n  i/)  O  n  a.  n  n  ii  » 

HUIW0£O0£O«f**lljartUJ  il  II  UJ  O  H  IN  N  N  of  M  |l  Of  N  J  ^  3  < 

Maocau.ju.i'a.xxixajKiiu.u.hi-i-^uaziu.ku.L.k 

-•  ^  04 


u 


214 


K0H*KE«RR*M»M*H«12.*«74/2 
$BEAR*FA/A2 
SL02»*0M/ IC+F5/AI 
SB*FA/A3 


a 

u. 

• 

■4- 

U. 


< 

U. 


O 

X 


•  o 

co  • 

i/l  o 

*> 

«M  U. 

Q  O' 

_l  - 

t/1  «*» 

•  • 

0 C  O 
<HN 

iu  te¬ 
as  •  o 
*A  O  *- 
•  • 

OC  O  O 
H-  O 

•  u.  — 

H  »  lA 
_  •  M 
A  H 
^  »  •  • 

<  o  ~  ^ 

^  - - » 

CD  41  <  • 

U.  t-  X  >-  O. 

n  «—  at  —  O  O 
fifiC  OU.H  Z 
^  2  u.  -  vuu 

m 


215 


Simulation  of  System  1 


Input  Dot  a 


Quantity 

Value 

Dimensions 

Count  (inte¬ 
gration  steps 

1 ine  of  output ) 

10 

none 

I 

u 

.  709 

slug-ft2 

!r 

3.  888 

slug-f t“ 

c 

s 

2.32 

f  t-1  l>-sec 

K 

s 

212. 1-1696.  8 

ft -lb  rad 

C 

e 

. 000382 

ft-lb-scc2 

CL 

.  000956-1 

f t-lb-sec2 

T 

e 

2801.6 

rt-n> 

Start  time 

Variable 

see 

End  time 

Variable 

see 

i 

Variable 

rad  sec 
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SIMULATION  OF  SYSTEM 
IMPLICIT  RL-AL -£ ( A-Z  I 
I  N'T  CCER  f..N,  Nf.N,  COUNT 
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Simulation  of  Systea  2 


Input  Data 

Quality 

Value 

Dimensions 

N 

2 

none 

Count  (integration  steps/ 

line  of  output 

10 

none 

r 

3 

in. 

L 

2.61 

in. 

• 

m 

B 

.  2385 

lb 

m 

* 

0 

lb 

m 

ll 

.0003 

slug-ft2 

CD 

0 

ft-lb-sec 

TT 

(reverse  stop  contact) 

1 

rad 

kl 

(reverse  stop) 

66,800 

ft-lb/rad 

TTT 

(forward  stop  contact) 

1.97 

rad 

kl 

(forward  stop) 

66,800 

ft-lb/rad 

RFL 

** 

(free  length  of  spring  ) 

0 

in. 

KE 

** 

(linear  spring  constant  ) 

0 

lb  'in. 

C 

8 

2.32 

ft-lb-sec 

K. 

212.1-1696.8 

ft-lb/rad 

Any  cosfolnatlon  of  ■  and 

B 

which  sums 

to  m  can  be  used, 

**A  coupling  spring  can  be  added  between  the  link  pins  (0)  and  the 
center  of  the  coupling  (not  used  In  analysis). 
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Quality 

Valua 

Dimensions 

I 

• 

.5 

aluf-ft2 

*R 

3.888 

aluf-ft2 

*fe 

.015 

slu*-ft2 

‘fR 

.209 

alu*-ft2 

C 

.000382 

ft-lb-aec2 

e 

CL 

. 0009564 

ft-lb-aec2 

T 

2804 

ft-lb 

e 

Start  Tiae 

Variabla 

sec 

End  Tlae 

Variable 

sec 

U) 

Variable 

rad/sec 
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SIMULATION  OF  SYSTEM  2 
IMPLICIT  REAL*8(A-Z) 
INTEGER  NNtNNNfCOUNTiNM 
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FUNCTION  TLARGE<T,TT,TTT,X»Y| 
IMPLICIT  R£AL*8 ( A-Z  ) 
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APPENDIX  V 

DERIVATION  OF  DRIVE  TRAIN  TRANSFER  FUNCTION 


The  following  Is  the  derivation  of  the  drive  train  transfer  function  used 
in  References  1,  6,  and  9,  and  given  by  equation  209.  This  derivation  wa3 
adapted,  in  part,  from  a  derivation  by  R.  M.  Swick  and  C.  A.  Skarvan. 


Equations  144  through  149  are  substituted  into  equation  150,  and  the  shaft 

damping  coefficient,  C  ,  is  set  equal  to  zero,  to  give  equation  371. 

s 
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AT  (s) 
e 


rSWRl  „2  rK?(IR+Ieg)  +  CDCfl 

+  Lk0(cd+ct)  J  +  L  vw  J 


S+l 


(371) 


The  denominator  of  equation  371  is  factored  by  long  division  as  follows: 

.2 


[Veg  1  2  f^eg^^R 
LK9(IR+Ietf)J  -I„)2  - 


S+l 


9  eg  R 
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(Vcd> 

LwvJ 

Lvw  J  L 

Ke<V  V  J  + 

f1*1**  1  s3  , 

I'**'*  1  s2 

|k«(C  +C  )  j  s 

Ik. (i  +  i„ ) 

L  0  T  D  J 

L  0  eg  R  J 

n  ^Wr  1,2  iVWV'Wl 

+  LWV'VVJ  +LK9(VCD>  J 

2  2  2  2 
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K£(CVCn)(I  *ID> 

6  T  D  eg  R 


kQ(I  +I„)3  *  2C_C  I  I  -  c2,2  -C“l~ 
6  eg  R  D  T  eg  R  D  eg  T  R 


Ke(CT+CD)(Ieg+IR)' 


.  VW3  -  'Vc  -  ct'ii>: 

WV'VV2 


Under  the  assumption  that 


u  eg  r  R  o  eg  R 


the  above  equation  reduces  to 
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For  th.  values  of  the  p»rai«,tcr  uaoO  la  .1.1.  report,  .ho  !«..  »•■*« 

of  the  inequality  becomes 


f  (2.768)  (0. 69'l  i)  -  (1 . 106)  (3. 888)  1  -  5.65 


while  the  right  side  becomes 


2 

#00 . 0(-0. 6944 f  3. 88)  =  7.68  >  10 


Therefore,  the  inequality  is  true  and  the  assumption  based  on  it  is  a 
reasonable  one. 


After  the  denominator  is  factored,  the  drive  train  transfer  iunction  is 
given  by  equation  372. 


i.W  (s) 
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(<VCJ  [*8  S  '  K8  '  *- 


AT  (s) 
e 


7v»v  .i  T) s2. r-” * Ct'^) s. .I 

LvvcJ'  J  LlvvvJ  v  J 


(372) 


The  following  substitutions  are  made,  leaving  the  drive  train  transfer 
function  in  its  final  form,  given  by  equation  373. 
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APPENDIX  VI 

BODE  ANALYSIS  COMPUTER  PROGRAM 

The  plotting  subroutines  used  in  this  program  are  available  on  the 
FORTRAN  IV  C  Compiler  at  the  Northeast  Florida  Regional  Data  Center. 

Input  Data 


Quantity 

Value 

(and  Meaning) 

Units 

SSD 

5.0 

percent 

T 

B 

0.1869 

seconds 

T1 

0.02 

seconds 

T2 

0.03 

seconds 

T 

e 

0.3 

seconds 

T 

I 

0.773 

seconds 

T 

r 

0.2 

seconds 

kb 

0.941 

ft-lb 

Ib/hr 

Ke 

800.0 

ft-lb 

rad 

!r 

3.  888 

ft-lb 

.  .  2 
rad/sec 

I  or 
e 

(with 

ZTS) 

0. 5  or 

ft-lb 

*eg 

(without 

ZTS) 

0.6944 

rad/sec2 

CD 

2.768 

ft-lb 

rnd/sec 

CT 

1.106 

ft-lb 
rad  sec 

kr 

0.  5 

none 

C 

a 

variable 

ft-lb 
rad  sec 

234 


variable 


0.015 

0.2094 


rad/aec 

ft-lb 

2 

rad/sec 

ft-lb 

— — 5 

rad/aec 


1  =  CDP/LGR  speed  governor 

2  =  FE/LGR  speed  governor  none 

3  =  DFF  speed  governor 

1  =  Standard  Drivetrain  none 

2  =  Decoupled  Rotor  Model  of  ZTS  Coupling 

3  =  Linearized  Design  of  ZTS  Coupling 


1  =  High  Power  none 

2  =  Low  Power 

1  =  use  Same  Drivetrain 

Parameters  Done 

2  =  Read  New  Drivetrain 
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APPENDIX  VII 

LINEAR  TIME  RESPONSE  COMPUTER  PROGRAM 


The  subroutine  DPRQD  Is  available  In  the  IBM  Scientific  Subroutine  Package. 
The  plotting  subroutines  used  in  this  program  are  the  same  as  those  listed 
in  Appendix  VI.  The  units  for  all  input  data  in  this  program  may  be  found 
in  the  data  list  in  Appendix  VI. 
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APPENDIX  VIII 

COMPUTER  SIMULATION  OF  NONLINEAR  SYSTEMS 

Any  input  data  whoa*  dimensions  art  omitted  b*lov  may  b*  found  in 
Appendix  VI. 

I .  Simulation  of  System  Without  ZT8  Coupling 

Input  Data 


Quantity 

Value  (or  Meaning) 

Unlta 

Count  (integration 
steps  per  line  of 
output) 

10 

none 

‘eg 

0.6944 

- 

!r 

3.888 

- 

Ke 
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- 

c. 

2.24 

- 
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0.000382 
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* 

C, 

0.0009584 
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L 

Start  time 
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seconds 

End  time 

variable 

seconds 

T1 

0.02 

- 

T2 

0.03 

- 

KP 

40.27 

- 

&T 

'S? 

f 

-3. 3333 

ft-lb 

rad/sec 

dr 

0. 3523 

ft-lb 

lb/hr 
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Quantity 

Valua  (or  Meanlnc) 
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dT 

e 

W 
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0.  37 
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lfeAr 

dT 

• 

ST 

c 

3.2563 

ft-lb 
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KEYl 

0  for  frequency  raaponaa 

1  for  tine  raaponaa 

nona 

01 

variable 

rad/aac 
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Start  time 

variable 

seconds 

End  time 

variable 

seconds 

T1 

0.02 

- 

T2 

0.03 

- 

KP 

40.27 

- 

dT 

ST* 

g 

-3.3333 

ft-lb 

rad/sec 

dT 

0. 3523 

ft-lb 

lb/hr 

3t 

0. 57  • 

ft-lb 

lb/hr 

dr 

e 

3.  2563 

ft-lb 

ui 

g 

rad/sec 

KEY1 

0  for  frequency  response 

1  for  time  response 

none 

<u 

variable 

rad/sec 
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LIST  OF  SYMBOLS 


T 

E(s) 

hp 

I 


I 


eg 


gb 


mechanical  advantage  of  set  speed  control 

cross-sect  tonal  area  of  speed  governor  fuel  regulating  piston 

viscous  ZTS  coupling  damping  coefficient 

linearized  aerodynamic  rotor  damping  coefficient 

aerodynamic  power  turbine  damping  coefficient 

speed  governor  fuel  regulating  piston  viscous  damping 
coef  f iclent 

aerodynamic  rotor  damping  coefficient 
viscous  shaft  damping  coefficient 

linearized  aerodynamic  power  turbine  damping  incff tcient 
speed  error  Laplace  function 
shaft  horsepower 

moss  polar  moment  of  Inertia  of  power  turbine 

mass  polnr  moment  of  inertia  of  power  turbine  and  transmission 

mass  polar  moment  of  inertia  of  gas  producer 

mass  polar  moment  of  inertia  of  transmission 

mass  polnr  moment  of  Inertia  of  a  ZTS  coupling  link  about 
its  center  of  gravity 

mass  polnr  moment  of  i  ertia  of  helicopter  rotor 
engine  transfer  function  gain  torque  to  fuel  flow- 
speed  governor  piston  spring  stiffness 
helicopter  load  gain 

governor  proportional  gain,  fuel  flow  to  speed 

power  turbine  governor  reset  gain,  speed  request  to 
actual  speed 
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LIST  OF  SYMBOLS  (Continued) 

governor  proportional  gain,  fuel  flow  to  speed 
engine  torquemeter  shaft  stiffness 
helicopter  mast  shaft  stiffness 
leng*  of  ZTS  coupling  link 

concentrated  mass  at  end  of  ZTS  coupling  link 
mass  of  speed  governor  fuel  regulatinr  piston 
distributed  mass  of  ZTS  coupling  link 
compressor  discharge  pressure 

pressure  proportional  to  engine  power  turbine  speed 
distance  from  center  of  ZTS  coupling  to  interior  link  pin 
Laplace  operator 
steady-state  speed  droop 
torque 

torque  across  ZTS  coupling 

engine  torque  applied  to  power  turbine 

Laplace  function  for  engine  torque  applied  to  power  turbine 

torque  applied  to  gas  producer 

load  torque  on  turbine  engine 

power  turbine  drag  torque 

Laplace  function  for  powor  turbine  speed 

fuel  flow 

Laplace  function  for  fuel  flow 

Laplace  function  for  gas  producer  speed 

speed  governor  fuel  regulating  piston  displacement 

power  turbine  speed  error 
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1,1  ST  OF  SYMBOLS  (Continued) 
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9(s) 
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r.  <  s  > 
2>i  (s) 
f(t) 
f  (t) 
f(t) 


helicopter  rotor  damping  r:itm 
power  turbine  damping  r.itio 
/TS  coupling  angle 

Laplace  function  for  /,TS  coupling  angle- 
power  turbine  angle  of  rotation 
helicopter  rotor  angle  of  rotation 
angle  of  twist  in  shaft 

Laplace  function  for  angle  of  twist  in  shaft 
engine  gas  producer  lead  time  constant 
engine  gas  producer  lag  time  constant 
load  lag  time  constant 

power  turbine  governor  reset  lag  time  constant 
speed  governor  lag  time  constant 
speed  governor  lag  time  constant 
power  turbine  speed 
gas  producer  speed 

undamped  natural  frequency  of  po-.er  turbine  and  rotor  on 
flexible  shaft 

natural  frequency  of  rotor  drive  system 
helicopter  rotor  speed 
speed  governor  set  speed 

Laplace  function  for  speed  governor  set  speed 
incremental  change  in  a  Laplace  function 
incremental  change  in  a  time  domain  function 
steady-state  value  of  a  time  domain  function 
time  derivative  of  a  time  domain  function 
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